Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,0510AT24.2% 217.4 / 21.1 128intergenic (+30/‑166)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (42/52);  new base T (1/29);  total (44/84)
Fisher's exact test for biased strand distribution p-value = 1.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 131 total aligned reads displayed.

TTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGCTTTTGCGCCTGAA  >  CP009273/160902‑161199
                                                                                                                                                     |                                                                                                                                                      
   tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                                     <  1:1103864/150‑1 (MQ=255)
       aaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                                 >  1:868334/1‑150 (MQ=255)
       aaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                                 >  1:868641/1‑150 (MQ=255)
            gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                            <  2:1690023/150‑1 (MQ=255)
             cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACAGGTCTGGCCGGACGACCCGGCAAAGACTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                           <  2:5208809/150‑1 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                          >  2:4510165/1‑150 (MQ=255)
               tACCGTGTGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                         <  1:499722/150‑1 (MQ=255)
               tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTAGCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                         <  1:5407612/150‑1 (MQ=255)
                aCCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTAACCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                                        >  2:256739/1‑150 (MQ=255)
                 ccGTGAGGTGCAAAATGAGATGAAAGGGCGTTATCCCAAACTTGTCTGGCCGGTCGACCCGGCAAATACTGCACCGACGCGACGGACGATCAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgggag                                                                                                                                       >  2:131282/1‑148 (MQ=39)
                 ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACAAGTCTGGCCGGTCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTTCGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgggga                                                                                                                                       >  2:131324/1‑148 (MQ=255)
                   gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                     >  1:662967/1‑150 (MQ=255)
                    tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATcc                                                                                                                                    >  1:2662021/1‑149 (MQ=255)
                    tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                    <  2:2804096/150‑1 (MQ=255)
                     gAGGTGAAAAAAGTTTTGAAAGGGCGTTATCCCAAAAATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAg                                                                                                                                   <  1:429318/150‑1 (MQ=255)
                      cggggCACAAAGAGTGTAAGGGGCGTTATCCCACAAAGTTCTGGCCGGACGACCCGGCACAAACGGCACCGACGCGACGTACTAAAAAGTATTCGTACGGTTGGTTTTACCCCTCACCCTGGGGGAGCGGGCCGGGGTGTGGGCCTCAgc                                                                                                                                  <  2:4146031/146‑1 (MQ=255)
                             aaaaGAGATGAATGGGCGTTATCCAAAACATGTCTGGCCGGACGACCGGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCGGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACg                                                                                                                           <  2:379618/150‑1 (MQ=255)
                               aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCggg                                                                                                                         >  2:305392/1‑148 (MQ=255)
                               aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                         >  2:2174846/1‑149 (MQ=255)
                               aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTAGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAGGtg                                                                                                                         >  2:304855/1‑149 (MQ=255)
                               aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTCACCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGtt                                                                                                                         >  2:2174501/1‑150 (MQ=255)
                                agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGtc                                                                                                                        >  1:4825507/1‑150 (MQ=255)
                                 gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTGTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGtgcg                                                                                                                       >  1:5394942/1‑147 (MQ=255)
                                   gaTTAATGGGGGGTATCCCCACCCTTTCTGGCGGTCCCACCGGGCAAATACGGCACCGCCGGGCCGGACGAAAAAGTTTTGGTATGGTGGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                     <  2:4621051/150‑1 (MQ=16)
                                    aTGAACGGGGTTTATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACTCCACCGACGCGACGTACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCATGTGGTATTGGGCGGGGTTTAGGGACTCAGCGCGCACGTtcaccc                                                                                                                    <  2:528466/150‑1 (MQ=255)
                                    aTGAACGGGCGTTATCCCAAACATGTCGGGCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAACGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                    <  2:1327459/150‑1 (MQ=38)
                                    aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccc                                                                                                                    >  1:2545041/1‑150 (MQ=255)
                                       aatgggCGTTATCCCTAACTTGTCTGGCCGGACGACCCGGCAAATACTGCCCCTACGCTCCGGACTAAAAAGTATTGGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTtcaccctca                                                                                                                 <  2:744752/147‑1 (MQ=255)
                                        acgggttttATCCCAAAAATGTGTGGCGGGACTACCCGGCAAAAACTCCCCCGACGCGACGTACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCGCCCTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                <  2:633747/143‑1 (MQ=255)
                                        acgggCTTTATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACTGCCCCGACGCGCCGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                <  2:3163016/148‑1 (MQ=38)
                                        aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccgcac                                                                                                                >  2:5045680/1‑150 (MQ=255)
                                        aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcactctcac                                                                                                                >  2:5045660/1‑150 (MQ=255)
                                           ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccacacccc                                                                                                             >  1:1813764/1‑149 (MQ=255)
                                           ggCGTTAACCAAAACATGTCGGGCCGGACGACCCGGCAAATAATTCCCGTAGGCTAGGGACGATAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTTTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                             <  2:4629137/150‑1 (MQ=25)
                                             cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGTGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                           <  2:3118287/150‑1 (MQ=255)
                                              gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCACCTAAc                                                                                                          >  2:1514378/1‑150 (MQ=255)
                                               ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCAACCaaaac                                                                                                         >  1:3856441/1‑148 (MQ=38)
                                                aagcccAAACAGGACGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAACTGTATTCGTAAGGTTGGTTTTCGCCCCCACCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                        <  2:4733679/147‑1 (MQ=255)
                                                    ccAAACATGTCTGGCCGGCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                    <  2:4000790/150‑1 (MQ=255)
                                                       aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCCCACCCAAAACCACAcccc                                                                                                 >  1:3479830/1‑149 (MQ=38)
                                                         cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCACACCCATAACCACTcccgcc                                                                                               >  2:4994340/1‑147 (MQ=37)
                                                          aTGTCTGGCCGGACGCCCCGGCAAAAACCGCACCGACGCGACGTACTAAAAATTATACGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGTTGAGGTCATCAGCTCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                              <  1:3577261/150‑1 (MQ=255)
                                                           tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAATTTTTCGTAAGGTTGGTTTACACCCACTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                             <  1:1302925/150‑1 (MQ=39)
                                                            gTCTGGCCGGACGACCCGCCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                            <  2:5442434/150‑1 (MQ=255)
                                                              cgggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                          <  2:555510/148‑1 (MQ=255)
                                                              cTGGCCGGACGCCCCGCAAAATACTCCACCGCCGGCACGGAATAAAAATTATTCTTACGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                          <  1:3747466/150‑1 (MQ=25)
                                                               tGGCCGTACGACCCGGAAAATACTGCACCGACGCGCCGGACTAAAATTTATTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                         <  2:1522129/150‑1 (MQ=37)
                                                               tGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAAGTATTCTTAAGTTTGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                         <  2:1522721/150‑1 (MQ=255)
                                                               gggCCGGACGACCCGGAAAATCATGCACCGACGCGAGTGACTAAACAGTATTCGTACGGTTGTTTTTCTCCCTCTCCCTGTGGTATTGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                         <  2:2000159/149‑1 (MQ=37)
                                                                 gcggcacgccCCGCAAAATACGGCCCCGACGCGACGTACTAAAATTTTTTCTTATGGTTGTTTTTCCCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                       <  1:3863633/145‑1 (MQ=255)
                                                                 gccgcccgacCCGCAAAAAACTGCACCGACGCCCCGGACGAAAATTTATTCGTACGGTTGTTTTTCCCCCTCCCCCTGTGGTAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                       <  1:3863005/144‑1 (MQ=255)
                                                                 gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCGCACCCCAAACCACCCCCGCAAgggtggg                                                                                       >  1:4911226/1‑146 (MQ=37)
                                                                  ccgccccacCCGCAAAATACTCCCCCCACGCCACGCACCAAAATTTATTCTTATGGTTGTTTTTCTCCCTCCCCTTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                      <  1:2168833/143‑1 (MQ=255)
                                                                   cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAACCCAATACCACCAGGGAGAggg                                                                                     >  2:3332018/1‑150 (MQ=37)
                                                                    ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCACTCCCGCCAgggggggggg                                                                                    >  2:4458536/1‑150 (MQ=37)
                                                                     gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCGAACCCTCACCCCGAAGGGGGAggggg                                                                                   >  2:3622248/1‑149 (MQ=37)
                                                                         acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGACCCGCACCCGAACCCACACCCACAAGGGGGAGGGGACCGa                                                                               >  1:4549992/1‑150 (MQ=37)
                                                                             ggCAAATACTGAACCTCCGCGACGGACGACACAGTAGACCTAAGGTTGGTTTTCTCCCTCTCCCGCTGGGCGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCACCCTCAAGGGAGAGGGGACCgcacga                                                                           <  2:730872/150‑6 (MQ=255)
                                                                                 aaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCACACCCACAAGGGGGAGGGGGCCGGTCGAGccc                                                                       >  1:3474165/1‑148 (MQ=25)
                                                                                    acggCACCGACGCGCCGGCCGAAAAAGTTTTCGTATGTTTGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACaaa                                                                   <  1:4877933/147‑1 (MQ=255)
                                                                                        cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAAACCACACCCACAAGGGTGAGGGGACCGATAGAGC‑CCAAAattt                                                               >  2:4254579/1‑150 (MQ=37)
                                                                                          ccGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAACCCGCACCCAAAACGGGGGGGGGACCCAGCCCGC‑ACAAAATTTgt                                                             >  2:481998/1‑149 (MQ=21)
                                                                                           cgacgctccgGACGAAAATTTATTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTATCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTgag                                                            <  2:365953/150‑1 (MQ=37)
                                                                                             acgcgacgGACGAAAATGTATTCTTATGGTTGTTTTTCCCCCTCTCCCTGTGGTTTTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTgagag                                                          <  1:4421275/150‑1 (MQ=255)
                                                                                                cgacgGACGAAAAAGTATTCGAAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGata                                                       <  2:256113/150‑1 (MQ=255)
                                                                                                 gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCAAACCCACAAGGGGGAGGGGGACCAAAGAGCAACAAATTTT‑ATAGatat                                                      >  1:828362/1‑150 (MQ=21)
                                                                                                   cggactAAAAAGTATTCTTATGTTTGGTTTACTCCCACGCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATAtct                                                    <  1:154495/150‑1 (MQ=255)
                                                                                                   cggacgAAAAAGTATTAGTAAGGTTGGTTTACGCCCTCTCCCTGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGTGGGGACCGATCGAGC‑ACAAATTTTGAGAGATAtct                                                    <  2:4420100/150‑1 (MQ=255)
                                                                                                          aaaaGTTTTCGTAAGGTGGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTg                                             <  2:4258134/150‑1 (MQ=255)
                                                                                                           aaaGTACCCTTATGGTTGGTTTCCTCCCACCCCCTGGGGGAGAGGCCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGt                                            <  1:307768/150‑1 (MQ=255)
                                                                                                             aGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCt                                          <  2:5370452/150‑1 (MQ=255)
                                                                                                             aGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCGAACCCACACCCACAAGGGAGAGGGGACCGATAGAGC‑CCCAAAATTGAGGGATATATTCGTATGTCt                                          >  2:2710878/1‑150 (MQ=21)
                                                                                                              gTATTCGACAGGTTGGTTTTCTCCCTCTCCCTGTGGGACAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCTCTCCCTAACCCTCTCCCTCAAGGGAGAGGGCACCGATCGAGC‑ACAAAATTTGAGAGATATCTTCTTCTGTCtt                                         <  2:2553621/150‑1 (MQ=38)
                                                                                                                aTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                       <  2:351163/150‑1 (MQ=255)
                                                                                                                 gtCGTATGGTTGGTTTTCTCCCTCTCCGGGGGGTAGAGGTCCGGGGTTAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                      <  1:18016/149‑1 (MQ=255)
                                                                                                                  tCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                     <  2:667210/150‑1 (MQ=255)
                                                                                                                   cttatgtttgtttTTCTCCCTCTCCCTGTGGTTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                    <  1:1486828/143‑1 (MQ=38)
                                                                                                                      acggagggttTTCCCCCACGCCCTGCGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTACACCCTCACCCTACCCCACTCCCTCACGGGAGAGGGGACCGATCGAGC‑ACAACTGTTGAGAGATCTCTTCTACTGTCTTGTAACaga                                 <  2:1175003/144‑1 (MQ=255)
                                                                                                                        ggttggttTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                               <  2:483223/150‑1 (MQ=255)
                                                                                                                           tggttTTCGCCCTCTCCGTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AAc                           <  2:3282317/150‑1 (MQ=255)
                                                                                                                           tggtgTTCACCCTCGCCCTGTGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AAc                           <  2:3282263/150‑1 (MQ=255)
                                                                                                                            ggttTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACa                          <  2:1160720/150‑1 (MQ=255)
                                                                                                                            ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCACACCCGAACCCTCTCCCTCAAGGGTGTGGGGACCGGTCGGGC‑GCAAATTTTGAGAGGTATGTTCTTGTGTCGTGTAACagtag‑acaa                          >  2:4076983/1‑147 (MQ=25)
                                                                                                                            ggttTTCTCCCTCACCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACa                          <  2:1160844/150‑1 (MQ=255)
                                                                                                                                  ctccctctccctGTGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAC‑AACAGAAAAt                    <  1:4895233/150‑1 (MQ=255)
                                                                                                                                   tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATc                   <  2:761876/150‑1 (MQ=255)
                                                                                                                                    ccctctccctGGGGGAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCg                  <  1:526437/150‑1 (MQ=255)
                                                                                                                                    ccccctccctGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCg                  <  2:3969965/150‑1 (MQ=255)
                                                                                                                                     cctctccctGGGGGTGAGGGCGGGTGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCgg                 <  1:314207/150‑1 (MQ=255)
                                                                                                                                       tctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGc               <  2:2802372/150‑1 (MQ=255)
                                                                                                                                       tcgccgggtgGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGc               <  1:2954762/143‑1 (MQ=255)
                                                                                                                                        ctccctGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGCt              <  1:771485/150‑1 (MQ=255)
                                                                                                                                          ccctGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGCttt            <  1:1039628/150‑1 (MQ=255)
                                                                                                                                          ccctGTGGGAGAGGTCGGGGGTTAGTGCATCAGCGCGCACTTTCACCCTCACCCTAACCCTCTCCATCATGGGAGAGGGGACCGATCGAGC‑ACAAATTTTTAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGCttt            <  2:4949516/150‑1 (MQ=255)
                                                                                                                                          ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGGTCACCCTCACCCTAACCCTCTCCCGCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGCttt            <  2:4457578/150‑1 (MQ=255)
                                                                                                                                          ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGCttt            <  2:4458449/150‑1 (MQ=255)
                                                                                                                                            ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCCAACCCACACCCCCCAGGGCGAGGGGACCCATAGAGC‑CCAAAATTTGTGAGGTATCTTATTATGTGTTTTTAAAAAAGAAACATAAAA‑CGGGCGTTTg          >  2:3756438/1‑150 (MQ=255)
                                                                                                                                             tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGCTTTTgc         <  2:4907972/150‑1 (MQ=255)
                                                                                                                                                  gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCTAACCCTCACCCACCAGGGAGAGGGGACCGCTCGCGC‑ACAAATTTTGAGAGATATATTCTTCTGTGTTGTAACAGAAT‑AACAGAAAATAGGGCTTTTTCGccgg    >  2:4509141/1‑148 (MQ=255)
                                                                                                                                                    gagGGCCGGGGTGAGGGCAACCGCGCGCACGTACACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATATTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGCTTTTGCGCCTGaa  <  2:3302862/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAG‑AACAGAAAATCGGGCTTTTGCGCCTGAA  >  CP009273/160902‑161199

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: ‑