Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 10.5% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,0680AC10.5% 460.5 / 7.5 232intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (94/107);  new base C (9/22);  total (103/129)
Fisher's exact test for biased strand distribution p-value = 8.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160919‑161217
                                                                                                                                                     |                                                                                                                                                       
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                        >  1:384670/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                        <  1:398353/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                        <  1:458216/150‑1 (MQ=255)
  gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                       >  1:1280050/1‑150 (MQ=255)
      ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGCCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcg                                                                                                                                                   >  1:4947278/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                              >  2:1078352/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                              >  2:1454580/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                             >  2:2063172/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGt                                                                                                                                            >  1:764267/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGt                                                                                                                                            >  1:2468029/1‑150 (MQ=255)
               agtgaTGAAAGGGCTTTATCCCAAACCTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGAGGGACGAAAACGTTTTGGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTtc                                                                                                                                          <  2:1003157/147‑1 (MQ=37)
                  gaTGAACGGGCGTTACCCCAAACTTTTCTGGCCGGACGACCCGGCAAATACTGCACCTCCGCGACGGACTAAAATTTTTTCGTATGGTGGTTTTGCTCCCTCCCCCTGTGGGTGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTtcacc                                                                                                                                       <  2:2589864/150‑1 (MQ=255)
                       aaGGGCGTTATCCCACACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                  <  2:1361093/150‑1 (MQ=255)
                       aaGGGCGTTATCCCAACCATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                  <  2:1993987/150‑1 (MQ=255)
                       aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                  >  2:3659929/1‑150 (MQ=255)
                         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                                <  2:3986020/150‑1 (MQ=255)
                             gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACAACCAAAc                                                                                                                            >  2:1850506/1‑150 (MQ=255)
                               tATCCAAAACATGTCGGGCCGGACGACCCGGCAATTACTTCACGGACGCGAGGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                                          <  2:4884438/150‑1 (MQ=37)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                         <  2:148741/150‑1 (MQ=255)
                                 cccaaaaCAGGTCGGGCCGGACCACCCGGCAAATACTGCACCGACGCGACGGACTAAAATGTATGCGAACGGTTGGTTTACCCCATCCCCCTGGGGTAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                        <  2:1956475/146‑1 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCACACCCGATCCcgct                                                                                                                       >  2:5290021/1‑150 (MQ=38)
                                  ccaaaaCATGTCGGGCCTGACGACCCGTCAAATACTGCCCCGACGCGACGTACGAAATTTTATTCGTAAGGTTGGTTTCCTCCCTCTCCCTGTGGTTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                       <  1:1624610/147‑1 (MQ=37)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                     <  2:3341283/150‑1 (MQ=255)
                                      aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCAACCtaaccctaacccacaacca                                                                                                                   >  2:3523774/1‑143 (MQ=38)
                                       aCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCAACCGCACCCAAACCCACTcccgc                                                                                                                  >  1:4442248/1‑148 (MQ=38)
                                        cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGACCCACACCCGAACCCACTaccgcc                                                                                                                 >  1:3481808/1‑147 (MQ=37)
                                           gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTTGACCCTCACCCGAAACCACTCCCGCAAgg                                                                                                              >  1:2102151/1‑150 (MQ=38)
                                           gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCCACCAAAACCTCACCCACaacg                                                                                                              >  2:2754286/1‑148 (MQ=38)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCACTCCCACACggg                                                                                                             >  1:3404057/1‑150 (MQ=255)
                                              tGGCCGGACGACCCGGCAAAGACTGCACCGCCGCGACGGACGTAAAAGTATTCGTAAGGTTGTTTTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                           <  2:3444158/150‑1 (MQ=38)
                                               ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCTCACCCCCAAGGgtga                                                                                                          >  1:3213470/1‑150 (MQ=37)
                                                 ccGGACGCCCCGTCAAATACGGCACCGACGCGACGGACGCAAAAGTATTCGTATGGTTGGTTTTCTCCCTCACCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                        <  2:3309430/150‑1 (MQ=38)
                                                    gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCGCAGCGCGCCCGTGCACCCACACCCACACCCACACCCGCCAGGGACAGGGGa                                                                                                     >  2:5178491/1‑150 (MQ=37)
                                                     acgacCCGGCAAATCCGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                    <  2:4115134/150‑1 (MQ=255)
                                                      cgccccGGCAAAAACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                   <  1:4469267/147‑1 (MQ=255)
                                                      cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCACCCGAAACCGCTACCACAAGGGAGTGGGGTcc                                                                                                   >  1:3403780/1‑150 (MQ=37)
                                                      cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCCCCCTCACCCGAAACCACTCCCTGACGGGAGCggggcgc                                                                                                   >  2:4401564/1‑147 (MQ=37)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCACCCACCCCCAAACCCACTCCCACAAGGGAGAGGGGAccc                                                                                                  >  1:1245811/1‑149 (MQ=37)
                                                       gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCCCCCCCAAACCCACTCCCGCAAGGGGGAGGGGGCCg                                                                                                  >  1:2208119/1‑150 (MQ=37)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGTGGACCGATc                                                                                               <  2:4690702/150‑1 (MQ=255)
                                                          ccGGAAAATCCTGCACCGACGCGCCGGACGAAAAAGTATTCTTAAGGTTGTTTTTATCCCTCTCCCTGTGGTTGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                               <  1:1374922/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                            <  2:1834380/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                            <  2:4487192/150‑1 (MQ=255)
                                                               aaaaACTGCACCCACGCGACGGACGAAAAAGTATTCGTACGTTTGGTTTTCTCCCACTCCCGGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                          <  1:402464/150‑1 (MQ=38)
                                                                aaTACGGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGATTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac                                                                                         <  1:3919822/150‑1 (MQ=255)
                                                                  tACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCAAAACCACTCCCTCCAGGGAGAGGGGACCGATCGCGCAC‑‑aa                                                                                     >  2:3188666/1‑150 (MQ=38)
                                                                  tACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGTCACCCACACCCAAACCCACACCCGCAAGGGGGCGGGGCCCCATAGAGcacca                                                                                       >  1:3835837/1‑148 (MQ=25)
                                                                    cggCACCGACGCGACGGACGAAACCGTATTCGTCCGGTTGGTGTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAAt                                                                                   <  2:5457134/148‑1 (MQ=255)
                                                                    cTGCACCTCCGCGAGGTACTAAAATTTATTGTTATGTTGGGTTTTCACCCGCCCCCTGGGGTTGTGGGCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAAt                                                                                   <  1:1067867/150‑1 (MQ=255)
                                                                    cTGCACCGACGCGAGGGACTAAAAATTATTGGTATGGTTGGTTTTCTCCCTCCCCCGGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAAt                                                                                   <  1:1572616/150‑1 (MQ=255)
                                                                     tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAAtt                                                                                  >  2:1694039/1‑150 (MQ=255)
                                                                     tGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGCACCCACACCCGAACCCTCTCCCACCCGGGGGGGGGGCCCGATCGTGCAC‑‑AAAtg                                                                                  >  2:4187379/1‑149 (MQ=25)
                                                                     ggCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAAtt                                                                                  <  2:4354251/149‑1 (MQ=255)
                                                                               cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCACACCCAAACCCACAACCACAAGGGAGTGGGGACCGAAAGAGCAC‑‑AAATTTTGAGAGata                                                                        >  1:4740931/1‑150 (MQ=37)
                                                                                gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCACCCTCACCCTAACCCTCTCCCTCAAGGGCGAGGGGACCCCTCCAGCAC‑‑AAATTTTGAGAGatat                                                                       >  2:3360586/1‑150 (MQ=38)
                                                                                gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCACACCCAAAACCACTACCACCAGGGAGAGGGGACCGATCGAGCAC‑‑CAAATTTGAGAGatat                                                                       >  1:4081985/1‑150 (MQ=37)
                                                                                 acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATAtc                                                                      <  2:3629294/150‑1 (MQ=255)
                                                                                  cggacgAAAATGTTTTCTTATGGTTGGTTTTCTCCAGCTCCGTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATAtct                                                                     <  1:2908391/150‑1 (MQ=255)
                                                                                  cggacgAAAAAGGATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATAtct                                                                     <  2:3619857/150‑1 (MQ=255)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCACACCCTAACCCTCTCCCGCCAGGGAGAGGGGACCGATCGAGCCC‑‑AAATTTTTGGAGATAtctt                                                                    >  1:3506559/1‑150 (MQ=37)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTCCCACAACGGCGAGGGGACCCATCGACCCC‑‑AAAATTTGATAGATAtatt                                                                    >  2:4682675/1‑150 (MQ=25)
                                                                                     acgAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATAtcttct                                                                  <  2:2983442/150‑1 (MQ=255)
                                                                                     acgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCAAACCCAAAAAGGGGAGGGGGCCCAAAGAGCACAAAAATTTTGAGAtataatat                                                                    >  2:2549394/1‑146 (MQ=16)
                                                                                      cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCAAACCCACCAGGGGGTGGGGGCCCCTAGAGCAC‑‑AAAATTTGTTAGATAtcttatt                                                                 >  2:408501/1‑150 (MQ=21)
                                                                                       gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATAtcttcttc                                                                <  1:3059139/150‑1 (MQ=255)
                                                                                         aaaaGTATTCTTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTg                                                              <  2:1183428/150‑1 (MQ=255)
                                                                                               aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCAAACCCACTCCCTACAGGGAGCGGGGACCGATAGAGCAC‑‑CAATGTTGAGAGCTATCATCTTCGGTCTGGt                                                        >  2:3151707/1‑150 (MQ=25)
                                                                                               aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGAACGTTCAACCTCAACCTAAACCAATACCTCAAGGGAGAGGGGAACGATAGAGCAC‑‑AAATTTTGAGACATATATTCATATGTCTTGt                                                        >  2:2859286/1‑150 (MQ=25)
                                                                                                 gcGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  1:3451290/149‑1 (MQ=255)
                                                                                                    tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACCACCGAAACCACTCCCTCACGGGTGAGGGGTCCCAAAGAGCAC‑‑AAAATTTTAGGGATATATTCTTATTTCTTGttcacg                                                   >  1:4840503/1‑145 (MQ=12)
                                                                                                      aGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCTCTCCCACAAGGGAGAGGGGACCGATCGAGCAC‑‑AAAATTTTAGAGATATATTCTTATGTCTTGTTAAagaa                                                 >  2:1910178/1‑150 (MQ=37)
                                                                                                       ggttgtttTCCTCCCTCGCCCTGGGGGAGTGGGCGGGGGGGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                                <  1:519342/150‑1 (MQ=255)
                                                                                                          tggttTTCTCCCTCGCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  2:1136944/150‑1 (MQ=255)
                                                                                                           ggttTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                            <  2:2031043/150‑1 (MQ=255)
                                                                                                           ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCACACCCAAAACCACACCCGCAAGGGGGAGGGGCGCGCTCGGGCAC‑‑AAAATTTTTGTGATATATTAATATGTCTTGTAAAAGAAGAACa                                            >  1:3611735/1‑150 (MQ=16)
                                                                                                           ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCTCACCCAAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGTGAGATATCTTATTCTGTCTTGTAACAGAAGTACa                                            >  2:4445516/1‑150 (MQ=38)
                                                                                                                  tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  2:794337/150‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  2:3461219/150‑1 (MQ=255)
                                                                                                                     ctcgcccgGGGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                  <  2:1501316/150‑1 (MQ=255)
                                                                                                                      tctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCAACCACACCCCAACCCACTACCACCAGGGGGGGGGGGCCCCACGAGGAC‑‑AAAATTTGTGAGGTATATTCTTGTGTGGTGGAACACGAGAACCGAAAATAgggg                                 >  2:235550/1‑149 (MQ=12)
                                                                                                                       ttccgtGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                <  1:3313339/149‑1 (MQ=255)
                                                                                                                       ctccctGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                <  2:1935406/150‑1 (MQ=255)
                                                                                                                       ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCAAACCCAAAACCACTCCCTCAAGGGCGAGGGGCCCGGTCCAGCAC‑‑AAATATTGAGAGATATATTATTATGTGTTTTTACAGCAGTACGGAACATCGGGcg                                >  2:3816992/1‑149 (MQ=19)
                                                                                                                          cctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  2:4611840/150‑1 (MQ=255)
                                                                                                                          ccggtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  2:2153571/147‑1 (MQ=255)
                                                                                                                           ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCACTACCACAAGGGGGTTGGGACCGATCGAGCAC‑‑AAATTTTGTGTGATTTCTTATTAGGTCGTGTGACATTAGAACAGAAAATCGGGGTTTTg                            >  2:1298193/1‑150 (MQ=255)
                                                                                                                            tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCGAAAACACACCCTCAAGGGGGTGGGGACCGTTCGAGCAC‑‑AAATTTTGGGAGATATATTATTATGTGTTGTAACAGAAGAACAGAAAATAgggggttttc                           >  2:4108730/1‑143 (MQ=255)
                                                                                                                               gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                        <  2:4931145/150‑1 (MQ=255)
                                                                                                                                ggAGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:1906952/150‑1 (MQ=255)
                                                                                                                                ggAGTGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:5014315/150‑1 (MQ=255)
                                                                                                                                  agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTCCCGCAAGGGAGAGGGGACCCATCGAGCAC‑‑AAAATTTGATAGATATATTATTATGTATTGTAAAAGAAGAACAGAAAATAGGGCTTTTTCGCCTGa                     >  2:132717/1‑150 (MQ=255)
                                                                                                                                    agGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCAAACCCTCTACCTCCATGGAGAGGGGCCCGATCGAGCAC‑‑AAATTTTGAGAGATAAATTATTATTTCTTGTAAAAGAATAACAGAAAAACGGGCTTTTGCGGCTGaaa                   >  2:4325234/1‑149 (MQ=255)
                                                                                                                                      ggCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                 <  1:4536251/150‑1 (MQ=255)
                                                                                                                                       gcgggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  1:2549394/147‑1 (MQ=255)
                                                                                                                                        ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg               <  1:2035290/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCAACCTAACCCGAAACCACTACCACAATGGAGAGGGGACCCATAGAGCAC‑‑AAATTTTTATAGATATATTATTATGTCGTGTAAAAGAAAAACAGAAAAACGGGGTTTTGAGCCAGAATATTggg             >  2:4970939/1‑148 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAAACCAAAACCACAAAGGAGAGGGGGCCGAACGAGCAC‑‑AAAATTTTTGATTTATAATATTAAGTAGTTTAAAAGAATAACAAAAAATAGGGCTGTTTCGCCGGaaaaatgag             >  2:1913834/1‑143 (MQ=255)
                                                                                                                                            ggTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  1:1910198/150‑1 (MQ=255)
                                                                                                                                                  cgCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:2791872/149‑1 (MQ=255)
                                                                                                                                                    cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc   <  2:4225970/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAC‑‑AAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160919‑161217

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: