Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,202,110 T→C 81.5% F611S (TTT→TCT)  yehI → DUF4132 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,202,1100TC81.5% 215.2 / 38.0 98F611S (TTT→TCT) yehIDUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (13/5);  new base C (50/28);  total (64/34)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01

Only 100 of 101 total aligned reads displayed.

ACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGGCGCAAC  >  NC_000913/2202014‑2202209
                                                                                                |                                                                                                   
aCCTGTATTTATTGGCACAGTTGCCGCGTGAGGCGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAACAGCATCTGTTCGCGGGTGAGGACTATtctc                                                                                                   >  2:254509/1‑99 (MQ=255)
 ccTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATtctctc                                                                                                 <  1:1372379/100‑1 (MQ=255)
   tGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAg                                                                                               >  1:1342471/1‑100 (MQ=255)
   tGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAg                                                                                               <  1:1670762/100‑1 (MQ=255)
   tGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCATATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAg                                                                                               >  2:1944170/1‑100 (MQ=255)
    gTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGCGGACTATTCTCTCAGt                                                                                              >  1:2129276/1‑100 (MQ=255)
    gTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGt                                                                                              >  2:1608089/1‑100 (MQ=255)
    gTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGt                                                                                              >  2:1593067/1‑100 (MQ=255)
     tatttattGGCACAGTTGCCGCGTGCGATGGCTTTGACCTTCTGGATGCGTATCAATGAGAACAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTa                                                                                             >  1:1010632/1‑100 (MQ=255)
     tatttattGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTa                                                                                             >  1:1261735/1‑100 (MQ=255)
      atttattGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTAt                                                                                            >  1:2863812/1‑100 (MQ=255)
      atttattGGCACAGTGGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTa                                                                                             <  1:2771617/99‑1 (MQ=255)
        ttattGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATcc                                                                                          >  2:76567/1‑100 (MQ=255)
          attGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTc                                                                                        >  2:406447/1‑100 (MQ=255)
             ggCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGa                                                                                     <  2:391817/100‑1 (MQ=255)
             ggCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGa                                                                                     <  1:379682/100‑1 (MQ=255)
               cacaGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGAtt                                                                                   >  1:798615/1‑100 (MQ=255)
               cacaGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTATCTCAGTATCCTCGGCtt                                                                                   >  2:1995534/1‑100 (MQ=255)
                 caGTTGCCGCGTGAGATGGCTTTGACATTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTgg                                                                                 >  1:2839425/1‑100 (MQ=255)
                   gTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGAt                                                                               <  1:1979880/100‑1 (MQ=255)
                     tGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATgc                                                                             <  1:1662480/100‑1 (MQ=255)
                     tGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATgc                                                                             >  1:474434/1‑100 (MQ=255)
                      gCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATgcg                                                                            >  1:438968/1‑100 (MQ=255)
                      gCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATgcg                                                                            >  2:1894302/1‑100 (MQ=255)
                      gCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATgcg                                                                            <  1:1108931/100‑1 (MQ=255)
                       ccGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATgcgc                                                                           <  2:1159577/100‑1 (MQ=255)
                         gcgTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTa                                                                         <  2:312992/100‑1 (MQ=255)
                              agaTGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAgg                                                                    >  1:35492/1‑100 (MQ=255)
                               gaTGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAgg                                                                    >  1:1806840/1‑99 (MQ=255)
                               gaTGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGt                                                                   >  2:1044914/1‑100 (MQ=255)
                               gaTGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGt                                                                   >  2:1004938/1‑100 (MQ=255)
                                 tGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTct                                                                 >  1:969910/1‑100 (MQ=255)
                                  ggCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTctg                                                                >  2:335881/1‑100 (MQ=255)
                                   gCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTctgc                                                               >  2:2615114/1‑100 (MQ=255)
                                    cTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGGTTGGATGCGCTACCAGGTctg                                                                >  2:2707284/1‑98 (MQ=255)
                                    cTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTctgct                                                              >  1:1711749/1‑100 (MQ=255)
                                       tGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGtt                                                           <  2:2732040/100‑1 (MQ=255)
                                        gACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTg                                                          <  2:401536/100‑1 (MQ=255)
                                        gACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATACTCGGCTTGGATGCGCTACCAGGTCTGCTGTTg                                                          >  1:547345/1‑100 (MQ=255)
                                         aCCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTgg                                                         <  2:502162/100‑1 (MQ=255)
                                          ccTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGc                                                        >  2:2637161/1‑100 (MQ=255)
                                           cTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCt                                                       >  1:555587/1‑100 (MQ=255)
                                             tatGGCTGCGTATCAATGTGTAAAATCATCTGTTAGCGGGTGCGTACTATTTTTTCAGTATCCTCGGATTGGTTGCGCTACCAGGTCTGCTGTTGGCttt                                                     >  1:1875132/3‑100 (MQ=255)
                                             tCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCttt                                                     <  2:235259/100‑1 (MQ=255)
                                             tCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCttt                                                     <  1:2136439/100‑1 (MQ=255)
                                               tGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCttttt                                                   >  2:17054/1‑100 (MQ=255)
                                               tGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCttttt                                                   >  1:80219/1‑100 (MQ=255)
                                                ggCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTc                                                  >  1:1863227/1‑100 (MQ=255)
                                                ggCTGCGTATCAATGAGAAAAAGCATCCGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTc                                                  >  1:2609104/1‑100 (MQ=255)
                                                 gCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTca                                                 <  1:2824643/100‑1 (MQ=255)
                                                  cTGCGTATCAATGTGAAAAAGCATCTGTTCGCGGGTGAGGACTAGTCTCTCAGTATGCTCGGATTGAATGCGCTAACAGGGCTGATATTGGCtttgtacc                                                >  1:2514417/1‑95 (MQ=255)
                                                  cTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTcac                                                <  2:798618/100‑1 (MQ=255)
                                                    gCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACAt                                              >  1:1167387/1‑100 (MQ=255)
                                                     cGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATc                                             >  2:197134/1‑100 (MQ=255)
                                                     cGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATc                                             <  1:1873196/100‑1 (MQ=255)
                                                     cGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATc                                             >  2:1688545/1‑100 (MQ=255)
                                                     cGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATc                                             >  1:1344251/1‑100 (MQ=255)
                                                       tATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGCCTGCTGTTGGCTTTTTCACATc                                             >  2:2804349/1‑98 (MQ=255)
                                                        aTCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTc                                          <  1:1044911/100‑1 (MQ=255)
                                                        aTCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTc                                          >  1:753138/1‑100 (MQ=255)
                                                         tCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTcc                                         >  1:1664069/1‑100 (MQ=255)
                                                         tCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATATTCTCAGTATCCTCGGATTGGACGCGCTACCAGGTCTACTGTTGGCATTTTCACATCGccc                                         >  1:2770497/1‑100 (MQ=255)
                                                          cAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCa                                        >  1:1972158/1‑100 (MQ=255)
                                                          cAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCa                                        >  2:2168165/1‑100 (MQ=255)
                                                          cAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCa                                        >  2:2087897/1‑100 (MQ=255)
                                                            aTGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCaaa                                      <  2:1471701/100‑1 (MQ=255)
                                                            aTGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCaaa                                      >  1:1657834/1‑100 (MQ=255)
                                                            aTGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCaa                                       >  1:646696/1‑99 (MQ=255)
                                                             tGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCaaaa                                     >  2:1238615/1‑100 (MQ=255)
                                                             tGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGAGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTaaaaaa                                     <  1:1498428/100‑1 (MQ=255)
                                                              ttccgatctcATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAg                                    <  1:2427273/91‑1 (MQ=255)
                                                              gagaAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCaaaa                                     <  2:2432157/99‑1 (MQ=255)
                                                               agaAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCCAAAGa                                   >  2:919146/1‑100 (MQ=255)
                                                               agaAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGa                                   >  1:935453/1‑100 (MQ=255)
                                                                gaAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGaa                                  >  1:2036467/1‑100 (MQ=255)
                                                                gaAAAAGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATACTCTGATTGGATGCGCTACCAGGTCTGCTCTTGGCTTTTTCACATCGTCCAAAAGaa                                  >  1:1293300/1‑100 (MQ=255)
                                                                  aaaaGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAAc                                <  1:1658313/100‑1 (MQ=255)
                                                                   aaaGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACa                               <  2:2696939/100‑1 (MQ=255)
                                                                     aGCATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACAtt                             <  1:1000448/100‑1 (MQ=255)
                                                                       cATCTGTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGagatcggaa                           >  2:2427288/1‑92 (MQ=255)
                                                                         tCTGTTCGCGGGTGAGGCCTATTATCTCAGTATCCTCGGATTGGATGCGCTACCCGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCg                         <  1:1995523/100‑1 (MQ=255)
                                                                            gTTCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTa                      >  2:2478806/1‑100 (MQ=255)
                                                                             ttCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTaa                     >  1:2069562/1‑100 (MQ=255)
                                                                              tCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTGACATCGTCCAAAAGAAACATTTCCGTTAAt                    >  1:116374/1‑100 (MQ=255)
                                                                              tCGCGGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAAt                    >  2:231755/1‑100 (MQ=255)
                                                                                gcgGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAAttt                  <  2:2200833/100‑1 (MQ=255)
                                                                                gcgGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAAttt                  <  2:2099026/100‑1 (MQ=255)
                                                                                 cgGGTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAAtttt                 >  1:841361/1‑100 (MQ=255)
                                                                                    gTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTaaa              >  2:557848/1‑100 (MQ=255)
                                                                                    gTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTaaa              <  2:1565132/100‑1 (MQ=255)
                                                                                    gTGAGGACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTaaa              >  2:2849174/1‑100 (MQ=255)
                                                                                       aGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAAttt           >  2:1416848/1‑100 (MQ=255)
                                                                                         gACTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTg         <  1:1396366/100‑1 (MQ=255)
                                                                                          aCTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTgg        <  2:2721023/100‑1 (MQ=255)
                                                                                          aCTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTgg        >  2:2236055/1‑100 (MQ=255)
                                                                                           cTATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGgc       <  1:1944161/100‑1 (MQ=255)
                                                                                            tATTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGgcg      <  2:1261740/100‑1 (MQ=255)
                                                                                             aTTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGgcg      >  1:529454/1‑99 (MQ=255)
                                                                                             aTTCTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGgcgc     >  1:1854097/1‑100 (MQ=255)
                                                                                                ttCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGGCGCAAc  <  1:1416843/100‑1 (MQ=255)
                                                                                                |                                                                                                   
ACCTGTATTTATTGGCACAGTTGCCGCGTGAGATGGCTTTGACCTTCTGGCTGCGTATCAATGAGAAAAAGCATCTGTTCGCGGGTGAGGACTATTTTCTCAGTATCCTCGGATTGGATGCGCTACCAGGTCTGCTGTTGGCTTTTTCACATCGTCCAAAAGAAACATTTCCGTTAATTTTAAATTTTGGCGCAAC  >  NC_000913/2202014‑2202209

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