Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5690CG28.2% 160.3 / 30.4 85intergenic (+35/‑161)hrpB/mrcBputative ATP‑dependent RNA helicase HrpB/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB
Reads supporting (aligned to +/- strand):  ref base C (26/35);  new base G (0/24);  total (26/59)
Fisher's exact test for biased strand distribution p-value = 3.86e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAT  >  NC_000913/164420‑164715
                                                                                                                                                       |                                                                                                                                                    
gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGCGAGGGcc                                                                                                                                                      >  1:1627762/1‑150 (MQ=255)
gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                      <  2:2456614/150‑1 (MQ=255)
gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCACTGTGGGAGAGGGcc                                                                                                                                                      >  1:788287/1‑150 (MQ=255)
 aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                                     <  1:785176/150‑1 (MQ=255)
   aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                   <  1:336055/150‑1 (MQ=255)
   aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                   <  1:678767/150‑1 (MQ=255)
   aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                   <  2:1179291/150‑1 (MQ=255)
   aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                   <  2:1680823/150‑1 (MQ=255)
   aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                   <  2:1680915/150‑1 (MQ=255)
       gCGTACCGTGCGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAGGAATTCGTAAGGTTGGTGTTCTCCCTCTCCCTCTGTGAGAGGGCCGGGGTGa                                                                                                                                              >  2:826330/1‑150 (MQ=255)
         ttACCGTGAGGTGAAAAAAGAGATTAAAGGGCTTTATCCCAAACTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                            <  1:1934027/149‑1 (MQ=39)
         gTACCGTTAGGTGCAAATAGTGATGAATGGGGGTTATCCCAAACTTGTCTGGCCGTACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                            <  2:2841310/150‑1 (MQ=38)
         gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                            <  2:44915/150‑1 (MQ=255)
          tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                           >  2:2720647/1‑150 (MQ=255)
           aCCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                                          >  1:592727/1‑150 (MQ=255)
             cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                        <  2:2828673/150‑1 (MQ=255)
             cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                        <  2:1126266/150‑1 (MQ=255)
              gTGAGGTGCAAAAAGAGATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                       <  2:1487042/150‑1 (MQ=255)
               tGAGGTGCAAACAGATATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                      <  2:297531/150‑1 (MQ=255)
               tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                      >  2:2268990/1‑150 (MQ=255)
                gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                     >  1:2457050/1‑150 (MQ=255)
                   gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                  >  1:681707/1‑150 (MQ=255)
                    tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCGCAgcgcg                                                                                                                                 >  1:530307/1‑150 (MQ=255)
                    tGCAAAAAGAGATGAAAGGGCGTTAACCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                 >  1:1308420/1‑150 (MQ=255)
                       aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                              >  2:1261090/1‑150 (MQ=255)
                       aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                              >  1:502308/1‑150 (MQ=255)
                        aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                             >  2:1465690/1‑150 (MQ=255)
                         aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                            >  2:1711096/1‑150 (MQ=255)
                              tttGAATGGGCGTTATCCCAAACAGGTCTGGCGGGACGACCGGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCTTACGGTTGGTTTTCTCCCTCTCCCTGTGGGATCGGGCGGGGTTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                       <  2:2969609/148‑1 (MQ=255)
                               aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACggggggcc                                                                                                                      >  1:1086738/1‑143 (MQ=255)
                                tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCACGTggacccc                                                                                                                     >  2:1892816/1‑149 (MQ=255)
                                  aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCAGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccccc                                                                                                                   >  1:2365960/1‑147 (MQ=255)
                                    aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccgcacc                                                                                                                 >  2:1408401/1‑150 (MQ=255)
                                            aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCGCCCCaaacac                                                                                                         >  1:1495721/1‑146 (MQ=38)
                                                aaaccaggTCTGGCCGGACGACCGGGAAAATCCTGCACCGACGCTACGGCGGAAAAAGTTTTCTAAAGGTTGGTTTGCTCCCTCTCCCGGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                    <  2:2576309/143‑1 (MQ=255)
                                                     aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCACACCCACTCCCACaa                                                                                               >  2:2743000/1‑150 (MQ=37)
                                                      ggTCGGGCCGGCCGCCCGGGCAAATCCTCCCCCGCCGCGAGGCACTAAAATTTATCCGTATGGTTGTTTTCCGCCCTCTCCGTGTGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                              <  2:1493594/149‑1 (MQ=255)
                                                       gacgggCCGGACGCCCCGGAAAATACTGCACCTACGCGACGGACTAAAAAGTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                             <  2:2748185/146‑1 (MQ=37)
                                                          gggCCGCACGACCGGTCAAATAGTGCCCCGCCGCCAGGGACTAAAATGTTTTCGAATGGTTGGTTTCCGCCCGCCCCCTGTGGGTGCGGGCCGGGGTGTGGCCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                          <  2:1151501/149‑1 (MQ=255)
                                                          ggcccGGACCACCGGGAAATACCTCCCCCGCCGCGCGGGCCTAAAAATTTTTCGTTGGGTGGTTTTACCCCCCCCCCCGGGGGGAGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                          <  1:2134228/147‑1 (MQ=255)
                                                           ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGTTTGGTTTTCTCCCTCCCCCTGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                         <  2:2372280/150‑1 (MQ=255)
                                                            gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGATGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                        <  1:1029469/150‑1 (MQ=255)
                                                             ccccccgacCCGCAAAATACTGCACCGACGCCACGGACGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                       <  2:2136420/145‑1 (MQ=38)
                                                             ccGGACGCCCCGGAAAATACGGCCCCGACGCGACGGACTAAAATTTTTTCTTAGGTTTGTTTTACTCCCTCCCCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                       <  1:2901503/150‑1 (MQ=255)
                                                             ccGCACGCCCCCGCAAATCCGCCCCCGACGGGCGGCACGAAAAAGTATTCGAACGGTTGGTTTTCTCCCTCACCCTGGGGGAGAGGCCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                       <  1:2901611/150‑1 (MQ=255)
                                                              cGGACGCCCCGGCAATTACGGCCCCGCCGCGACGGACGAAAAACTTTTCGTAAGGTTGGTTTTCCCCCACACCCGGGGGGAGTGGCCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                      <  2:2555070/150‑1 (MQ=255)
                                                              cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAACCCAAACCCACCACGGAGAggg                                                                                      >  2:1877227/1‑150 (MQ=37)
                                                              cGGACGACCCGCCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCACCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                      <  2:2955863/150‑1 (MQ=255)
                                                               gtaccccccGGCAAATACTGCACCGACGCCCCGCACTAAATTTTTTTCTAATGTTTGTTTTTCACCACCGCCCTGGGGTTGTGGGCGGGGTTTGGCGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggg                                                                                     <  2:681707/144‑1 (MQ=255)
                                                                  ccgaccGGCAAAATACGCCCCCGACGCGACGGACGAAAAAGTATTCGAACGGTTGTTTTACCCCCACCCCCGGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                   <  2:739798/146‑1 (MQ=255)
                                                                      ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTGGCCCCCACCCCAACCCACAACCCCAACGGGGAGGGGGCCgcac                                                                              >  1:2821670/1‑147 (MQ=25)
                                                                         gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCCAACCCGCGCCCCCCAGGGGGGTGGGGCCGCTCgcg                                                                           >  2:2882355/1‑148 (MQ=25)
                                                                          cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCCCCCAAAACCCCACCCGCCAGGGGGAGGGGGCCGATCGAGc                                                                          >  2:2717309/1‑150 (MQ=25)
                                                                            aaaTACTGCACCGACGCGCCGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                         <  2:141207/150‑1 (MQ=255)
                                                                                     cACCGAGGCGACGTACTAAAATTTTTTGTTATGTTTGGTTTTCGCCCGCCCCGGGGGGGAGTGGGCGGGGGTTCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                 <  1:2514112/150‑1 (MQ=255)
                                                                                     cACCGACGCGCCGGACTAAAATGTATTCTTATGGTGGTTTTTCACCCTCGCCGTGTGGGAGGGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                 <  1:2497164/150‑1 (MQ=255)
                                                                                     cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCCACCCAAACCCCCACCCCACAAGGGGGGGGGGCCCCACGAGCAGaaaggga                                                                 >  2:2338962/1‑146 (MQ=21)
                                                                                       ccgccgCGACGGACGAAAAATTCTTCGTACGGTGGGTTTTCTCCCTCTCCCTGGGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                               <  2:1945176/150‑1 (MQ=38)
                                                                                       cccacgcgacgGACGAAAAAGATTTCGCAAGGTTGGTTTACTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                               <  2:1946079/147‑1 (MQ=255)
                                                                                       ccGACGCGACGGACGAAAAAGTATTCTTAAGGTTGGTTTTCTCCCACCCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                               <  2:2580007/150‑1 (MQ=255)
                                                                                       ccGACGCGACGGACGAAAAACTATTCGTATGGTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                               <  2:2960556/150‑1 (MQ=255)
                                                                                       ccGACCCTACGGACCAAAACTTATTCTTATGTTTGTTTTCCCCCCGCCCCCGGTGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                               <  2:530307/150‑1 (MQ=255)
                                                                                         gacgcgacgGACTAAAAAGTATTGGTATGGTTGTTTTTCTCCCTCTCCCTGGGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                             <  1:2771884/150‑1 (MQ=255)
                                                                                           cgcgccgGACGAAAACGTATTCGTACGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                           <  2:1636411/150‑1 (MQ=38)
                                                                                           cgcgacgGACGAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                           <  2:1636418/150‑1 (MQ=255)
                                                                                              gacggacgAACAAGTATTCGAAAGGTTGGTTTACTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGGGCACAAATTTTGAGAGatat                                                        <  2:2648782/150‑1 (MQ=38)
                                                                                                      aaaatttATTCTTATGGTTGGTTTTCTCCCTCGCCGTGTGGGTTCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACTGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                <  1:2062950/144‑1 (MQ=255)
                                                                                                      aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCACCGTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACTGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                <  1:2063462/150‑1 (MQ=255)
                                                                                                           tttttggTATGGTTGTTTTCCCCCCACCCCGGGGGGGTGTGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                           <  1:2268990/144‑1 (MQ=255)
                                                                                                           tgaTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                           <  2:1851735/148‑1 (MQ=255)
                                                                                                             aTTCTTATGTTTGGTTTTCGCCCACCCCCTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                         <  2:592727/150‑1 (MQ=255)
                                                                                                               tCGTATGGTTGTTTTTCCCCCTCCCCGTGGGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                       <  1:2707455/150‑1 (MQ=255)
                                                                                                                    tgtttgtttTTCTCCCTCTCCCTGGGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                  <  1:1711096/147‑1 (MQ=39)
                                                                                                                     ggttggttTTCTCCCTCGCCCTGTGGGTGTGGGTGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                 <  1:1892816/150‑1 (MQ=255)
                                                                                                                          gttTTCTCCCTCTCCGTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                            <  1:2161598/150‑1 (MQ=255)
                                                                                                                              tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTGGGCCCCCACCCCAACCCACACCCCCACGGGGGGGGGGGCCGAGCGAGGCCAAAATTTGTGAGAGATATTAATGTGTGTTGTAACACAAGCACAGaaaa                        >  2:2692559/1‑150 (MQ=12)
                                                                                                                                tccctctccgtGGGGTTGTGGGCGGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                      <  2:2502915/150‑1 (MQ=255)
                                                                                                                                cccctctccctGTGGGAGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                      <  2:1014640/149‑1 (MQ=255)
                                                                                                                                cccctcgccctGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                      <  2:2503106/149‑1 (MQ=255)
                                                                                                                                cccctcccccgGGGGTAGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                      <  2:1015332/149‑1 (MQ=255)
                                                                                                                                  cctctccgtGGGGTTGTGGGCGGGGGTGGGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                    <  1:561415/150‑1 (MQ=255)
                                                                                                                                        cgtgtgGTAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt              <  2:502308/148‑1 (MQ=255)
                                                                                                                                         ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCATAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg             >  1:766155/1‑150 (MQ=255)
                                                                                                                                              gttgcgggCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt        <  1:1864600/145‑1 (MQ=255)
                                                                                                                                              ggtgtgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt        <  1:1864605/145‑1 (MQ=255)
                                                                                                                                               tgtggggCCGGGGTTGGGGGCATCTGCTCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGTGGGGACCTATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt        <  1:826330/146‑1 (MQ=255)
                                                                                                                                                tgcgggCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa      <  1:2055802/147‑1 (MQ=255)
                                                                                                                                                    ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAATCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCGGAatat  >  2:819045/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAT  >  NC_000913/164420‑164715

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