Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 416,819 A→G V45A (GTT→GCT)  sbcD ← ATP dependent, structure specific DNA nuclease ‑ SbcD subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913416,8190AG99.6% 1062.1 / ‑5.3 245V45A (GTT→GCT) sbcDATP dependent, structure specific DNA nuclease ‑ SbcD subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (119/125);  total (120/125)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

CAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTG  >  NC_000913/416722‑416916
                                                                                                 |                                                                                                 
caGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGcagcaa                                                                                                >  1:617287/1‑101 (MQ=255)
 aGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCaat                                                                                               <  1:617142/101‑1 (MQ=255)
  gATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCaata                                                                                              <  1:616862/101‑1 (MQ=255)
  gATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCaata                                                                                              <  1:617050/101‑1 (MQ=255)
   aTGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACTGGCAGCaataa                                                                                             <  1:616677/101‑1 (MQ=255)
   aTGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCaataa                                                                                             >  1:617291/1‑101 (MQ=255)
      aCAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCg                                                                                          >  1:617295/1‑101 (MQ=255)
      aCAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCg                                                                                          >  2:617296/1‑101 (MQ=255)
         gCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCAt                                                                                       >  2:617300/1‑101 (MQ=255)
          ccAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATc                                                                                      <  1:616828/101‑1 (MQ=255)
            aGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATcc                                                                                     <  1:616972/100‑1 (MQ=255)
             gTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCAc                                                                                   >  1:617304/1‑101 (MQ=255)
              tttGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCAcc                                                                                  >  2:617307/1‑101 (MQ=255)
               ttGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCt                                                                                 <  1:616934/101‑1 (MQ=255)
               ttGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCt                                                                                 >  1:617309/1‑101 (MQ=255)
                 gCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGa                                                                               >  2:617311/1‑101 (MQ=255)
                    gTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATgg                                                                            >  2:617315/1‑101 (MQ=255)
                     tAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATggg                                                                           >  2:617317/1‑101 (MQ=255)
                      aaaTTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGt                                                                          >  1:617318/1‑101 (MQ=255)
                       aaTTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGtt                                                                         >  2:617323/1‑101 (MQ=255)
                        aTTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGttt                                                                        <  1:616943/101‑1 (MQ=255)
                        aTTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGttt                                                                        >  2:617325/1‑101 (MQ=255)
                          tGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTtgt                                                                      <  1:616857/101‑1 (MQ=255)
                             caacaaAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCt                                                                   >  1:617329/1‑101 (MQ=255)
                              aacaaAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTtgt                                                                      <  1:616804/97‑1 (MQ=255)
                              aacaaAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTg                                                                  >  1:617333/1‑101 (MQ=255)
                              aacaaAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTg                                                                  >  2:617332/1‑101 (MQ=255)
                                 aaaaCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGtct                                                               >  1:617338/1‑101 (MQ=255)
                                  aaaCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCt                                                                   >  2:617341/1‑96 (MQ=255)
                                   aaCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTcc                                                             <  1:616904/101‑1 (MQ=255)
                                   aaCGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTcc                                                             >  1:617342/1‑101 (MQ=255)
                                      ggTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCcagc                                                          >  2:617347/1‑101 (MQ=255)
                                      ggTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCcagc                                                          <  1:616921/101‑1 (MQ=255)
                                      ggTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCcagc                                                          >  1:617348/1‑101 (MQ=255)
                                       gTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCcagca                                                         <  1:617005/101‑1 (MQ=255)
                                           tataACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCa                                                     >  2:617351/1‑101 (MQ=255)
                                            ataACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAg                                                    >  1:617352/1‑101 (MQ=255)
                                                 gTGCGGGCGTAACTGGGGGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCa                                                                                    <  1:617210/64‑1 (MQ=255)
                                                   gcgcgcGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATc                                                                                           <  1:616743/55‑1 (MQ=255)
                                                   gCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGaa                                             >  1:617358/1‑101 (MQ=255)
                                                    cGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGaaa                                            >  2:617359/1‑101 (MQ=255)
                                                     gggCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGaaaa                                           >  2:617361/1‑101 (MQ=255)
                                                         gTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGgcagca                                                                                                 <  2:616775/43‑1 (MQ=255)
                                                          tAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTg                                      >  2:617364/1‑101 (MQ=255)
                                                          tAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTg                                      >  2:617366/1‑101 (MQ=255)
                                                            aCTGGGGGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATc                                                                                           <  2:616928/46‑1 (MQ=38)
                                                            aCTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCAt                                                                                       <  2:616895/50‑1 (MQ=255)
                                                            aCTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGt                                                                 <  2:616911/72‑1 (MQ=255)
                                                             cggggCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGtctc                                                              <  1:616951/72‑1 (MQ=255)
                                                               ggggggggAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTg                                      <  2:617107/96‑1 (MQ=255)
                                                               ggggggCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGtt                                 <  2:616878/101‑1 (MQ=255)
                                                               gggCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGtt                                 >  2:617370/1‑101 (MQ=255)
                                                               gggCGGCGAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCAAGCAGCCAg                                                    <  2:616935/82‑1 (MQ=255)
                                                                 gcggcgAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCa                               >  2:617372/1‑101 (MQ=255)
                                                                   ggggAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTtgtg                                                                     <  2:617122/58‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCa                                                                                    <  2:617022/43‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATg                                                                             <  2:616997/50‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGt                                                   <  2:616836/76‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTc                                                  <  2:617047/77‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCa                                                 <  2:616950/78‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCa                                                 <  1:616974/78‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGc                                         <  2:617097/86‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGa                                     <  2:616979/90‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGtt                                 <  2:616812/94‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCa                               <  1:616996/96‑1 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTc                          >  2:617376/1‑101 (MQ=255)
                                                                      gAGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTc                          >  1:617374/1‑101 (MQ=255)
                                                                       aGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGtt                                 <  2:617124/93‑1 (MQ=255)
                                                                       aGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTcg                         <  2:617189/101‑1 (MQ=255)
                                                                       aGCCGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTcg                         >  2:617379/1‑101 (MQ=255)
                                                                         ccGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTcgcg                       <  2:617070/101‑1 (MQ=255)
                                                                         ccGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTcgcg                       >  1:617383/1‑101 (MQ=255)
                                                                          cGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCgg                      >  1:617385/1‑101 (MQ=255)
                                                                          cGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCgg                      >  2:617388/1‑101 (MQ=255)
                                                                          cGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCgg                      <  1:617000/101‑1 (MQ=255)
                                                                          cGGTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCgg                      <  1:617010/101‑1 (MQ=255)
                                                                           ggTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGc                     >  1:617389/1‑101 (MQ=255)
                                                                           ggTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGc                     <  2:617013/101‑1 (MQ=255)
                                                                            gTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCt                    <  1:617173/101‑1 (MQ=255)
                                                                            gTATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCt                    <  1:617119/101‑1 (MQ=255)
                                                                             tATCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCtt                   >  2:617391/1‑101 (MQ=255)
                                                                              aTCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGt                                  >  2:617397/1‑85 (MQ=255)
                                                                              aTCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCttt                  >  2:617395/1‑101 (MQ=255)
                                                                               tCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCtttt                 >  1:617398/1‑101 (MQ=255)
                                                                               tCGAAAACATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCtttt                 <  1:616910/101‑1 (MQ=255)
                                                                                  aaaaCATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACt              <  1:616925/101‑1 (MQ=255)
                                                                                  aaaaCATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACt              >  2:617403/1‑101 (MQ=255)
                                                                                  aaaaCATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACt              >  2:617401/1‑101 (MQ=255)
                                                                                    aaCATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGt            <  1:616955/101‑1 (MQ=255)
                                                                                    aaCATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGt            <  2:617196/101‑1 (MQ=255)
                                                                                      cATCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGc                     >  1:617407/1‑90 (MQ=255)
                                                                                       aTCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGa         >  1:617409/1‑101 (MQ=255)
                                                                                        tCACCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCtttt                 >  2:617411/1‑92 (MQ=255)
                                                                                          aCCGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGt      >  2:617413/1‑101 (MQ=255)
                                                                                           ccGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGtt     <  2:617033/101‑1 (MQ=255)
                                                                                           ccGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGtt     <  1:617058/101‑1 (MQ=255)
                                                                                            cGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAg       <  2:616986/98‑1 (MQ=255)
                                                                                            cGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTc    >  1:617417/1‑101 (MQ=255)
                                                                                            cGGCAGCAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTc    <  2:617027/101‑1 (MQ=255)
                                                                                              gcagcaATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTg  <  2:617067/101‑1 (MQ=255)
                                                                                                 |                                                                                                 
CAGATGACAGCCAGTTTGCTGTAAATTGACAACAAAACGGTTGTATAACGTGCGGGCGTAACTGGGCGGCGAGCCGGTATCGAAAACATCACCGGCAACAATAATCGCATCCACCTGATGGGTTTGTGCTGTCTCCAGCAGCCAGTCAAGAAAAGCCTGATGTTCAGCTTCGCGGCTTTTACTGTAGAAGTTCTG  >  NC_000913/416722‑416916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: