Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,529,483 A→T 13.4% intergenic (‑165/‑1291) nanC ← / → fimB N‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,529,4830AT13.4% 575.9 / 68.1 187intergenic (‑165/‑1291)nanC/fimBN‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (83/79);  new base T (10/15);  total (93/94)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 192 total aligned reads displayed.

AGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGATCTCATTT  >  CP009273/4529336‑4529631
                                                                                                                                                   |                                                                                                                                                    
aGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTcaca                                                                                                                                                    >  2:2605382/1‑150 (MQ=255)
 ggATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACat                                                                                                                                                   >  2:2504428/1‑150 (MQ=255)
 ggATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACat                                                                                                                                                   >  2:3358162/1‑150 (MQ=255)
  gATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACatt                                                                                                                                                  <  2:2568688/150‑1 (MQ=255)
   aTACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACatta                                                                                                                                                 <  2:1779867/150‑1 (MQ=255)
   aTACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACatta                                                                                                                                                 <  2:4129917/150‑1 (MQ=255)
      cAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTAttt                                                                                                                                              >  2:960065/1‑150 (MQ=255)
        gAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCACCCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTtaa                                                                                                                                            >  1:243567/1‑150 (MQ=255)
          aaGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCTCATTATTtaaat                                                                                                                                          >  2:2340494/1‑150 (MQ=255)
          aaGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTtaaat                                                                                                                                          >  1:4048512/1‑150 (MQ=255)
           aGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCTCATTATTtaaata                                                                                                                                         >  1:2819911/1‑150 (MQ=255)
               aTACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACa                                                                                                                                     >  1:450380/1‑150 (MQ=255)
                 aCAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACaga                                                                                                                                   >  1:1488689/1‑150 (MQ=255)
                   aaaGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCACCCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAgaga                                                                                                                                 >  1:218834/1‑150 (MQ=255)
                   aaaGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAgaga                                                                                                                                 >  2:2335985/1‑150 (MQ=255)
                   aaaGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAgaga                                                                                                                                 >  2:2337389/1‑150 (MQ=255)
                   aaaGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAgaga                                                                                                                                 >  2:1713523/1‑150 (MQ=255)
                     aGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATc                                                                                                                               <  1:4087900/150‑1 (MQ=255)
                     aGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATc                                                                                                                               <  1:1035219/150‑1 (MQ=255)
                           cATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGAtt                                                                                                                         >  1:519687/1‑150 (MQ=255)
                           cATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGAtt                                                                                                                         >  2:588813/1‑150 (MQ=255)
                           cATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGAtt                                                                                                                         >  2:4438/1‑150 (MQ=255)
                           cATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATATCTTTTCACATTATTTAAATAAACAGAGATCCAGAtt                                                                                                                         >  2:4077/1‑150 (MQ=255)
                             tAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTaa                                                                                                                       >  1:1075022/1‑150 (MQ=255)
                             tAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTaa                                                                                                                       >  1:1819026/1‑150 (MQ=255)
                              aaaCAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTaaa                                                                                                                      >  1:3596828/1‑150 (MQ=255)
                              aaaCAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTaaa                                                                                                                      >  1:179795/1‑150 (MQ=255)
                                    ttAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCTCATTATTTAAATAAACAGAGATCCAGATTAAATACCTg                                                                                                                <  1:3619422/150‑1 (MQ=255)
                                      aGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAg                                                                                                              <  2:53385/150‑1 (MQ=255)
                                        cATTTATGTTGTTTAATATTTTTTTGCCAGTCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGta                                                                                                            <  2:2672501/150‑1 (MQ=255)
                                         aTTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGtat                                                                                                           >  2:1263141/1‑150 (MQ=255)
                                         aTTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGtat                                                                                                           >  2:1666385/1‑150 (MQ=255)
                                          tttATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGtata                                                                                                          <  1:4225418/150‑1 (MQ=255)
                                             aTGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATaaaa                                                                                                       >  2:3299302/1‑150 (MQ=255)
                                              tgttgtGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAAt                                                                                                      <  2:1371220/150‑1 (MQ=255)
                                               gttgtGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAAtc                                                                                                     >  2:2321744/1‑150 (MQ=255)
                                                ttgtGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAAtct                                                                                                    <  2:854098/150‑1 (MQ=255)
                                                ttgtGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAAtct                                                                                                    <  2:3701792/150‑1 (MQ=255)
                                                   tgtAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCtctt                                                                                                 <  2:3741761/150‑1 (MQ=255)
                                                    gtAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCTCATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCtcttc                                                                                                >  2:4037489/1‑150 (MQ=255)
                                                    gtAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCtcttc                                                                                                <  1:688656/150‑1 (MQ=255)
                                                       atatTTTTTTGCCAGGCTTATAGTGTCTTTGGCCACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGa                                                                                             >  1:2279935/1‑150 (MQ=255)
                                                       atatTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGa                                                                                             >  2:931186/1‑150 (MQ=255)
                                                        tatTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGAt                                                                                            <  2:1946558/150‑1 (MQ=255)
                                                           ttttttGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGtt                                                                                         <  2:547731/150‑1 (MQ=255)
                                                           ttttgtGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGtt                                                                                         <  2:3877579/150‑1 (MQ=255)
                                                             ttttGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTa                                                                                       <  2:68857/150‑1 (MQ=255)
                                                             ttttGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTa                                                                                       <  1:2265212/150‑1 (MQ=255)
                                                               ttGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAAt                                                                                     <  2:4069825/150‑1 (MQ=255)
                                                                tGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAAtt                                                                                    <  2:2937277/150‑1 (MQ=255)
                                                                 gCCAGGCGTATAGGGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTGGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTg                                                                                   <  2:4220671/150‑1 (MQ=255)
                                                                   cAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCTCATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGAt                                                                                 <  2:2213355/150‑1 (MQ=255)
                                                                   cAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGAt                                                                                 >  2:1068189/1‑150 (MQ=255)
                                                                        ttATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGaa                                                                            <  2:1144583/150‑1 (MQ=255)
                                                                              tgtCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCTCATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCg                                                                      <  1:1165872/150‑1 (MQ=255)
                                                                              tgtCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCg                                                                      <  2:3945022/150‑1 (MQ=255)
                                                                                tCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTa                                                                    <  1:3676332/150‑1 (MQ=255)
                                                                                  tttGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAg                                                                  <  2:2172370/150‑1 (MQ=255)
                                                                                    tGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCt                                                                <  2:179795/150‑1 (MQ=255)
                                                                                     ggCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCta                                                               <  1:1614211/150‑1 (MQ=255)
                                                                                       cAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAATCtata                                                             >  1:2043491/1‑150 (MQ=255)
                                                                                           cGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCAc                                                         <  2:2043491/150‑1 (MQ=255)
                                                                                            ggTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACt                                                        <  1:3884256/150‑1 (MQ=255)
                                                                                             gTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACtt                                                       <  1:1046982/150‑1 (MQ=255)
                                                                                             gTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACtt                                                       <  1:1133345/150‑1 (MQ=255)
                                                                                             gTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACtt                                                       >  1:3551895/1‑150 (MQ=255)
                                                                                                     aTTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATc                                               >  1:835678/1‑150 (MQ=255)
                                                                                                     aTTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATc                                               >  1:835573/1‑150 (MQ=255)
                                                                                                      ttttATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCTCATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATcc                                              >  1:3346374/1‑150 (MQ=255)
                                                                                                       tttATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCa                                             <  2:1770544/150‑1 (MQ=255)
                                                                                                       tttATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCa                                             >  2:2383594/1‑150 (MQ=255)
                                                                                                        ttATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCAt                                            >  1:2592463/1‑150 (MQ=255)
                                                                                                         tataTTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCTCATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCAtt                                           >  1:3672502/1‑150 (MQ=255)
                                                                                                          atatTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCAttt                                          <  1:2650365/150‑1 (MQ=255)
                                                                                                           tatTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCAtttt                                         >  1:3821316/1‑150 (MQ=255)
                                                                                                                 tttGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCag                                   >  2:3026219/1‑150 (MQ=255)
                                                                                                                   tGTATAAGGTCTCCTGTGAAAAATCTCTATTCTCATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCagag                                 <  2:2862921/150‑1 (MQ=255)
                                                                                                                     tataAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTc                               >  2:1851422/1‑150 (MQ=255)
                                                                                                                       taAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCAt                             >  2:243482/1‑150 (MQ=255)
                                                                                                                        aaGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGGTCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATa                            >  1:2713098/1‑150 (MQ=255)
                                                                                                                        aaGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATa                            <  1:420550/150‑1 (MQ=255)
                                                                                                                         aGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATaa                           <  1:2605382/150‑1 (MQ=255)
                                                                                                                         aGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATaa                           <  1:4003761/150‑1 (MQ=255)
                                                                                                                           gTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATaaaa                         >  2:1784783/1‑150 (MQ=255)
                                                                                                                               ccTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGa                     >  2:3508274/1‑150 (MQ=255)
                                                                                                                                cTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                    <  2:3919652/150‑1 (MQ=255)
                                                                                                                                cTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGaa                    <  2:206818/150‑1 (MQ=255)
                                                                                                                                 tgtgAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAAt                   <  1:2230769/150‑1 (MQ=255)
                                                                                                                                      aaaaTCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTg              <  2:98663/150‑1 (MQ=255)
                                                                                                                                       aaaTCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCTTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATATAATgtgtgt             >  2:3206014/1‑150 (MQ=255)
                                                                                                                                       aaaTCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGt             >  2:2118552/1‑150 (MQ=255)
                                                                                                                                       aaaTCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGt             <  2:1066679/150‑1 (MQ=255)
                                                                                                                                       aaaTCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGt             >  1:3138042/1‑150 (MQ=255)
                                                                                                                                         aTCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTg           >  2:4161563/1‑150 (MQ=255)
                                                                                                                                           ctctTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCAGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAt         <  2:388230/150‑1 (MQ=255)
                                                                                                                                             ctTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtct       <  2:2279935/150‑1 (MQ=255)
                                                                                                                                             ctTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtct       >  2:491406/1‑150 (MQ=255)
                                                                                                                                              ttttCTCATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGAtctc      >  1:4116522/1‑150 (MQ=255)
                                                                                                                                                 tCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGTGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGATCTCAtt   <  1:1713523/150‑1 (MQ=255)
                                                                                                                                                  cacaTTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGATCTCAttt  <  2:4006130/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAATCTCTTTTCACATTATTTAAATAAACAGAGATCCAGATTAAATACCTGAGTATAAAATCTCTTCTGATGTTTAATTGATTTGAATGTTCGTAAGCTATATCACTTACTCAATCCATTTTACCCAGAGTCATAAAATAGAATGTCTGTTGATCTCATTT  >  CP009273/4529336‑4529631

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: