Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,434,260 C→A 100% P111H (CCT→CAT)  tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,434,2600CA99.0% 655.3 / ‑2.0 198P111H (CCT→CAT) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (119/77);  total (119/79)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 211 total aligned reads displayed.

CGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCTCC  >  CP009273/4434114‑4434405
                                                                                                                                                  |                                                                                                                                                 
cgttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatc                                                                                                                                                 <  2:3366086/149‑1 (MQ=255)
 gttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcc                                                                                                                                                <  1:360252/149‑1 (MQ=255)
 gttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcc                                                                                                                                                >  2:2463358/1‑149 (MQ=255)
 gttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcc                                                                                                                                                <  2:3659882/149‑1 (MQ=255)
  ttcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGcatcct                                                                                                                                               >  1:1600097/1‑149 (MQ=255)
   tcgttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCtt                                                                                                                                              >  1:547122/1‑149 (MQ=255)
     gttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCt                                                                                                                                            <  2:2108205/149‑1 (MQ=255)
      ttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTg                                                                                                                                           <  1:3122475/149‑1 (MQ=255)
      ttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTg                                                                                                                                           >  1:550210/1‑149 (MQ=255)
      ttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTg                                                                                                                                           >  2:3929605/1‑149 (MQ=255)
      ttATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCGGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTg                                                                                                                                           <  1:1066854/149‑1 (MQ=255)
       tATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGa                                                                                                                                          >  2:2684378/1‑149 (MQ=255)
        aTGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGaa                                                                                                                                         <  2:378253/149‑1 (MQ=255)
        aTGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGaa                                                                                                                                         <  1:3880986/149‑1 (MQ=255)
         tGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAAc                                                                                                                                        <  2:2828046/149‑1 (MQ=255)
         tGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAAc                                                                                                                                        >  1:2924/1‑149 (MQ=255)
          gAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACa                                                                                                                                       >  1:2277210/1‑149 (MQ=255)
          gAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACa                                                                                                                                       <  1:3357249/149‑1 (MQ=255)
            gcagcCAGGCGTGGCGGTAAAAGCGGGCCATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAgg                                                                                                                                     >  2:1566579/1‑149 (MQ=255)
             cagcCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGAAAGTAAAAAAATGCATCCTTCTGAACAGGt                                                                                                                                    <  1:1021171/149‑1 (MQ=255)
              agcCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGtt                                                                                                                                   >  1:1660071/1‑149 (MQ=255)
                ccAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTga                                                                                                                                 >  1:3537977/1‑149 (MQ=255)
                   ggCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTgaaga                                                                                                                              >  1:2898739/1‑149 (MQ=255)
                    gCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTgaagaa                                                                                                                             <  1:3891917/149‑1 (MQ=255)
                        ggCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAgg                                                                                                                         <  1:4062360/149‑1 (MQ=255)
                             tAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGAt                                                                                                                    <  2:1793143/149‑1 (MQ=255)
                              aaaaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATa                                                                                                                   <  1:27170/149‑1 (MQ=255)
                              aaaaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATa                                                                                                                   <  2:3388902/149‑1 (MQ=255)
                              aaaaGCGGGCAATCGGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCGGGCGCTGAAAGACATTCAGGTCAACCTCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATa                                                                                                                   <  1:27079/149‑1 (MQ=255)
                                aaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGc                                                                                                                 >  1:1329940/1‑149 (MQ=255)
                                aaGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGc                                                                                                                 >  2:3176641/1‑149 (MQ=255)
                                    gggCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTc                                                                                                             <  1:1209905/149‑1 (MQ=255)
                                    gggCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTc                                                                                                             >  2:1571850/1‑149 (MQ=255)
                                         aTCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTg                                                                                                        >  2:366379/1‑149 (MQ=255)
                                            tGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAt                                                                                                     >  1:3603041/1‑149 (MQ=255)
                                            tGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAt                                                                                                     >  2:2537893/1‑149 (MQ=255)
                                             gCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAtc                                                                                                    <  1:1470441/149‑1 (MQ=255)
                                               aTCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGAtctc                                                                                                  <  2:550210/149‑1 (MQ=255)
                                                   ggCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCa                                                                                              >  2:2745786/1‑149 (MQ=255)
                                                   ggCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCa                                                                                              <  1:940926/149‑1 (MQ=255)
                                                    gCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCAc                                                                                             >  1:1183016/1‑149 (MQ=255)
                                                       gTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGcc                                                                                          <  1:3420331/149‑1 (MQ=255)
                                                         cGGGCTTGAGGGCCTGTGGAACAGCCTTGTTTGTATTAATGATCTGGCGCTTAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAATATAGCGGTCCGCTGGATCTCTCCACGCCGt                                                                                        <  1:2179229/149‑1 (MQ=255)
                                                              ttGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATccc                                                                                   >  1:2798915/1‑149 (MQ=255)
                                                                     ccTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCAccc                                                                            >  2:3561104/1‑149 (MQ=255)
                                                                     ccTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCAccc                                                                            <  2:1027860/149‑1 (MQ=255)
                                                                     ccTGTGGAACAGCAGTGTTTGTATTAATGAGCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCAccc                                                                            <  2:751826/149‑1 (MQ=255)
                                                                      cTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCt                                                                           >  2:1377660/1‑149 (MQ=255)
                                                                        gtgGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGa                                                                         <  2:3820699/149‑1 (MQ=255)
                                                                        gtgGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGa                                                                         <  1:2734996/149‑1 (MQ=255)
                                                                        gtgGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGa                                                                         <  2:2572913/149‑1 (MQ=255)
                                                                           gAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGacac                                                                      >  1:2582060/1‑149 (MQ=255)
                                                                           gAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGacac                                                                      <  2:615330/149‑1 (MQ=255)
                                                                               agcagTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGt                                                                  >  1:860939/1‑149 (MQ=255)
                                                                               agcagTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGt                                                                  <  1:3195578/149‑1 (MQ=255)
                                                                                gcagTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTg                                                                 >  1:2447735/1‑149 (MQ=255)
                                                                                gcagTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTg                                                                 <  2:2127763/149‑1 (MQ=255)
                                                                                    tgtTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCAc                                                             <  2:3940701/149‑1 (MQ=255)
                                                                                    tgtTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCAc                                                             <  2:1939999/149‑1 (MQ=255)
                                                                                       ttGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTgg                                                          >  2:1599354/1‑149 (MQ=255)
                                                                                       ttGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTgg                                                          >  2:1079352/1‑149 (MQ=255)
                                                                                          tATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACa                                                       >  1:1720608/1‑149 (MQ=255)
                                                                                          tATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACa                                                       >  2:2047425/1‑149 (MQ=255)
                                                                                           aTTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACaa                                                      >  1:3484628/1‑149 (MQ=255)
                                                                                               aTGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTc                                                  >  1:347416/1‑149 (MQ=255)
                                                                                                tGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCa                                                 >  2:1052708/1‑149 (MQ=255)
                                                                                                 gATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCaa                                                <  2:2993568/149‑1 (MQ=255)
                                                                                                 gATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCaa                                                <  2:2993627/149‑1 (MQ=255)
                                                                                                  aTCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAAc                                               >  2:1814486/1‑149 (MQ=255)
                                                                                                   tCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACa                                              <  1:2101547/149‑1 (MQ=255)
                                                                                                     tGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACAtt                                            <  1:843367/149‑1 (MQ=255)
                                                                                                          cTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGAt                                       <  2:476876/149‑1 (MQ=255)
                                                                                                           tGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGAtt                                      <  1:1922380/149‑1 (MQ=255)
                                                                                                            gAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtg                                     >  2:2357457/1‑149 (MQ=255)
                                                                                                            gAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtg                                     >  2:1376731/1‑149 (MQ=255)
                                                                                                            gAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtg                                     >  1:799120/1‑149 (MQ=255)
                                                                                                             aaaGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtga                                    <  2:3406681/149‑1 (MQ=255)
                                                                                                              aaGACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgat                                   >  1:2040485/1‑149 (MQ=255)
                                                                                                                gACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgatga                                 >  1:2601161/1‑149 (MQ=255)
                                                                                                                gACATTCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATtgatga                                 >  2:1269774/1‑149 (MQ=255)
                                                                                                                    ttCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGAcacc                             <  2:815509/149‑1 (MQ=255)
                                                                                                                    ttCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGAcacc                             <  2:1465811/149‑1 (MQ=255)
                                                                                                                    ttCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGAcacc                             <  2:1753403/149‑1 (MQ=255)
                                                                                                                    ttCAGGTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACAGTAAGATTGATGAcacc                             >  2:2964646/1‑149 (MQ=255)
                                                                                                                        ggTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACgg                         >  1:3952717/1‑149 (MQ=255)
                                                                                                                        ggTCAACATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCAAATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACgg                         >  2:2415163/1‑149 (MQ=255)
                                                                                                                            aaCATCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATc                     <  2:3070287/149‑1 (MQ=255)
                                                                                                                               aTCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGt                  >  2:2920928/1‑149 (MQ=255)
                                                                                                                                tCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTg                 >  2:411452/1‑149 (MQ=255)
                                                                                                                                tCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTg                 >  2:504512/1‑149 (MQ=255)
                                                                                                                                tCGACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTg                 >  1:113023/1‑149 (MQ=255)
                                                                                                                                  gACAGTAAAAAAATGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGAt               <  1:3646985/149‑1 (MQ=255)
                                                                                                                                        aaaaaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACtt         <  1:2684378/149‑1 (MQ=255)
                                                                                                                                        aaaaaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACtt         >  1:1281192/1‑149 (MQ=255)
                                                                                                                                         aaaaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTc        >  1:190440/1‑149 (MQ=255)
                                                                                                                                         aaaaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTc        >  1:2265466/1‑149 (MQ=255)
                                                                                                                                           aaaaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCAc      <  2:1727853/149‑1 (MQ=255)
                                                                                                                                             aaTGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCt    >  2:3407775/1‑149 (MQ=255)
                                                                                                                                               tGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCTcc  >  1:2621657/1‑149 (MQ=255)
                                                                                                                                               tGCATCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCTcc  >  1:2622419/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGTTCGTTATGAGCAGCCAGGCGTGGCGGTAAAAGCGGGCAATCTGCATCTGGCTGTCGGGCTTGAGTGCCTGTGGAACAGCAGTGTTTGTATTAATGATCTGGCGCTGAAAGACATTCAGGTCAACATCGACAGTAAAAAAATGCCTCCTTCTGAACAGGTTGAAGAAGAGGAAGATAGCGGTCCGCTGGATCTCTCCACGCCGTATCCCATCACCCTGACACGGGTGGCACTGGACAACGTCAACATTAAGATTGATGACACCACGGTATCGGTGATGGACTTCACCTCC  >  CP009273/4434114‑4434405

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: