Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,871,605 C→A V324F (GTC→TTC)  SAUSA300_RS09270 ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,871,6050CA98.5% 398.5 / ‑3.3 136V324F (GTC→TTC) SAUSA300_RS09270diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (1/132);  total (1/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 136 total aligned reads displayed.

TCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGACATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCATTATTATCA  >  NC_007793/1871508‑1871699
                                                                                                 |                                                                                              
     tctttCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTa                                                                                         <  1:1967460/97‑1 (MQ=255)
     tCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTa                                                                                         <  1:1976248/100‑1 (MQ=255)
     ccgatCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTa                                                                                         <  1:1285563/96‑1 (MQ=255)
      cgatCTGGGTGGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGGCACGTCACCCATGAAATCTGTAt                                                                                        <  1:4906800/97‑1 (MQ=255)
      ccTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTag                                                                                        <  1:4652842/100‑2 (MQ=255)
      ccTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTAt                                                                                        <  1:5540178/100‑1 (MQ=255)
      ccTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTAt                                                                                        <  1:1710694/100‑1 (MQ=255)
       cTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGtgct                                                                                       <  1:4207722/100‑4 (MQ=255)
        ttCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATgg                                                                                      <  1:4344266/100‑1 (MQ=255)
        ttCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATgg                                                                                      <  1:648209/100‑1 (MQ=255)
        ttCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATgg                                                                                      <  1:1000950/100‑1 (MQ=255)
         tCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATggt                                                                                     <  1:3612033/100‑2 (MQ=255)
         tCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGa                                                                                     <  1:211741/100‑1 (MQ=255)
           tGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGaag                                                                                   <  1:2924830/100‑2 (MQ=255)
           tGGGTAACGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAAt                                                                                   <  1:1539359/100‑1 (MQ=255)
                        ttGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCt                                                                      <  1:4375682/100‑1 (MQ=255)
                          gATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACCCCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCttt                                                                    <  1:2665347/100‑1 (MQ=255)
                             aCCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCTCCCATGAAATCTGTATGGAATTTCATTCCCATCttttat                                                                 <  1:635815/100‑3 (MQ=255)
                             aCCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCttttat                                                                 <  1:4077182/100‑3 (MQ=255)
                               tcttccgatcTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:2781594/91‑1 (MQ=255)
                               tcttccgatcTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:2616899/91‑1 (MQ=255)
                               tcttccgatcTGCGTTCTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:753666/91‑1 (MQ=255)
                               tccgatctacTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:4472576/91‑1 (MQ=255)
                               cgatctggccTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:432631/92‑1 (MQ=255)
                               ccttccgatcTGCGTCTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:1407653/91‑1 (MQ=255)
                               cGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCGCGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:4507101/100‑1 (MQ=255)
                               cGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:2574105/100‑1 (MQ=255)
                               cGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:1531080/100‑1 (MQ=255)
                               cGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCAcc                                                               <  1:425356/100‑1 (MQ=255)
                                gatctgACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:681721/95‑1 (MQ=255)
                                gAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:2266742/100‑1 (MQ=255)
                                gAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:5669340/100‑1 (MQ=255)
                                gAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:2630675/100‑1 (MQ=255)
                                gAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:2550640/100‑1 (MQ=255)
                                gAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:3466876/100‑1 (MQ=255)
                                gAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:3504546/100‑1 (MQ=255)
                                cttccgatcTGTGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:2437225/92‑1 (MQ=255)
                                cttccgatcTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:5664264/92‑1 (MQ=255)
                                cttccgatcTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:813002/92‑1 (MQ=255)
                                cgatctgccTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCa                                                              <  1:4153720/93‑1 (MQ=255)
                                 ttccgatcTGTGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCaa                                                             <  1:4407625/93‑1 (MQ=255)
                                 ttccgatcTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCaa                                                             <  1:443137/93‑1 (MQ=255)
                                 aaaaGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCCCCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCaa                                                             <  1:3078353/100‑1 (MQ=255)
                                 aaaaGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAAGTTCATTCCCATCTTTTCACCaa                                                             <  1:4801595/100‑1 (MQ=255)
                                  aaaGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCagc                                                            <  1:5036087/100‑3 (MQ=255)
                                  aaaGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCaat                                                            <  1:1017470/100‑2 (MQ=255)
                                   aaGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCaagt                                                           <  1:2749647/100‑3 (MQ=255)
                                       tcTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTg                                                       <  1:3943249/99‑1 (MQ=255)
                                       tcTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTg                                                       <  1:31823/99‑1 (MQ=255)
                                       tcTGCGTCTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAGATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTg                                                       <  1:4755436/99‑1 (MQ=255)
                                                   ttATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGt                                           <  1:4979413/100‑1 (MQ=255)
                                                   ttATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGt                                           <  1:2206407/100‑1 (MQ=255)
                                                   ttATCATATGTAATAACACCAATGTTAGTTGTCACGCCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGt                                           <  1:1173822/100‑1 (MQ=255)
                                                    tATCATATGTAACAACACCAATGTTAGTTGTCACGTCTCCCATGACATCTGTATGGAATTTCACTCCCATCTTTTCACCAAAATCTGAATTAAATAAGtt                                          <  1:1250262/100‑2 (MQ=255)
                                                      tCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGc                                        <  1:5370600/100‑1 (MQ=255)
                                                      tCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGc                                        <  1:2991837/100‑1 (MQ=255)
                                                       cATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCg                                       <  1:3048903/100‑1 (MQ=255)
                                                       cATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCg                                       <  1:2667354/100‑1 (MQ=255)
                                                       cATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCg                                       <  1:4083411/100‑1 (MQ=255)
                                                       cATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCg                                       <  1:5131977/100‑1 (MQ=255)
                                                        atatGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGa                                      <  1:3715112/100‑1 (MQ=255)
                                                        atatGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGa                                      <  1:3609743/100‑1 (MQ=255)
                                                         tatGTAATAACACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGAt                                     <  1:3039735/100‑1 (MQ=255)
                                                              aataaCACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTa                                <  1:2257896/100‑1 (MQ=255)
                                                               ataaCACCAATGTTAGTTGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTag                               <  1:391163/100‑2 (MQ=255)
                                                                               ttGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTtg               <  1:5671700/100‑1 (MQ=255)
                                                                               ttGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTtg               <  1:4869128/100‑1 (MQ=255)
                                                                                tGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTtgt              <  1:5618552/100‑1 (MQ=255)
                                                                                tGTCACGTCACCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTtgt              <  1:1271321/100‑1 (MQ=255)
                                                                                 gtcacgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTtgtg             <  1:295691/100‑1 (MQ=255)
                                                                                  tccgatctcccTTGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTtgtgt            <  1:1771472/92‑2 (MQ=255)
                                                                                  tcccgtctcCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:4584768/97‑1 (MQ=255)
                                                                                  tcacgtcgcCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:625032/100‑1 (MQ=255)
                                                                                  tcacgtcacCCATGAAATCTGTATGGAATTTTATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:3484601/100‑1 (MQ=255)
                                                                                  tcacgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:1727840/100‑1 (MQ=255)
                                                                                  tcacgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:4748662/100‑1 (MQ=255)
                                                                                  tcacgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:5435594/100‑1 (MQ=255)
                                                                                  tcacgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:5394015/100‑1 (MQ=255)
                                                                                  tcacgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:3437285/100‑1 (MQ=255)
                                                                                  tcacgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:4282087/100‑1 (MQ=255)
                                                                                  cgatctcacCCGTGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:51073/95‑1 (MQ=255)
                                                                                  ccgatctacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:4925208/93‑1 (MQ=255)
                                                                                  ccgatctacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGc            <  1:3772257/93‑1 (MQ=255)
                                                                                   cacgtcccCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCa           <  1:4707096/100‑1 (MQ=255)
                                                                                   cacgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCa           <  1:4944557/100‑1 (MQ=255)
                                                                                   cacgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCa           <  1:2249278/100‑1 (MQ=255)
                                                                                   atctccacTCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCa           <  1:27423/95‑1 (MQ=255)
                                                                                    acgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCag          <  1:3843385/100‑2 (MQ=255)
                                                                                    acgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCag          <  1:393896/100‑2 (MQ=255)
                                                                                    acgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCac          <  1:492468/100‑2 (MQ=255)
                                                                                     cgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCatt         <  1:3715544/100‑1 (MQ=255)
                                                                                     cgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCatt         <  1:4776906/100‑1 (MQ=255)
                                                                                     cgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCatt         <  1:2008382/100‑1 (MQ=255)
                                                                                     cgtcacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCatg         <  1:3240506/100‑2 (MQ=255)
                                                                                      tcttcccATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACACTTGTGCattt        <  1:3710803/96‑2 (MQ=255)
                                                                                        cacCCATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCattatt      >  1:4173038/1‑100 (MQ=255)
                                                                                            cATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTATATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCATTATTATCa  <  1:825713/100‑1 (MQ=255)
                                                                                            cATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCATtattataa  <  1:3465896/100‑3 (MQ=255)
                                                                                            cATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCATTATTATCa  <  1:692105/100‑1 (MQ=255)
                                                                                            cATGAAATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCATTATTATCa  <  1:5258337/100‑1 (MQ=255)
                                                                                                 |                                                                                              
TCAAATCCTTCTGGGTAGCGTAAGTTGATACCGAAAAGACCTGCGTTTTCATTATCATATGTAATAACACCAATGTTAGTTGTCACGTCACCCATGACATCTGTATGGAATTTCATTCCCATCTTTTCACCAAAATCTGAATTAAATAAGTAGCGATTACTAAATGCTACAAACGCTTGTGCATTATTATCA  >  NC_007793/1871508‑1871699

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑