Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,041,662 G→T 10.5% F580L (TTC→TTA gfcD ← putative O‑antigen capsule production periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,041,6620GT10.5% 820.5 / 46.8 225F580L (TTC→TTAgfcDputative O‑antigen capsule production periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (116/85);  new base T (13/11);  total (129/96)
Fisher's exact test for biased strand distribution p-value = 8.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

CTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACATCGTTTCCAG  >  CP009273/1041513‑1041811
                                                                                                                                                     |                                                                                                                                                     
 tGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGa                                                                                                                                                      >  1:3726969/1‑150 (MQ=255)
   aTTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGaaa                                                                                                                                                    >  2:2844944/1‑150 (MQ=255)
    ttGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAAc                                                                                                                                                   >  2:1828748/1‑150 (MQ=255)
     tGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCCTTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAAcc                                                                                                                                                  >  2:3072454/1‑150 (MQ=255)
     tGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAAcc                                                                                                                                                  <  1:1671011/150‑1 (MQ=255)
      gccgccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCt                                                                                                                                                 >  1:1051873/1‑150 (MQ=255)
      gccgccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCt                                                                                                                                                 <  1:3813041/150‑1 (MQ=255)
      gccgccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGACACTAAACCt                                                                                                                                                 >  1:2984771/1‑150 (MQ=255)
        cgccCATACCCCTACCGCCACGCCGCTGTTAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGt                                                                                                                                               <  1:1743377/150‑1 (MQ=255)
        cgccCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGt                                                                                                                                               >  1:1341142/1‑150 (MQ=255)
            cATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGc                                                                                                                                           <  1:1557256/150‑1 (MQ=255)
            cATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGc                                                                                                                                           >  2:1430353/1‑150 (MQ=255)
            cATACCCCGACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGc                                                                                                                                           <  2:3081862/150‑1 (MQ=255)
             aTACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCg                                                                                                                                          >  2:1204018/1‑150 (MQ=255)
              tACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGt                                                                                                                                         >  1:2706196/1‑150 (MQ=255)
               aCCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTg                                                                                                                                        <  1:4069078/150‑1 (MQ=255)
               aCCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTg                                                                                                                                        >  2:1316300/1‑150 (MQ=255)
                ccccTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTgg                                                                                                                                       <  2:2609692/150‑1 (MQ=255)
                 cccTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGa                                                                                                                                      <  1:730996/150‑1 (MQ=255)
                 cccTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGa                                                                                                                                      <  1:731003/150‑1 (MQ=255)
                 cccTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGa                                                                                                                                      >  1:2932668/1‑150 (MQ=255)
                     aCCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGaataa                                                                                                                                  >  1:4042966/1‑150 (MQ=255)
                      ccGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGaataat                                                                                                                                 <  1:4019021/150‑1 (MQ=255)
                      ccGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGaataat                                                                                                                                 <  1:2926542/150‑1 (MQ=255)
                         ccACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCgg                                                                                                                              >  2:850036/1‑150 (MQ=255)
                           aCGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTg                                                                                                                            <  1:1636462/150‑1 (MQ=255)
                           aCGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTg                                                                                                                            <  2:3662813/150‑1 (MQ=255)
                                      aaaGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGcatcct                                                                                                                 >  2:1931241/1‑148 (MQ=255)
                                      aaaGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGcatcat                                                                                                                 >  2:77098/1‑150 (MQ=255)
                                      aaaGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGcatcat                                                                                                                 >  2:1040158/1‑150 (MQ=255)
                                       aaGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCAtg                                                                                                                <  1:982025/150‑1 (MQ=255)
                                            tttGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTc                                                                                                           <  1:3106672/150‑1 (MQ=255)
                                            tttGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTc                                                                                                           >  1:1257526/1‑150 (MQ=255)
                                              tGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGACGGCGGGTTCCAGTAAGCCGTCACGAAACCTGGAGGCGTGGAATAATCGGTGAAACGCATCATGTTGTccc                                                                                                         >  2:2081284/1‑150 (MQ=255)
                                                gCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAg                                                                                                       <  2:1985629/150‑1 (MQ=255)
                                                 cGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGt                                                                                                      <  2:2932663/150‑1 (MQ=255)
                                                 cGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGt                                                                                                      >  1:1059759/1‑150 (MQ=255)
                                                   aCGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCa                                                                                                    <  2:1383357/150‑1 (MQ=255)
                                                    cGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCAc                                                                                                   >  1:421077/1‑150 (MQ=255)
                                                      tCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGt                                                                                                 >  2:2298375/1‑150 (MQ=255)
                                                       cGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGtt                                                                                                >  2:41815/1‑150 (MQ=255)
                                                       cGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGtt                                                                                                >  1:462820/1‑150 (MQ=255)
                                                         aTCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGc                                                                                              >  1:3547501/1‑150 (MQ=255)
                                                         aTCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGc                                                                                              >  1:2664175/1‑150 (MQ=255)
                                                           cGTTGCGCCTTTATCTTTTGCCATATATTGCCCAACGCAAAGTTTCATCAGTACGCCATAGAGCGTCGGCGGGTTCAAGTAAGCCGTCACGAAACCTGTTGGCGTTGAATAATCTGTGAAACGATTCAAGTTTTCAAAGTCACGTTGAtt                                                                                            >  1:2112589/1‑150 (MQ=255)
                                                           cGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCtt                                                                                            >  2:4148090/1‑150 (MQ=255)
                                                           cGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCtt                                                                                            <  2:4275861/150‑1 (MQ=255)
                                                             ttGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCa                                                                                          <  1:1909967/150‑1 (MQ=255)
                                                             ttGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAACAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCa                                                                                          <  1:3330598/150‑1 (MQ=255)
                                                             gtGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCa                                                                                          <  2:1604858/149‑1 (MQ=255)
                                                              tGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCAc                                                                                         >  1:2555168/1‑150 (MQ=255)
                                                                     ttATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTa                                                                                  <  1:2563204/150‑1 (MQ=255)
                                                                      tATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTaa                                                                                 <  1:2811200/150‑1 (MQ=255)
                                                                        tCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACg                                                                               <  1:1111533/150‑1 (MQ=255)
                                                                         cTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGt                                                                              >  2:2788157/1‑150 (MQ=255)
                                                                          ttttGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTc                                                                             <  1:3090132/150‑1 (MQ=255)
                                                                            ttGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCca                                                                           <  2:3609404/150‑1 (MQ=255)
                                                                              gCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCcaca                                                                         >  2:605884/1‑150 (MQ=255)
                                                                              gCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCcaca                                                                         <  1:3631676/150‑1 (MQ=255)
                                                                                 aGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACAccc                                                                      <  2:361650/150‑1 (MQ=255)
                                                                                  gATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCa                                                                     >  2:3865870/1‑150 (MQ=255)
                                                                                      ttGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAgcgc                                                                 >  1:715198/1‑150 (MQ=255)
                                                                                        gCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGccc                                                               >  1:4480792/1‑150 (MQ=255)
                                                                                        gCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGccc                                                               <  2:176289/150‑1 (MQ=255)
                                                                                        gCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGccc                                                               >  1:2472713/1‑150 (MQ=255)
                                                                                           cAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCcagc                                                            >  1:2244993/1‑150 (MQ=255)
                                                                                            aaCGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCcagca                                                           >  2:4397138/1‑150 (MQ=255)
                                                                                               gCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGc                                                        >  1:2404995/1‑150 (MQ=255)
                                                                                               gCTAAGTTTCATCAGTACGCCAATGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGc                                                        >  2:101139/1‑150 (MQ=255)
                                                                                                   aGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACTAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCt                                                    <  1:2338045/150‑1 (MQ=255)
                                                                                                   aGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCt                                                    >  1:3522331/1‑150 (MQ=255)
                                                                                                    gTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTa                                                   <  2:1780647/150‑1 (MQ=255)
                                                                                                     tttCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAg                                                  >  2:3334484/1‑150 (MQ=255)
                                                                                                          tcaGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGc                                             >  1:3053600/1‑150 (MQ=255)
                                                                                                            aGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGa                                           >  2:2514214/1‑150 (MQ=255)
                                                                                                             gTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGAt                                          <  1:4317909/150‑1 (MQ=255)
                                                                                                                cGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATaaa                                       >  2:4194744/1‑150 (MQ=255)
                                                                                                                 gCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAag                                      >  2:2159099/1‑150 (MQ=255)
                                                                                                                  ccATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAagc                                     >  1:3149296/1‑150 (MQ=255)
                                                                                                                      tGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCt                                 <  1:629360/150‑1 (MQ=255)
                                                                                                                       gAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTc                                <  1:1432420/150‑1 (MQ=255)
                                                                                                                       gAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTc                                <  2:2045789/150‑1 (MQ=255)
                                                                                                                          cGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGa                             >  1:1498088/1‑150 (MQ=255)
                                                                                                                            tCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAAc                           <  1:407706/150‑1 (MQ=255)
                                                                                                                               gcggGTTCCAGTAAGCCGTCACTAAAACTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCATTCAAGTTGCTTCACGTAGTTAACTTCCACACCCAGCGCCCAGCAGTCATCTAGCGGGCGATAAAGCAGCTCTGAACCGa                        >  2:3479554/1‑150 (MQ=255)
                                                                                                                                cggGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGac                       >  2:1078210/1‑150 (MQ=255)
                                                                                                                                 gggTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGaca                      >  2:2426728/1‑150 (MQ=255)
                                                                                                                                 gggTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGaca                      <  2:3157460/150‑1 (MQ=255)
                                                                                                                                  ggTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGacac                     >  2:352863/1‑150 (MQ=255)
                                                                                                                                    ttCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCg                   <  2:4123281/150‑1 (MQ=255)
                                                                                                                                     tCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCgg                  <  2:2216276/150‑1 (MQ=255)
                                                                                                                                       cAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCg                >  1:1669121/1‑150 (MQ=255)
                                                                                                                                        aGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGt               >  2:4598652/1‑150 (MQ=255)
                                                                                                                                          tAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTAc             >  1:585442/1‑150 (MQ=255)
                                                                                                                                          tAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTAc             >  1:372013/1‑150 (MQ=255)
                                                                                                                                           aaGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACa            >  1:3144698/1‑150 (MQ=255)
                                                                                                                                           aaGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACa            <  1:1163267/150‑1 (MQ=255)
                                                                                                                                            aGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACAt           >  1:1566741/1‑150 (MQ=255)
                                                                                                                                                gTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACATCGtt       >  2:2242625/1‑150 (MQ=255)
                                                                                                                                                     gAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACATCGTTTCCAg  >  1:27732/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGATTGCCGCCCATACCCCTACCGCCACGCCGCTGTCAAAGCGTTTGGCGACGTCGATCGTTGCGCCTTTATCTTTTGCCAGATATTGCCCAACGCTAAGTTTCATCAGTACGCCATTGAGCGTCGGCGGGTTCCAGTAAGCCGTCACGAAACCTGTTGGCGTGGAATAATCGGTGAAACGCATCATGTTGTCCCAGTCACGTTGCTTCACGTAGTTAACGTCCACACCCAGCGCCCAGCAGGCATCTAGCGGGCGATAAAGCAGCTCGGAACCGACACCGGCGTACATCGTTTCCAG  >  CP009273/1041513‑1041811

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