Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,303 A→T 9.6% E717D (GAA→GAT sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273557,3030AT 9.6% 996.1 / 74.8 314E717D (GAA→GATsfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base A (135/148);  new base T (13/17);  total (148/166)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 322 total aligned reads displayed.

GGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCA  >  CP009273/557154‑557452
                                                                                                                                                     |                                                                                                                                                     
gggTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaa                                                                                                                                                       <  1:4145629/150‑1 (MQ=255)
gggTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaa                                                                                                                                                       >  2:2395902/1‑150 (MQ=255)
gggTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaa                                                                                                                                                       <  1:1521894/150‑1 (MQ=255)
     taACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCa                                                                                                                                                  <  1:3387748/150‑1 (MQ=255)
     taACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCa                                                                                                                                                  >  1:2887374/1‑150 (MQ=255)
       aCTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCaaa                                                                                                                                                >  1:3179234/1‑150 (MQ=255)
       aCTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCaaa                                                                                                                                                <  2:4083116/150‑1 (MQ=255)
       aCTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCaaa                                                                                                                                                >  1:4059600/1‑150 (MQ=255)
         tACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACt                                                                                                                                              <  2:3152536/150‑1 (MQ=255)
         tACGACCGCGATCAACATGACGTGAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACt                                                                                                                                              <  2:2013731/150‑1 (MQ=255)
          aCGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTg                                                                                                                                             >  2:2980077/1‑150 (MQ=255)
          aCGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTg                                                                                                                                             >  1:1566557/1‑150 (MQ=255)
          aCGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAAAGGCATAACGCTGAGCCAGCATTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGAAGGCGGTGTACGCAAACAAAATCAAAcag                                                                                                                                             >  2:2978468/1‑148 (MQ=255)
           cGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAGCTgg                                                                                                                                            <  2:1537582/150‑1 (MQ=255)
            gACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGc                                                                                                                                           <  1:2601626/150‑1 (MQ=255)
             aCCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCa                                                                                                                                          <  2:859024/150‑1 (MQ=255)
              ccGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCAt                                                                                                                                         <  2:3471655/150‑1 (MQ=255)
               cgcgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtt                                                                                                                                        >  1:350307/1‑150 (MQ=255)
               cgcgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtt                                                                                                                                        >  2:3011451/1‑150 (MQ=255)
                   aTCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTaa                                                                                                                                    <  2:716294/150‑1 (MQ=255)
                      aaCATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCg                                                                                                                                 >  2:1115257/1‑150 (MQ=255)
                      aaCATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCg                                                                                                                                 >  2:2679340/1‑150 (MQ=255)
                      aaCATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCg                                                                                                                                 <  1:1426674/150‑1 (MQ=255)
                       aCATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGa                                                                                                                                >  1:2528517/1‑150 (MQ=255)
                        cATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGAc                                                                                                                               >  1:209777/1‑150 (MQ=255)
                           gACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTgg                                                                                                                            <  1:3105401/150‑1 (MQ=255)
                            aCGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGgc                                                                                                                           <  1:4138588/150‑1 (MQ=255)
                             cGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGgcg                                                                                                                          <  1:2011379/150‑1 (MQ=255)
                                tAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCgg                                                                                                                       >  1:972681/1‑150 (MQ=255)
                                          cTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTg                                                                                                             >  2:3799043/1‑150 (MQ=255)
                                          cTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTg                                                                                                             >  2:3292332/1‑150 (MQ=255)
                                           tGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGa                                                                                                            <  1:3139202/150‑1 (MQ=255)
                                             tCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATg                                                                                                          >  1:3660035/1‑150 (MQ=255)
                                             tCTGGCGGTGTGGTCGGGCAGGAAAATGGCATAACGCTGAGCCAGCCTTGAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATg                                                                                                          <  2:4201606/150‑1 (MQ=255)
                                                          tCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACgg                                                                                             >  1:298969/1‑150 (MQ=255)
                                                          tCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACgg                                                                                             <  1:1574337/150‑1 (MQ=255)
                                                           cGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGt                                                                                            <  1:1087276/150‑1 (MQ=255)
                                                            gggCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGtt                                                                                           <  1:96925/150‑1 (MQ=255)
                                                             ggCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGttt                                                                                          >  1:3712883/1‑150 (MQ=255)
                                                             ggCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGttt                                                                                          <  2:2578516/150‑1 (MQ=255)
                                                              gCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTa                                                                                         >  1:3423243/1‑150 (MQ=255)
                                                              gCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTa                                                                                         >  1:2301844/1‑150 (MQ=255)
                                                                aTGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATc                                                                                       >  2:2615858/1‑150 (MQ=255)
                                                                   aaaaTGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGt                                                                                    >  2:1167829/1‑150 (MQ=255)
                                                                      aTGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGtata                                                                                 >  1:3233458/1‑150 (MQ=255)
                                                                          cATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCg                                                                             <  1:4378772/150‑1 (MQ=255)
                                                                           aTAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGt                                                                            >  2:2034456/1‑150 (MQ=255)
                                                                            tAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTa                                                                           <  1:1144308/150‑1 (MQ=255)
                                                                             aaCGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTAt                                                                          >  1:3766380/1‑150 (MQ=255)
                                                                              aCGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATc                                                                         <  1:1295453/150‑1 (MQ=255)
                                                                                  tGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCgcgc                                                                     <  1:4469781/150‑1 (MQ=255)
                                                                                  tGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCgcgc                                                                     <  1:2233728/150‑1 (MQ=255)
                                                                                   gAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCt                                                                    >  1:2543566/1‑150 (MQ=255)
                                                                                   gAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGAATCGCGCt                                                                    >  2:3998435/1‑150 (MQ=255)
                                                                                    agccagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCtt                                                                   >  1:2292589/1‑150 (MQ=255)
                                                                                     gccagccTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTg                                                                  <  1:353941/150‑1 (MQ=255)
                                                                                       cagccTTTAGGGGATACCAATGTTTTGATTAAAGCTCCTGGCGCAGGCGGTGTACGCATTGAAAATCAATCTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGAt                                                                >  1:1206889/1‑150 (MQ=255)
                                                                                       cagccTTTAGGGGATACCAATGTTTTGATTAAAGCTCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGAt                                                                >  1:1207303/1‑150 (MQ=255)
                                                                                           cTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCa                                                            <  2:43106/150‑1 (MQ=255)
                                                                                           cTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCa                                                            <  2:3394918/150‑1 (MQ=255)
                                                                                           cTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCa                                                            <  1:1618168/150‑1 (MQ=255)
                                                                                            tttAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCaa                                                           <  2:3317146/150‑1 (MQ=255)
                                                                                             ttAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAAt                                                          <  2:59687/150‑1 (MQ=255)
                                                                                             ttAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAAt                                                          <  1:805542/150‑1 (MQ=255)
                                                                                             ttAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAAt                                                          <  2:3892064/150‑1 (MQ=255)
                                                                                              tAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATa                                                         <  2:2097907/150‑1 (MQ=255)
                                                                                                 gggATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGa                                                      >  2:202132/1‑150 (MQ=255)
                                                                                                    aTACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATgg                                                   >  2:1709744/1‑150 (MQ=255)
                                                                                                    aTACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATgg                                                   >  2:1853406/1‑150 (MQ=255)
                                                                                                    aTACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATgg                                                   <  1:1209913/150‑1 (MQ=255)
                                                                                                     tACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATggg                                                  <  1:4363593/150‑1 (MQ=255)
                                                                                                     tACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATggg                                                  >  1:1466944/1‑150 (MQ=255)
                                                                                                      ctcAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACTCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATgggg                                                 <  2:4408082/148‑1 (MQ=255)
                                                                                                      aCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATgggg                                                 <  2:2692397/150‑1 (MQ=255)
                                                                                                       ccAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGa                                                <  1:3327564/150‑1 (MQ=255)
                                                                                                       ccAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGa                                                <  1:3591240/150‑1 (MQ=255)
                                                                                                        cAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGaa                                               >  2:922030/1‑150 (MQ=255)
                                                                                                         aaTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAAt                                              >  1:1437775/1‑150 (MQ=255)
                                                                                                         aaTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAAt                                              >  2:1464319/1‑150 (MQ=255)
                                                                                                          aTGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAAtt                                             >  2:4360519/1‑150 (MQ=255)
                                                                                                          aTGTTTTGATTAAAGCGCCTGGCGCAGGAGGTGTACGCATTGTAAATCAATCTGTCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAAtt                                             <  1:2686229/150‑1 (MQ=255)
                                                                                                           tGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTc                                            >  1:660626/1‑150 (MQ=255)
                                                                                                             ttttGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCa                                          <  1:40999/150‑1 (MQ=255)
                                                                                                                tGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCg                                       <  1:390750/150‑1 (MQ=255)
                                                                                                                  aTTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGAt                                     <  1:2454871/150‑1 (MQ=255)
                                                                                                                     aaaGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGtt                                  >  1:2417106/1‑150 (MQ=255)
                                                                                                                       aGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGa                                >  1:762035/1‑150 (MQ=255)
                                                                                                                          gcCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaaa                             >  2:2187283/1‑150 (MQ=255)
                                                                                                                           ccTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaaaa                            <  2:2743056/150‑1 (MQ=255)
                                                                                                                            cTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaaaaa                           >  1:563072/1‑150 (MQ=255)
                                                                                                                             tGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGaaaaaaa                          >  2:2647727/1‑150 (MQ=255)
                                                                                                                               gcgcAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAata                        <  1:4358274/150‑1 (MQ=255)
                                                                                                                                     gCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTagca                  <  1:548047/150‑1 (MQ=255)
                                                                                                                                         tgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGt              >  1:4041458/1‑150 (MQ=255)
                                                                                                                                         tgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGt              <  2:1437775/150‑1 (MQ=255)
                                                                                                                                         tgtACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGt              >  1:1789042/1‑150 (MQ=255)
                                                                                                                                           tACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTtg            <  2:1661242/150‑1 (MQ=255)
                                                                                                                                             ataatagAAAATCAAACTTGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCAACGTTTTATCTGTATTACCGTATCGCGATTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAATAATATTTGCTGCGTtgtg          <  2:375498/144‑1 (MQ=255)
                                                                                                                                                  tGATAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGAc     <  2:3179234/150‑1 (MQ=255)
                                                                                                                                                  tGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGACTTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGAc     >  2:1824916/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTATGCCACGGTTTATCGGTATAACCGTATCGCGCTTGATACCAATACGATGGGGAATTCCATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCA  >  CP009273/557154‑557452

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: