Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,356 Δ1 bp coding (196/993 nt) rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,3550C.98.4% 295.1 / ‑1.1 62coding (197/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (31/30);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCAT  >  NC_000913/2867208‑2867499
                                                                                                                                                   |                                                                                                                                                
gCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCCGACGCAAGTTACTCTCGATCATCCGGCGGCGGGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGGGAATAAcaa                                                                                                                                                >  1:848571/1‑148 (MQ=255)
  cAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCCAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGGGAATAAcaat                                                                                                                                               >  1:884160/1‑146 (MQ=255)
    gACCACGATTGCCATAACGGCGGGCAATTTTTGCCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATACACTTCTTCTGCGGCCGTTAACAGTGGTGAATAAcaacctca                                                                                                                                           <  2:327848/151‑8 (MQ=255)
    gACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCa                                                                                                                                          <  2:193973/151‑1 (MQ=255)
     aCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACGATCTCAc                                                                                                                                         <  1:20010/151‑1 (MQ=255)
      ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCa                                                                                                                                          <  2:1089507/149‑1 (MQ=255)
      ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTAATCTCGATCATCCGGCGGCGAGAAGCGACATCTCCCCGCAGCGGGCGGCGCGCCAAATAAACTTCTTCTTCGGCCGTTAACAGGGGTGAATAACAAACTCAcc                                                                                                                                        >  1:217703/1‑151 (MQ=255)
           aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAgg                                                                                                                                    <  1:19867/150‑1 (MQ=255)
             tGCCATAACGGCGGGCAATTTTTACCACCCGACGCAAGTTAATCTCGATCATCCGGCGGCGGGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAAAAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAgg                                                                                                                                    >  2:1033042/1‑148 (MQ=255)
              gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTaaa                                                                                                                                <  2:395527/151‑1 (MQ=255)
                     cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGt                                                                                                                         <  1:408058/151‑1 (MQ=255)
                        cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCTGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTcgc                                                                                                                      >  2:880267/1‑151 (MQ=255)
                             aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCa                                                                                                                 <  2:121521/151‑1 (MQ=255)
                                      caccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                                                                                                        <  2:229330/151‑1 (MQ=255)
                                        ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgt                                                                                                      >  1:355405/1‑151 (MQ=255)
                                        ccaGACGCAAGTTACTCGCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgc                                                                                                      <  2:379362/151‑2 (MQ=255)
                                         caGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtg                                                                                                     >  1:1072495/1‑151 (MQ=255)
                                            aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTgg                                                                                                  <  2:632404/151‑1 (MQ=255)
                                               cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTc                                                                                               >  1:633536/1‑151 (MQ=255)
                                                 aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTccc                                                                                             <  2:185287/151‑1 (MQ=255)
                                                      ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCg                                                                                         >  2:579885/1‑150 (MQ=255)
                                                      ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCg                                                                                         <  2:175191/150‑1 (MQ=255)
                                                        ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                                        <  1:936838/149‑1 (MQ=255)
                                                           gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                                    >  1:404089/1‑150 (MQ=255)
                                                           gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                                    >  2:96696/1‑150 (MQ=255)
                                                            atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                                    >  1:822494/1‑149 (MQ=255)
                                                              catcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGtt                                                                                >  1:809496/1‑151 (MQ=255)
                                                                  cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCct                                                                             >  2:767772/1‑150 (MQ=255)
                                                                     cggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                                         <  1:342478/151‑1 (MQ=255)
                                                                     cggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                                         <  2:362513/151‑1 (MQ=255)
                                                                      ggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                                         >  1:60131/1‑151 (MQ=255)
                                                                       gcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCgg                                                                       >  1:253510/1‑151 (MQ=255)
                                                                         gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGcc                                                                     >  2:597482/1‑151 (MQ=255)
                                                                            agGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCaa                                                                   <  2:598309/150‑1 (MQ=255)
                                                                              gCGGCATCACCACGCAGTGCTCGACGCGCAAAATAAACTTCTTCTTCGAACGTAAACAGTGGTGAATAACAAACTCACCAATGTAAAGCTGGGGCGCGTCCAAAACGCGCTGTGTGGGTCCCTGAGATAACAGATCCTCTTc                                                                         <  2:735667/142‑1 (MQ=255)
                                                                              gCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATc                                                                <  1:248518/151‑1 (MQ=255)
                                                                                gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGt                                                              >  2:81201/1‑151 (MQ=255)
                                                                                   aTCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAt                                                           >  2:961139/1‑151 (MQ=255)
                                                                                       ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACt                                                       <  1:676384/151‑1 (MQ=255)
                                                                                        cACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTg                                                      <  2:813731/151‑1 (MQ=255)
                                                                                                gcgcGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTg                                              >  1:674368/1‑151 (MQ=255)
                                                                                                gcgcGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTg                                              <  1:650916/151‑1 (MQ=255)
                                                                                                    gACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTct                                          >  2:911595/1‑151 (MQ=255)
                                                                                                      cgcgCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTcta                                         >  2:994380/1‑150 (MQ=255)
                                                                                                       gcgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctact                                       >  2:34815/1‑151 (MQ=255)
                                                                                                        cgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctacta                                      <  1:913782/151‑1 (MQ=255)
                                                                                                        cgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctact                                       >  1:160876/1‑150 (MQ=255)
                                                                                                         gcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctacta                                      >  2:786741/1‑150 (MQ=255)
                                                                                                          cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctacta                                      >  1:1039071/1‑149 (MQ=255)
                                                                                                           aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAgg                                   <  2:344289/151‑1 (MQ=255)
                                                                                                           aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAgg                                   >  2:293345/1‑151 (MQ=255)
                                                                                                            aaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGc                                  >  1:804077/1‑151 (MQ=255)
                                                                                                             aaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGcc                                 >  2:758504/1‑151 (MQ=255)
                                                                                                                aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCttt                              >  1:697981/1‑151 (MQ=255)
                                                                                                                aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCtt                               <  1:1066151/150‑1 (MQ=255)
                                                                                                                 aaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCtttt                             >  2:338931/1‑151 (MQ=255)
                                                                                                                          ttcGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaaa                     >  2:847258/1‑150 (MQ=255)
                                                                                                                           tcGGCCGTTAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaaa                     <  2:467785/149‑1 (MQ=255)
                                                                                                                                   tAACAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt            <  2:871830/150‑1 (MQ=255)
                                                                                                                                    aaCAGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc          <  1:187586/151‑1 (MQ=255)
                                                                                                                                       aGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt       <  1:579892/151‑1 (MQ=255)
                                                                                                                                       aGTGGTGAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt       <  1:16832/151‑1 (MQ=255)
                                                                                                                                          ggtgAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCTAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtctc    <  1:6307/151‑1 (MQ=255)
                                                                                                                                           gtgAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtctc    >  2:31604/1‑150 (MQ=255)
                                                                                                                                            tgAATAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcat  >  2:582089/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCAT  >  NC_000913/2867208‑2867499

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