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

AAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTG  >  NC_000913/164421‑164710
                                                                                                                                                      |                                                                                                                                              
aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                >  1:450068/1‑149 (MQ=255)
aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                <  1:698027/149‑1 (MQ=255)
aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                >  2:1050182/1‑149 (MQ=255)
   ggAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                             <  2:273188/149‑1 (MQ=255)
    gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgggg                                                                                                                                            <  1:354821/149‑1 (MQ=255)
     aGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                           >  1:1043256/1‑149 (MQ=255)
      gcttACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGTACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTATCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                                          <  2:1033978/146‑1 (MQ=255)
      gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                                          >  1:1513989/1‑149 (MQ=255)
      gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                                          <  1:813383/149‑1 (MQ=255)
      gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTg                                                                                                                                          >  2:28612/1‑149 (MQ=255)
       cGTACCGTGAGGTGCAAAAAGAGATGCAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                         <  2:1265746/149‑1 (MQ=255)
        gTACGGTGTGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  1:550717/149‑1 (MQ=255)
        gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  2:1975523/149‑1 (MQ=255)
        gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  1:1505776/149‑1 (MQ=255)
        gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  2:331198/149‑1 (MQ=255)
        gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  1:213792/149‑1 (MQ=255)
        gTACCGTGAGGTGAAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAAGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                        <  1:550730/149‑1 (MQ=39)
         tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                       >  2:1842447/1‑149 (MQ=255)
           ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                                     >  2:973531/1‑149 (MQ=255)
           ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGAAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                                     >  1:1468256/1‑149 (MQ=255)
                    gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                            >  2:1413109/1‑149 (MQ=255)
                    gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGCCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                            >  2:2014494/1‑149 (MQ=255)
                      aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATTCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                          >  2:320965/1‑149 (MQ=255)
                      aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                          >  1:1364871/1‑149 (MQ=255)
                        aacgacaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGAGGCGAGGGACTAAAAAGTATCCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        <  2:1249123/143‑1 (MQ=38)
                        aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACg                                                                                                                        >  2:669674/1‑149 (MQ=255)
                        aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        >  1:166398/1‑149 (MQ=255)
                        aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        >  1:921891/1‑149 (MQ=255)
                        aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        >  1:1306677/1‑149 (MQ=255)
                        aaaGAGATGAAAGGGCGTTATCCAAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                        >  1:719629/1‑149 (MQ=255)
                         aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                       >  1:47066/1‑149 (MQ=255)
                          agagATCCAAGGGCGTTAACCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATAGTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGtt                                                                                                                     <  2:2083344/149‑1 (MQ=38)
                            ataTGAATGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGCCGCGCCGGACGAAAACGTTTTCTTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtca                                                                                                                   <  2:166398/147‑1 (MQ=25)
                            agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGtgga                                                                                                                    >  2:1134151/1‑146 (MQ=255)
                               tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggcccc                                                                                                                 >  2:1785939/1‑143 (MQ=255)
                                gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccg                                                                                                               >  2:1475562/1‑148 (MQ=255)
                                 aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccac                                                                                                              >  1:2130270/1‑147 (MQ=255)
                                  atgggCTTTATCCAAAACTTGTCTGGCGGGACGCCCGGTCAAATACTCCACCGACGCGACGGACTAAAATTTATTGTAATGGTTGTTTTTCTCCCTCTCCCTGTGGTTTTGGGCGGGGTTGTGGGCATCAGCGCGCACGTtcaccctca                                                                                                             <  2:1865866/147‑1 (MQ=255)
                                    gggCGTTATCCCAAACATGTCTGGCCGGGCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccgcacc                                                                                                           >  1:1837934/1‑149 (MQ=255)
                                        gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCaaa                                                                                                       >  2:1184534/1‑149 (MQ=38)
                                           aTCCCAAACATTTCGGGCCGGACGCCCGGCAAAATCCTGCACCGACGCGACGGCCGAAAACTTTTACGTAAGGTTGTTTTTCCCCCCCTCCGTGTGGGAGCGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                    <  2:1382279/149‑1 (MQ=255)
                                           aTCCCAAACATGTCTGGCGGGACGACCCGTAAATTACTCCACCGCGGCGCCTGACTAAAAATTTTTCTTTTTTTTGTTTTTCCCCCCCCCGGTGGGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                    <  2:544305/149‑1 (MQ=255)
                                             cccAAACATGTCGGGCGGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAATTTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                  <  2:1937136/149‑1 (MQ=37)
                                               cAAACTTGTGGGGCGGGCCGACCCGCAAATTAGTGCACCGACGCGCGGGACTAAAAATTATTGGTATGTTTGTTTTTCCCCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                <  2:355189/149‑1 (MQ=255)
                                               cAAACAGGTCGGGCCGCACCCCCCGAAAAAAACTCCCCCGACGCGACGGACTAAAAATTATTCTTATGTTTGGTTTTCCCCCTCCCCGTGGGGTAGCGGGCGGGGGGGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                <  2:244395/149‑1 (MQ=255)
                                                aaccATGTCTGGCCGTACCCCCCGGCAAATACTGCACCGCCGCCCCGTACTAAAAATTATTCGTAAGTTGGGTTTTCGCCCACTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                               <  2:942493/146‑1 (MQ=255)
                                                 aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGGCCCccacccccacccacaacc                                                                                              >  1:1862891/1‑143 (MQ=37)
                                                     ggcggggCGGGCCGACCCGCAAACTACGCCACCGACGCGCCGCACTAAAATTTATTCTTACGGTTGGTTTCCCCCCTCCCCCTGGGGGCGGGGGCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                          <  1:2003351/144‑1 (MQ=255)
                                                       tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCGGCGCGCACGGTCACCCGCCCCCGGAGCCGCTCCCTCACgg                                                                                        >  2:1312122/1‑149 (MQ=37)
                                                        cTGGCCGGCCCCCCCGGCAAATACGGCCCCGCCGCGCCGGACTAAAAAGTATACGTATGGTTGGGTTGCTCCCTCCCCCTGGGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                                                                       <  2:1084543/149‑1 (MQ=255)
                                                           gCCGGACGCCCCGGAAACTACTGCACCGAGGCGCCGGACTAAAAAGTATTCTAACGGTTGGTTGACCCCCACCCCCGGGGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                    <  2:236371/149‑1 (MQ=255)
                                                           gCCGGACGCCCCGCAAAATACGGCCCCGCCGCGAGGCACTAAAAATTATTCGAAAGTTTGGTTTTCCCCCCCCCCCTGTGGGTGTGGGCGGGGTTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                    <  1:1506663/149‑1 (MQ=255)
                                                            ccGGACGACCCCGAAACTACTGCACCGACGCGCCGCACGAAAAAGTATCCGTAAGTTTGGTTTCCTCCCTCGCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                   <  2:1413724/149‑1 (MQ=37)
                                                                accacCCGGCAAATCCTCCCCCGACGCGACGCACGAAAAATTATTCTTACGGTTGTTTTACCCCCCCACCGTGTGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                               <  1:1568941/146‑1 (MQ=255)
                                                                   acCCGCAAAATACTCCCCCCCCGCGCCGTACAAAAAAGTATTCTTATGGTGGTTTTCCCCCCCCCCCCGGGGGGCGCGGGCCGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                            <  2:562866/149‑1 (MQ=255)
                                                                     cccgcaaaaTACGCCACCGCCGCGCCGGCCTAAAATTTTTTGGTATGGTTGTTTTTCCCCCTCTCCGTGGGGGTGGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                           <  1:1732796/143‑1 (MQ=255)
                                                                           aaaGACTGCACCGACGCGACGGACCAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCACTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCTCTACCTCAAGGGATATGGGAGCGATCGTGc                                                                     <  2:2056087/149‑1 (MQ=37)
                                                                              tACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCAACCAAACCCCCTCCCGCACGGGGGAGGGGGCCGCTCGAGcgca                                                                  >  1:2092573/1‑149 (MQ=25)
                                                                                    ccctccgcgacgTAATAAAATTTTTTCTTATGGTTGGTTTTCCCCCTCCCCGTGGGGGTGTGGGCGGGGTTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                            <  2:638020/144‑1 (MQ=255)
                                                                                    ccccacgcgccgGACTAAAATGTTTTCGTAAGTTTGTTTTTCTCCCTCCCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                            <  2:526290/145‑1 (MQ=255)
                                                                                    aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCCCCCCAACCCAAACCCACACCCGCAAGGGGGAGGGGCCCGCTCGCGCCCCAAtgtt                                                            >  2:2054880/1‑149 (MQ=21)
                                                                                    aCCGACGCGACGGAAGAAAATGTTTTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                            <  2:1864433/149‑1 (MQ=255)
                                                                                     ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCTCCCCCGAACCCACTCCCGCGAGGGGGCGGGGCCCCAGCCAGCGCACAATTTg                                                           >  2:1358849/1‑149 (MQ=21)
                                                                                      cgccgcgaggGACTAAAAATTTTTGTTATGGTTGGTTTTCCCCCTCCCCGGGGGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                          <  2:548003/146‑1 (MQ=255)
                                                                                          gcgacgGACTAAAAAGTTTTCGTTAGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                      <  2:1653948/149‑1 (MQ=255)
                                                                                             acggactAAAATGTTTTCTTATGTTTGGTTTTCGCCCACTCCGGGTGGGAGTGGGCCGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                   <  2:1364871/149‑1 (MQ=255)
                                                                                                    aaattttATTCGTATGGTTGGTTTTCACCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                            <  1:2014494/143‑1 (MQ=38)
                                                                                                    aaaaaGTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                            <  2:450068/149‑1 (MQ=38)
                                                                                                       aaGTATTGTTATGGTTGTTTTTCTCCCTCCCCGTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                         <  2:1306677/149‑1 (MQ=255)
                                                                                                           aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                     <  2:1169680/149‑1 (MQ=255)
                                                                                                            atCGAAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGTTGAGGGCATCAGCCCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                    <  1:2115819/148‑1 (MQ=255)
                                                                                                                  atggttggttTTCTCACACGCCCTGGGGCAGTGGGCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTGTGAGAGATATCTTCTACTGTCTTGTAACag                               <  2:1406442/147‑1 (MQ=255)
                                                                                                                   tggttgtttTTCTCCCACTCCCTTTGGTAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                              <  2:112222/148‑1 (MQ=255)
                                                                                                                   tggttggttATCCCCCTCACCCTGTGGGATCGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                              <  1:2085287/148‑1 (MQ=255)
                                                                                                                   agttttgttTTCCCCCTCGCCGTGGGGGTGGGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                              <  1:1031121/143‑1 (MQ=255)
                                                                                                                    ggttgtttTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                             <  1:1511225/149‑1 (MQ=255)
                                                                                                                    ggttggttTACCCCCTCCCCGTGTGGGCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                             <  2:1043256/149‑1 (MQ=255)
                                                                                                                           tttCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                      <  2:1779887/149‑1 (MQ=255)
                                                                                                                                ccctctccctGTGGGTGTGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                 <  1:1109585/149‑1 (MQ=255)
                                                                                                                                  ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCAACCAAAACCCCACCCGCAAGGGCGTGGGGTCCCCTCGAGCACCAATTTTGAGAGATATATTATTCAGTCTTTTTACAGAAAAAAAAAAAATAgg               >  2:1554526/1‑149 (MQ=16)
                                                                                                                                     tccctGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt            <  2:1989918/149‑1 (MQ=255)
                                                                                                                                       ccggtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt          <  2:1468256/146‑1 (MQ=255)
                                                                                                                                         tgtgggtgtgggCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg        <  1:1475562/149‑1 (MQ=255)
                                                                                                                                             ggtgtgggCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc    <  2:1797774/144‑1 (MQ=255)
                                                                                                                                             ggAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc    <  2:213210/149‑1 (MQ=255)
                                                                                                                                               agagGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCCACCGCCACCCAACCCACACCCCCAAGGGAGGGGGGGCCCATCGGGCACCAATTTTGGGAGATAAATTGTTATGGTGTGTAAGCGAAGAACAGAAAAAAGGGCgttttggcccg  >  2:684559/1‑142 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
AAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTG  >  NC_000913/164421‑164710

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: