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)
    gACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCa                                                                                                                                          <  2:193973/151‑1 (MQ=255)
     aCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CGATCTCAc                                                                                                                                         <  1:20010/151‑1 (MQ=255)
      ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCa                                                                                                                                          <  2:1089507/149‑1 (MQ=255)
      ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTAATCTCGATCATCCGGCGGCGAGAAGCGACATCTCCCCGCAGCGGGCGGCGCGCCAAATAAACTTCTTCTTCGGCCGTTAACAGGGGTGAATAA‑CAAACTCAcc                                                                                                                                        >  1:217703/1‑151 (MQ=255)
           aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAgg                                                                                                                                    <  1:19867/150‑1 (MQ=255)
             tGCCATAACGGCGGGCAATTTTTACCACCCGACGCAAGTTAATCTCGATCATCCGGCGGCGGGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAAAAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAgg                                                                                                                                    >  2:1033042/1‑148 (MQ=255)
              gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTaaa                                                                                                                                <  2:395527/151‑1 (MQ=255)
                     cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGt                                                                                                                         <  1:408058/151‑1 (MQ=255)
                        cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCTGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTcgc                                                                                                                      >  2:880267/1‑151 (MQ=255)
                             aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCa                                                                                                                 <  2:121521/151‑1 (MQ=255)
                                      caccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                                                                                                        <  2:229330/151‑1 (MQ=255)
                                        ccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgt                                                                                                      >  1:355405/1‑151 (MQ=255)
                                        ccaGACGCAAGTTACTCGCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgc                                                                                                      <  2:379362/151‑2 (MQ=255)
                                         caGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtg                                                                                                     >  1:1072495/1‑151 (MQ=255)
                                            aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTgg                                                                                                  <  2:632404/151‑1 (MQ=255)
                                               cAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTc                                                                                               >  1:633536/1‑151 (MQ=255)
                                                 aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTccc                                                                                             <  2:185287/151‑1 (MQ=255)
                                                      ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCg                                                                                         >  2:579885/1‑150 (MQ=255)
                                                      ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCg                                                                                         <  2:175191/150‑1 (MQ=255)
                                                        ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa                                                                                        <  1:936838/149‑1 (MQ=255)
                                                           gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                                    >  1:404089/1‑150 (MQ=255)
                                                           gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                                    >  2:96696/1‑150 (MQ=255)
                                                            atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc                                                                                    >  1:822494/1‑149 (MQ=255)
                                                              catcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGtt                                                                                >  1:809496/1‑151 (MQ=255)
                                                                  cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCct                                                                             >  2:767772/1‑150 (MQ=255)
                                                                     cggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                                         <  1:342478/151‑1 (MQ=255)
                                                                     cggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                                         <  2:362513/151‑1 (MQ=255)
                                                                      ggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTc                                                                         >  1:60131/1‑151 (MQ=255)
                                                                       gcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCgg                                                                       >  1:253510/1‑151 (MQ=255)
                                                                         gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGcc                                                                     >  2:597482/1‑151 (MQ=255)
                                                                            agGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCaa                                                                   <  2:598309/150‑1 (MQ=255)
                                                                              gCGGCATCACCACGCAGTGCTCGACGCGCAAAATAAACTTCTTCTTCGAACGTAAACAGTGGTGAATAA‑CAAACTCACCAATGTAAAGCTGGGGCGCGTCCAAAACGCGCTGTGTGGGTCCCTGAGATAACAGATCCTCTTc                                                                         <  2:735667/142‑1 (MQ=255)
                                                                              gCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATc                                                                <  1:248518/151‑1 (MQ=255)
                                                                                gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGt                                                              >  2:81201/1‑151 (MQ=255)
                                                                                   aTCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAt                                                           >  2:961139/1‑151 (MQ=255)
                                                                                       ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACt                                                       <  1:676384/151‑1 (MQ=255)
                                                                                        cACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTg                                                      <  2:813731/151‑1 (MQ=255)
                                                                                                gcgcGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTg                                              >  1:674368/1‑151 (MQ=255)
                                                                                                gcgcGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTg                                              <  1:650916/151‑1 (MQ=255)
                                                                                                    gACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTct                                          >  2:911595/1‑151 (MQ=255)
                                                                                                      cgcgCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTcta                                         >  2:994380/1‑150 (MQ=255)
                                                                                                       gcgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctact                                       >  2:34815/1‑151 (MQ=255)
                                                                                                        cgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctacta                                      <  1:913782/151‑1 (MQ=255)
                                                                                                        cgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctact                                       >  1:160876/1‑150 (MQ=255)
                                                                                                         gcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctacta                                      >  2:786741/1‑150 (MQ=255)
                                                                                                          cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTctacta                                      >  1:1039071/1‑149 (MQ=255)
                                                                                                           aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAgg                                   <  2:344289/151‑1 (MQ=255)
                                                                                                           aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAgg                                   >  2:293345/1‑151 (MQ=255)
                                                                                                            aaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGc                                  >  1:804077/1‑151 (MQ=255)
                                                                                                             aaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGcc                                 >  2:758504/1‑151 (MQ=255)
                                                                                                                aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCttt                              >  1:697981/1‑151 (MQ=255)
                                                                                                                aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCtt                               <  1:1066151/150‑1 (MQ=255)
                                                                                                                 aaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCtttt                             >  2:338931/1‑151 (MQ=255)
                                                                                                                          ttcGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaaa                     >  2:847258/1‑150 (MQ=255)
                                                                                                                           tcGGCCGTTAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCaaaa                     <  2:467785/149‑1 (MQ=255)
                                                                                                                                   tAACAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt            <  2:871830/150‑1 (MQ=255)
                                                                                                                                    aaCAGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc          <  1:187586/151‑1 (MQ=255)
                                                                                                                                       aGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt       <  1:579892/151‑1 (MQ=255)
                                                                                                                                       aGTGGTGAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt       <  1:16832/151‑1 (MQ=255)
                                                                                                                                          ggtgAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCTAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtctc    <  1:6307/151‑1 (MQ=255)
                                                                                                                                           gtgAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtctc    >  2:31604/1‑150 (MQ=255)
                                                                                                                                            tgAATAA‑CAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcat  >  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: ‑