Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,798,582 C→A 12.9% Q243H (CAG→CAT waaP ← kinase that phosphorylates core heptose of lipopolysaccharide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,798,5820CA12.9% 1010.4 / 78.6 318Q243H (CAG→CATwaaPkinase that phosphorylates core heptose of lipopolysaccharide
Reads supporting (aligned to +/- strand):  ref base C (151/126);  new base A (27/14);  total (178/140)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 319 total aligned reads displayed.

TAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGC  >  CP009273/3798435‑3798728
                                                                                                                                                   |                                                                                                                                                  
   ttCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGt                                                                                                                                                 >  2:1670848/1‑148 (MQ=255)
    tCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGtt                                                                                                                                                >  1:3412709/1‑148 (MQ=255)
      aTAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCa                                                                                                                                              <  1:1977348/148‑1 (MQ=255)
        aaTCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAag                                                                                                                                            >  1:3335523/1‑148 (MQ=255)
         aTCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTGTAATACGCTCGGCCTTAACGCTAGCCATGTTCAATAAATACTGTTCAaga                                                                                                                                           >  1:2367205/1‑148 (MQ=255)
          tCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAagag                                                                                                                                          <  2:3511349/148‑1 (MQ=255)
             gCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACa                                                                                                                                       <  1:3473034/148‑1 (MQ=255)
             gCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACa                                                                                                                                       >  2:66928/1‑148 (MQ=255)
               aaCCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAat                                                                                                                                     <  1:3002658/148‑1 (MQ=255)
               aaCCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAat                                                                                                                                     >  1:999200/1‑148 (MQ=255)
                 cccAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAatat                                                                                                                                   >  1:1228082/1‑148 (MQ=255)
                 cccAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAatat                                                                                                                                   <  1:332747/148‑1 (MQ=255)
                   cAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATAttt                                                                                                                                 <  2:704888/148‑1 (MQ=255)
                   cAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATAttt                                                                                                                                 <  2:302127/148‑1 (MQ=255)
                      tttATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGc                                                                                                                              >  2:2678336/1‑148 (MQ=255)
                      tttATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGc                                                                                                                              >  2:451944/1‑148 (MQ=255)
                          taaaataaaTAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTaa                                                                                                                          <  1:1377274/148‑1 (MQ=255)
                          taaaataaaTAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTaa                                                                                                                          <  2:1719333/148‑1 (MQ=255)
                           aaaataaaTAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTaaa                                                                                                                         >  1:3505352/1‑148 (MQ=255)
                            aaataaatAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAg                                                                                                                        >  2:3223631/1‑148 (MQ=255)
                                aaataGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCAt                                                                                                                    >  2:2734120/1‑148 (MQ=255)
                                  ataGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATcc                                                                                                                  <  1:1687601/148‑1 (MQ=255)
                                   taGTCTTAATAAAGTGAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAAGTCTGTTCAAGAGACAATATTTTGCGTATAGGCATccc                                                                                                                 <  1:3032256/148‑1 (MQ=255)
                                    aGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCa                                                                                                                >  1:2906911/1‑148 (MQ=255)
                                      tCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCaaa                                                                                                              >  2:3290251/1‑148 (MQ=255)
                                      tCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCaaa                                                                                                              <  1:962239/148‑1 (MQ=255)
                                       cATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCaaaa                                                                                                             >  2:2692279/1‑148 (MQ=255)
                                           ataaAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATaaa                                                                                                         <  2:2083121/148‑1 (MQ=255)
                                                gttagttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTc                                                                                                    <  1:1879747/148‑1 (MQ=255)
                                                 ttagttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCa                                                                                                   >  1:2782276/1‑148 (MQ=255)
                                                  tagttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCAt                                                                                                  <  2:1538161/148‑1 (MQ=255)
                                                  tagttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCAt                                                                                                  <  1:3202149/148‑1 (MQ=255)
                                                   agttCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATa                                                                                                 >  1:2326055/1‑148 (MQ=255)
                                                         aGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGc                                                                                           >  2:2922389/1‑148 (MQ=255)
                                                          gTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGcc                                                                                          <  2:505795/148‑1 (MQ=255)
                                                            aCATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAg                                                                                        >  2:869788/1‑148 (MQ=255)
                                                              aTACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGat                                                                                      <  1:2655482/148‑1 (MQ=255)
                                                              aTACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGat                                                                                      <  2:1409962/148‑1 (MQ=255)
                                                               tACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGata                                                                                     >  1:1077930/1‑148 (MQ=255)
                                                               tACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTATTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGAAAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGata                                                                                     <  1:1301979/148‑1 (MQ=255)
                                                                aCTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGatat                                                                                    >  1:1778066/1‑148 (MQ=255)
                                                                 cTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATAtc                                                                                   <  2:908135/148‑1 (MQ=255)
                                                                 cTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATAtc                                                                                   >  1:359863/1‑148 (MQ=255)
                                                                 cTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATAtc                                                                                   >  1:2241844/1‑148 (MQ=255)
                                                                  taataaATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTAAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATAtct                                                                                  >  1:2250833/1‑148 (MQ=255)
                                                                   aataaataTTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATAtctc                                                                                 >  1:2882840/1‑148 (MQ=255)
                                                                    ataaataTTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAATTTCTGTTCAATAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATctct                                                                                <  1:678207/148‑1 (MQ=255)
                                                                      aaataTTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCttt                                                                              >  2:1084291/1‑148 (MQ=255)
                                                                        atatTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGc                                                                            >  1:486872/1‑148 (MQ=255)
                                                                         tatTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCAAGATATCTCTTTGCg                                                                           >  2:2496393/1‑148 (MQ=255)
                                                                          atTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGt                                                                          <  1:1940015/148‑1 (MQ=255)
                                                                          atTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGt                                                                          >  2:2012289/1‑148 (MQ=255)
                                                                          atTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGt                                                                          >  1:248368/1‑148 (MQ=255)
                                                                              ttATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAgg                                                                      >  2:625328/1‑148 (MQ=255)
                                                                                  aaTCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCaa                                                                  >  1:3081602/1‑148 (MQ=255)
                                                                                  aaTCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCaa                                                                  >  2:2079601/1‑148 (MQ=255)
                                                                                  aaTCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCaa                                                                  >  2:1188723/1‑148 (MQ=255)
                                                                                    tCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAata                                                                >  1:2609539/1‑148 (MQ=255)
                                                                                    tCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAata                                                                <  2:1255996/148‑1 (MQ=255)
                                                                                    tCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAata                                                                <  2:2789981/148‑1 (MQ=255)
                                                                                      cTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATAtt                                                              >  1:1045368/1‑148 (MQ=255)
                                                                                             ttGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATtgat                                                       >  2:819757/1‑148 (MQ=255)
                                                                                              tgtgTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATtgatg                                                      <  1:2548037/148‑1 (MQ=255)
                                                                                                 gtTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCNTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATgaga                                                   >  1:1938943/1‑148 (MQ=255)
                                                                                                 gtTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATgaga                                                   <  2:1599921/148‑1 (MQ=255)
                                                                                                 gtTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATgaga                                                   >  1:1143994/1‑148 (MQ=255)
                                                                                                 gtTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATgaga                                                   >  1:2566665/1‑148 (MQ=255)
                                                                                                 gtTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATgaga                                                   >  2:378535/1‑148 (MQ=255)
                                                                                                   tCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGaaa                                                 >  1:1448386/1‑148 (MQ=255)
                                                                                                      ctctTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGAGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAatat                                              >  1:1664669/1‑148 (MQ=255)
                                                                                                      ctctTTAATACGCTCGGCCTTAAAGCTTGCAATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGTCATCCCAAAATAAACCTTCATAAAACGCCAGATATCTCTTAGCGTCAGGCCAATATTCATTGATGAGAAatat                                              >  2:3398077/1‑148 (MQ=255)
                                                                                                       tctTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAAtata                                             <  2:552834/148‑1 (MQ=255)
                                                                                                       tctTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAAtata                                             >  1:2743167/1‑148 (MQ=255)
                                                                                                        ctTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATaa                                            >  2:3096118/1‑148 (MQ=255)
                                                                                                         tttAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATaaa                                           <  2:2838035/148‑1 (MQ=255)
                                                                                                         tttAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATaaa                                           <  2:2224174/148‑1 (MQ=255)
                                                                                                              tACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAAt                                      >  1:2034965/1‑148 (MQ=255)
                                                                                                              tACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAAt                                      <  2:189390/148‑1 (MQ=255)
                                                                                                                cGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCa                                    >  1:892023/1‑148 (MQ=255)
                                                                                                                 gCTCGGCCTTAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAg                                   <  1:3237156/148‑1 (MQ=255)
                                                                                                                     ggCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCTTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCt                               <  1:1588387/148‑1 (MQ=255)
                                                                                                                       ccTTAACGATTGCCATGTTCAATAATTTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCttt                             <  2:3652643/148‑1 (MQ=255)
                                                                                                                         ttAACGCTTGCCATGTTCAATAAATTATGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTAt                           >  2:1565919/1‑148 (MQ=255)
                                                                                                                          tAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATc                          >  1:3750333/1‑148 (MQ=255)
                                                                                                                            aCGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATcgc                        >  2:3576398/1‑148 (MQ=255)
                                                                                                                               cTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGcc                     >  2:3536111/1‑148 (MQ=255)
                                                                                                                                 tGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAg                   >  2:3347631/1‑148 (MQ=255)
                                                                                                                                   ccATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAgcg                 >  1:272127/1‑148 (MQ=255)
                                                                                                                                    cATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAgcgg                <  2:2599679/148‑1 (MQ=255)
                                                                                                                                      tGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAgcggcg              >  2:2121643/1‑148 (MQ=255)
                                                                                                                                      tGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAgcggcg              <  1:1510097/148‑1 (MQ=255)
                                                                                                                                      tGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAgcggcg              <  1:2423299/148‑1 (MQ=255)
                                                                                                                                      tGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAgcggcg              <  2:2906880/148‑1 (MQ=255)
                                                                                                                                        ttCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGcgcg            >  1:714229/1‑148 (MQ=255)
                                                                                                                                         tCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCgg           <  2:3307435/148‑1 (MQ=255)
                                                                                                                                         tCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCgg           >  2:1116473/1‑148 (MQ=255)
                                                                                                                                          cAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGt          >  2:1325373/1‑148 (MQ=255)
                                                                                                                                           aaTAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTa         <  1:598339/148‑1 (MQ=255)
                                                                                                                                                  tCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGc  >  2:3746003/1‑148 (MQ=255)
                                                                                                                                                  tCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGc  >  2:1940254/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCCTTAACGCTTGCCATGTTCAATAAATTCTGTTCAAGAGACAATATTTTGCGTAAAGGCATCCCAAAATAAACCTTCATAAATCGCCAGATATCTCTTTGCGTCAGGCCAATATTCATTGATGAGAAATATAAACCAATCAGGTCTTTATCGCGCCAGCGGCGCGGTACTTTTGC  >  CP009273/3798435‑3798728

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: