Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 217,044 A→G 15.1% intergenic (‑41/+13) yaeF ← / ← proS putative lipoprotein/prolyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913217,0440AG15.1% 348.8 / 44.7 144intergenic (‑41/+13)yaeF/proSputative lipoprotein/prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (61/61);  new base G (12/10);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

GCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATT  >  NC_000913/216945‑217142
                                                                                                   |                                                                                                  
gCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                                    >  2:2008392/1‑100 (MQ=255)
gCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                                    >  2:101869/1‑100 (MQ=255)
 ccGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGAc                                                                                                   <  1:2684337/100‑1 (MQ=255)
     aaatAAAAGAAATCTAGGAAAAAAAATGCCACATTACACCTGTGGTTTATCTATTGTATAAAAAAAAAAAAAAAAAAGGGCCCTGAAATTCAGGAccc                                                                                                 <  2:1333835/95‑1 (MQ=255)
      aGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCttt                                                                                              <  1:1159156/100‑1 (MQ=39)
       gTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTc                                                                                             >  2:3011595/1‑100 (MQ=255)
        tAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                            >  2:1096505/1‑100 (MQ=255)
         aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttt                                                                                           <  1:162417/100‑4 (MQ=38)
         aaaaGAAAGATTGGGAGAAACAAGCGGCAGTACGCCTATGGTTTATCCATTTTATACAATCTATGTAAAAAAATGACCCTGAAATTCAGGACCCTATc                                                                                             <  2:2905715/98‑1 (MQ=255)
          aaagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCttttttt                                                                                          <  1:1064607/100‑5 (MQ=38)
            agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttac                                                                                        <  1:2589262/100‑7 (MQ=37)
            agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCa                                                                                        <  1:3119514/100‑1 (MQ=255)
            agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTGTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttac                                                                                        <  2:1889462/100‑7 (MQ=37)
             gaaagCTGGTGAGAAACAGGCGGCAGTAAGCCTTTGGTTTATCCATTTTGTACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCttttttt                                                                                          <  2:3448666/97‑5 (MQ=37)
             gaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCAt                                                                                       >  2:2543152/1‑100 (MQ=255)
              aaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttacat                                                                                      <  1:2581185/100‑9 (MQ=37)
               aagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttacatg                                                                                     <  1:1141942/100‑10 (MQ=25)
               aagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCa                                                                                     >  2:3680984/1‑100 (MQ=255)
                agcatggcagaAACAAGCGACAGTACGCGTTTGGTTTATCCATTTAAAAAAAAAAAAAAAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCa                                                                                     <  2:302164/91‑1 (MQ=255)
                 gCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGc                                                                                   >  2:1062715/1‑100 (MQ=255)
                 gCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGc                                                                                   >  1:3488022/1‑100 (MQ=255)
                  cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                                  <  2:3405877/100‑1 (MQ=255)
                  cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                                  <  1:3342531/100‑1 (MQ=255)
                   tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                                  <  1:1348002/99‑1 (MQ=255)
                   tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCt                                                                                 >  2:1882834/1‑100 (MQ=255)
                       gagaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTaa                                                                             >  2:3446813/1‑100 (MQ=255)
                       gagaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTaa                                                                             >  2:601241/1‑100 (MQ=255)
                         gaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATc                                                                           >  2:1440589/1‑100 (MQ=255)
                         gaAACAGGCGGCAGTACGCCTTTGGTTTATCCAGTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATc                                                                           >  2:2601219/1‑100 (MQ=255)
                          aaaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCt                                                                          >  2:1530548/1‑100 (MQ=255)
                           aaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                         <  1:3202589/100‑1 (MQ=255)
                           aaCAGGCGGAAGTACGCCTTTGGTTTATCCATTTTATAAAAAAAAAAAAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                         <  2:3393078/100‑1 (MQ=255)
                             cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                                       <  1:1729697/100‑1 (MQ=255)
                               ggcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTc                                                                     >  2:1350872/1‑100 (MQ=255)
                                 cggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                                   >  2:1682217/1‑100 (MQ=255)
                                   gcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCa                                                                 <  1:1936995/100‑1 (MQ=255)
                                         cGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATAtt                                                           >  2:1820536/1‑100 (MQ=255)
                                           ccTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcg                                                         >  1:3501692/1‑100 (MQ=255)
                                           ccTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcg                                                         <  1:2427727/100‑1 (MQ=255)
                                           ccTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcg                                                         <  2:484947/100‑1 (MQ=255)
                                           ccTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcg                                                         >  1:523711/1‑100 (MQ=255)
                                            cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                        >  1:3017507/1‑100 (MQ=255)
                                            cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                        >  1:3332328/1‑100 (MQ=255)
                                            cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                        >  1:3318140/1‑100 (MQ=255)
                                            cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                        >  2:2864616/1‑100 (MQ=255)
                                              ttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacg                                                      <  2:691418/100‑1 (MQ=255)
                                                 gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATg                                                   <  2:147168/100‑1 (MQ=255)
                                                    tATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCa                                                <  2:2995188/100‑1 (MQ=255)
                                                    tATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCa                                                >  2:884305/1‑100 (MQ=255)
                                                     aTCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACTCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAc                                               >  2:3115078/1‑100 (MQ=255)
                                                     aTCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAc                                               <  1:37153/100‑1 (MQ=255)
                                                      tCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAcc                                              <  1:3479622/100‑1 (MQ=255)
                                                           tttATACAATCCATGTAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtt                                        <  2:1663014/100‑1 (MQ=255)
                                                           tttATACAATCCATGGAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCt                                         <  1:3881694/100‑1 (MQ=255)
                                                           tgcataCAATCTATCTAAAAAAAGGGCTCTGAAATTCAGGACCCTTTCTAGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCATTCt                                         <  2:3619383/97‑1 (MQ=255)
                                                            ttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtt                                        <  2:39747/100‑1 (MQ=255)
                                                             tataCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtta                                       <  2:1923850/100‑1 (MQ=255)
                                                             tataCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtta                                       <  1:3608781/100‑1 (MQ=255)
                                                              ntaCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATGTGTTTCACCAGATATTCGACGATGTCACCAGTCttaa                                      >  2:2665687/2‑100 (MQ=255)
                                                              ataCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaa                                      >  1:3304465/1‑100 (MQ=255)
                                                               taCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaat                                     <  2:3594473/100‑1 (MQ=255)
                                                               taCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaat                                     >  2:731777/1‑100 (MQ=255)
                                                                  aaTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaatt                                    <  1:1755298/98‑1 (MQ=255)
                                                                   aTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAAc                                 >  2:2580033/1‑100 (MQ=255)
                                                                   aTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAAc                                 <  1:1820545/100‑1 (MQ=255)
                                                                    tCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACt                                <  1:2782641/100‑1 (MQ=255)
                                                                     cgatctAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                               <  1:3007538/95‑1 (MQ=255)
                                                                     ccATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                               <  1:3572403/100‑1 (MQ=255)
                                                                      cATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGt                              >  2:2997883/1‑100 (MQ=255)
                                                                       aTGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGtt                             <  2:2787469/100‑1 (MQ=255)
                                                                        tGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGttt                            >  2:1798248/1‑100 (MQ=255)
                                                                         gTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGGCGATGTCACCAGTCTTAATTAACTGttt                            >  2:2130740/1‑99 (MQ=255)
                                                                         gTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtc                           >  2:1086374/1‑100 (MQ=255)
                                                                         gTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtc                           >  1:1547983/1‑100 (MQ=255)
                                                                          tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct                          >  1:2306049/1‑100 (MQ=255)
                                                                          gccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct                          >  1:748626/6‑100 (MQ=37)
                                                                           ccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                         >  1:647315/5‑100 (MQ=38)
                                                                           ccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                         >  1:3189315/5‑100 (MQ=38)
                                                                           ccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                         >  2:204365/5‑100 (MQ=38)
                                                                           aaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                         >  2:1912804/1‑100 (MQ=255)
                                                                            cagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCg                        >  2:1126883/4‑100 (MQ=38)
                                                                            aaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCg                        <  1:1375335/100‑1 (MQ=255)
                                                                             aaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGc                       <  2:2001884/100‑1 (MQ=255)
                                                                               aaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGc                       >  2:3750679/1‑98 (MQ=255)
                                                                               aaaGGGCCTTGAAATTCAGGGCCCTTTCTGGCAGCAGCCTTTAACCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGcc                      <  1:908985/99‑1 (MQ=255)
                                                                                  gggTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                  >  2:983621/1‑100 (MQ=255)
                                                                                  gggTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                  >  1:701882/1‑100 (MQ=255)
                                                                                  gggTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                  >  2:1598451/1‑100 (MQ=255)
                                                                                  gggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                  >  1:34651/1‑100 (MQ=255)
                                                                                  gggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                  >  1:2317061/1‑100 (MQ=255)
                                                                                   ggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc                 <  1:1220510/100‑1 (MQ=255)
                                                                                     cccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGTCGTTGcga               >  2:605401/1‑100 (MQ=255)
                                                                                       cTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgacg             <  1:62924/100‑1 (MQ=255)
                                                                                       cTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgacg             >  2:3163693/1‑100 (MQ=255)
                                                                                        tGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgacga            >  1:76749/1‑100 (MQ=255)
                                                                                           aaTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGatat         <  1:101869/100‑1 (MQ=255)
                                                                                              tCAGGACTCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTta      <  1:3115098/100‑1 (MQ=255)
                                                                                                aGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtata    <  1:2159195/100‑1 (MQ=255)
                                                                                                 ggACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTatat   >  2:3356121/1‑100 (MQ=255)
                                                                                                  gACCCTTTCTGGCATCAGCCTTTAATCTGTCTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATAtt  <  1:1902017/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATT  >  NC_000913/216945‑217142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: