Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,1000TA20.6% 505.1 / 23.9 173intergenic (+79/‑117)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (16/119);  new base A (33/2);  total (51/122)
Fisher's exact test for biased strand distribution p-value = 9.38e-21
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 181 total aligned reads displayed.

ATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  CP009273/160951‑161249
                                                                                                                                                     |                                                                                                                                                     
gtCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCGGCGCGCACGGTCGACCCCACCCAAAACct                                                                                                                                                       >  1:607731/2‑150 (MQ=38)
 tCCAAAACATGTCGGGCCGGACGACCCGGAAAATACTGCACCGCCGCGACGGACGAAAATTTATTCGTATGGTTTTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGGGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                                      <  1:4004902/150‑1 (MQ=37)
  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                                                     <  2:3353606/150‑1 (MQ=255)
   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCAAAACCtatc                                                                                                                                                    >  2:1664857/1‑150 (MQ=37)
    cAACCAGGTCTGGCGGTCCGACCGGGAAAATACTCCCCCGACGCGACGGACCAAAAATTATTCTTTAGGTTGGTTTCCGCCATCTCCCTGTGGTTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                                                   <  2:4511849/150‑1 (MQ=255)
     aaaCAGGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                                                  <  2:3487436/150‑1 (MQ=255)
       acctggcTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                                                                <  2:2556834/144‑1 (MQ=38)
        cATGTCTGGCCGGACGACCCGGAAAAAAATGCACCGAGGCGAGGGACGAAAAAGTTTTCTAAAGGTTGGTTTTCACCCTCTCCCTGTGGTAGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                                                               <  2:1218819/150‑1 (MQ=255)
          tGTCGGGCCGGACGACCCGGAAAATACTGCCCCGAGGCGCCGGACTAAAAAGTATTCTTAAGGTTGGTTTACTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                             <  2:1003330/150‑1 (MQ=255)
               ggccggccGACCCGGAAAATACGGCACCGACGCGCCGGACGAAAAATTATTCGTAGGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                                        <  1:2923725/150‑1 (MQ=255)
                gCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTATACGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                       <  2:4062897/150‑1 (MQ=255)
                 gcGGACTCCCCGGAAAATACTGCACCGACGGGACGGACGAAAAAGTATTCGTACGGTTGGTTTGCTCCCTCGCCGTGTGGGAGAGGGCCGGGGTGAGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                                                      <  1:3507999/149‑1 (MQ=255)
                 ccGGCCGCCCCGCAAAATCCGCCACCGACGCGACGGACTAAAATTTTTTCTTATGGTTGTTTTTCCCCCTCCCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                                                      <  1:3366234/150‑1 (MQ=255)
                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCCCCCTAACCCACTCCCTCAAGGGCGAgg                                                                                                                                      >  2:994691/1‑150 (MQ=38)
                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCCAAACCAATCCCAAAACGGagagtg                                                                                                                                     >  2:4579917/1‑148 (MQ=37)
                   ggACGACCCGGCAAATACTGCCCCGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTGTCCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggg                                                                                                                                    <  2:151162/150‑1 (MQ=38)
                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGAACCCGCACCCAAAACCACACCCGCAAGGGGGAgggg                                                                                                                                    >  2:416090/1‑150 (MQ=37)
                    tccgacCCGGCAAATCCTGCCCCGACGCGACGGACTAAAATGTATGCGTACGGTTGTTTTTCTCCCCCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                                                   <  1:1595983/148‑1 (MQ=255)
                      cgacCCGGAAAATCCGCCCCCGACGCGCCGGACGAAAAAGTTTTCGTAAGGTTGGTTTCCACCCCCTCCCGGTGGGTTTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                                                 <  1:3015153/150‑1 (MQ=255)
                       gccccGTCAAATACTGCCCGGACGCGACGGACGAAATAGTATTCGTACGGTTGGTTTTCTCCACCCCCCTGTGGGAGAGGGCGGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                                                                <  2:4276402/148‑1 (MQ=255)
                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCAAACCCAAAACCACTACCTCAAGGGGGAGGGGAccc                                                                                                                                >  2:1343257/1‑149 (MQ=37)
                         cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                                                              <  2:3577229/150‑1 (MQ=255)
                         cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCCAACCCAATACCACAAGGGAGAGGGGAcccct                                                                                                                              >  2:2490120/1‑147 (MQ=37)
                          ccGGCAAAACCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                                             <  2:3232292/150‑1 (MQ=255)
                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCACACCCACAACCACAAAGGGGAGGGGGCCGATCg                                                                                                                            >  1:1088010/1‑150 (MQ=37)
                            ggCAAATACTGCACCGCCGCCCCGGACGAAACAGTATTCGTAAGGTTGGTTTTCGCCCACGCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                                                           <  2:27568/150‑1 (MQ=38)
                            ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGTCCCCCACCCAAAACCAAAACCACCAGGGGGAGGGGCCCGAACgg                                                                                                                           >  1:2246302/1‑149 (MQ=21)
                               acaTACTGCACCGACGCGACGGACGAAAACGTATACGTAAGGTTGGTTGCCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                                        <  2:929589/148‑1 (MQ=255)
                               aaaTCCTGCCCCGACGCGACGGACTAAAAATTATTCGTAAGGTTGGTTTTCGCCCTCTCCCTGGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                                        <  1:1148650/150‑1 (MQ=255)
                                  tccTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTACTCCCTCTCCCGGTGGGAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaa                                                                                                                     <  2:1279573/148‑1 (MQ=255)
                                   aCTGCACCGACGCGAGGGACTAAAAAGTATTCGTACGGTTGTTTTTCTCCCTCTCCGTGGGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                                                                                                    <  2:4610880/150‑1 (MQ=38)
                                   aCTCCACCGACGCGACGTACGAAAATGTATTGGTATGGTTGTTTTGCTCCCTCTCCCTGTGGGTGGGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACana                                                                                                                    <  1:4644061/150‑1 (MQ=255)
                                     tGCACCGACGCGACGGCCGAAAAAGTATACGTACGGTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                                                  <  2:2321338/150‑1 (MQ=255)
                                     tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCTCCACCGAAACCTCTCCCACAACGGGGGGGGGGCCGGTCGAGCCCCAAtt                                                                                                                  >  2:2137462/1‑150 (MQ=25)
                                     gtcACCGACGCGACGGACTAAAATTTATTCTTAAGTTTGGTTTACTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                                                  <  1:669842/148‑1 (MQ=255)
                                     ggcccccacgcgacgGACGAAAAAGTATGCGTCCGGTTGGTTTACTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                                                  <  2:4435450/143‑1 (MQ=38)
                                     ggCACCGACGCGACGTACGAAAAAGTATTCGAAAGGTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                                                  <  2:2122823/149‑1 (MQ=255)
                                        aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                                               <  1:3147900/150‑1 (MQ=255)
                                          cgacgcgacgGACTAAAAAGTATTCGAAAGGTGGGTTTTCTCCCTCCCCGTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                                             <  2:868962/150‑1 (MQ=38)
                                              gcgacgGACGAAAACGTATACGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                                                         <  2:1974242/150‑1 (MQ=255)
                                                  cggactAAAAATTAGTCGTCGGGTTGGTTTTCGCCCTCTCCGTGTGGGTGTGGGCCGGGGTGTGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                                                     <  1:2968332/150‑1 (MQ=255)
                                                  cggacgAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                                                     <  2:426810/150‑1 (MQ=255)
                                                  cggacgAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                                                     <  2:2282278/150‑1 (MQ=255)
                                                  cggacgAAAAAGTTTTCGTATGGTTGGTTTTCTCCCCCCCCCGGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                                                     <  1:3019962/150‑1 (MQ=39)
                                                  cggacgAAAAAGATTTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                                                     <  2:1973625/150‑1 (MQ=38)
                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCCCACCCCAAACCAAAACCGCCAGGGGGATGGGGCCGAAAGAGCCCAAATTTTGTGAGttataat                                                                                                    >  1:2009384/1‑147 (MQ=18)
                                                    gactAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCACCCTGTGGTAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                                                   <  2:571035/150‑1 (MQ=255)
                                                      cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                                                 <  2:386800/150‑1 (MQ=255)
                                                          aaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCCCCCCCCAACCCACAACCACAACGGGGATGGGGCCGCACGCGCACAAAAGTTGAGAGGTAtattaaaattt                                                                                             >  1:4526201/1‑143 (MQ=12)
                                                          aaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                                             <  2:1689001/150‑1 (MQ=255)
                                                             gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                          <  1:3555738/150‑1 (MQ=255)
                                                             gTATTCGTAAGGTTGGTTGGCTCCCTCGCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                          <  2:2956057/150‑1 (MQ=255)
                                                              tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                                         <  2:3934613/150‑1 (MQ=255)
                                                               aTTCGTACGGTGGGTTTTCTCCCTCTCCCTGGGGGTGTGGGCCGGGGGGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                                                        <  1:4229720/150‑1 (MQ=255)
                                                               aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCCCCCAAAACCACAACCTCAAGGGAGAGGGGCCCCATCGAGCACAAAATTTTAGAGATATATTATTCAGTATTGt                                                                                        >  2:3930806/1‑150 (MQ=25)
                                                                     aaGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACACCCACCAGGGCGCGGGGGCCGCGCGCGCACCAATTTTGTGTGAGATATTATTAGGTGTTGTGACaga                                                                                  >  1:4157445/1‑150 (MQ=12)
                                                                               ttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                                                        <  2:3805928/150‑1 (MQ=255)
                                                                                   ccatctccttTTGGGAGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACACTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                                    <  1:1289914/147‑1 (MQ=255)
                                                                                        tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGCGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                                                               <  2:586702/150‑1 (MQ=255)
                                                                                         ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                                              <  2:4336310/150‑1 (MQ=255)
                                                                                          cctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                                             <  1:3804224/150‑1 (MQ=255)
                                                                                           ctGTGGGTGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                                                            <  1:4580438/150‑1 (MQ=255)
                                                                                           cggtgGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                                                            <  1:4790707/148‑1 (MQ=255)
                                                                                            tgtgGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCTCACCCGAACCCGCTCCCACAAGGGAGAGGGGCCCGATCGAGCCCAAAATTTGTGTGCTATATTCTTCTGGCTTGTTAAAGGAGCACACGAAATAGGGCgttttg                                                           >  2:3521711/1‑148 (MQ=255)
                                                                                             gtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                                                          <  2:49875/150‑1 (MQ=255)
                                                                                                ggAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:3810206/150‑1 (MQ=255)
                                                                                                ggAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:1653307/150‑1 (MQ=255)
                                                                                                ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:3810289/150‑1 (MQ=255)
                                                                                                  agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                                                     <  1:3929401/150‑1 (MQ=255)
                                                                                                      ggCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                                                 <  1:2196859/150‑1 (MQ=255)
                                                                                                            ggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                           <  1:2567275/150‑1 (MQ=255)
                                                                                                            ggTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATCTCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                           <  2:2793401/150‑1 (MQ=255)
                                                                                                             gTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTACCACAAGGGAGAGGGGACCCAAAGAGCACAAATTTTGAGTTATATATTAATATGTATTGTAACACAAGAACACAAAAACGGGCGTTTTCGCCTGAAAATTTCGgag                                          >  2:3738996/1‑150 (MQ=255)
                                                                                                             gTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACACCCCCAACGGCGAGGGGCCCCCTCGAGCACCAAATTTGAGAGAGATATTCTTATGTGATTTAACACAAGAAAATAAAAAAGGGCGTTTGAGCCGGAATATTggggcg                                          >  2:150126/1‑148 (MQ=255)
                                                                                                               gAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCGTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa                                        <  1:2773378/150‑1 (MQ=255)
                                                                                                                aGGGCATCAGCGCGCACGTGCACCCACACCCAAACCCACAACCACAACGGTGAGGGGGCCGATCGAGCCCAAAATTTGTGTGTTATAATCTTCTGTCTTGTTAAAGAAGAACAGAAAAAACGGCTTTTTCGCCGGAATATTGAGGAGgaa                                       >  2:4085462/1‑150 (MQ=255)
                                                                                                                 gggCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCACAAAGGAGAGGGGACCCATAGACCACAAAATTTTAGATATATATTATTATGTCATGTAACAGAATAAAAGAAAAAACGGCTTTTTCGCCAGAAAAATTCGGAGaaaa                                      >  2:3221493/1‑150 (MQ=255)
                                                                                                                   gCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                                    <  2:1641703/150‑1 (MQ=255)
                                                                                                                     aTCAGCGCGCACGTTCACCCACACCCAAAACCACACCCACAAAGGAGAGGGGACCCAAAGAGCACAAAATTTTAGAGATATAATATTATGTATTGTAAAAAAAAAAAAAAAAAAAGGGCATTTTCGCCAGAAAATTTCGGAGAAAAAACa                                  >  2:176057/1‑150 (MQ=255)
                                                                                                                       cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATg                                <  1:2804512/150‑1 (MQ=255)
                                                                                                                        aGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATgg                               <  2:639437/150‑1 (MQ=255)
                                                                                                                        aGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATgg                               <  2:4098143/150‑1 (MQ=255)
                                                                                                                           gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                            <  1:1218338/150‑1 (MQ=255)
                                                                                                                           gcgcACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                            <  2:4689537/150‑1 (MQ=255)
                                                                                                                                cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt                       <  1:4679430/150‑1 (MQ=255)
                                                                                                                                    caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc                   <  2:1600243/150‑1 (MQ=255)
                                                                                                                                     accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                  <  2:2103226/150‑1 (MQ=255)
                                                                                                                                     accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGACCAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                  >  2:2050437/1‑150 (MQ=255)
                                                                                                                                     accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                  >  1:4516800/1‑150 (MQ=255)
                                                                                                                                          caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATGTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc             <  2:1187656/150‑1 (MQ=255)
                                                                                                                                           accctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGcc            >  2:3200177/1‑150 (MQ=255)
                                                                                                                                           accctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGcc            <  2:3563763/150‑1 (MQ=255)
                                                                                                                                                aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       <  2:152664/150‑1 (MQ=255)
                                                                                                                                                aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       >  2:1521093/1‑150 (MQ=255)
                                                                                                                                                 aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg      <  2:1797357/150‑1 (MQ=255)
                                                                                                                                                  cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     >  2:2859119/1‑150 (MQ=255)
                                                                                                                                                  cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     >  2:1724200/1‑150 (MQ=255)
                                                                                                                                                  cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     >  1:78236/1‑150 (MQ=255)
                                                                                                                                                    ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg   <  2:4723954/150‑1 (MQ=255)
                                                                                                                                                     tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  1:1233118/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  CP009273/160951‑161249

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: