Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,799,108 C→A 47.5% R68L (CGG→CTG)  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,799,1080CA47.5% 171.6 / inf 275R68L (CGG→CTG) waaPkinase that phosphorylates core heptose of lipopolysaccharide
Reads supporting (aligned to +/- strand):  ref base C (58/86);  new base A (67/64);  total (125/150)
Fisher's exact test for biased strand distribution p-value = 8.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

TGAGATCTTCGGTAATAATAAATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTG  >  CP009273/3798962‑3799255
                                                                                                                                                  |                                                                                                                                                   
tGAGATCTTCGGTAATAATAAATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCa                                                                                                                                                  <  1:3825454/150‑1 (MQ=255)
    aTCTTCGGTAATAATAAATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATga                                                                                                                                              >  1:2838911/1‑150 (MQ=255)
     tCTTCGGTAATAATAAATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATgag                                                                                                                                             <  1:776359/150‑1 (MQ=255)
      cTTCGGTAATAATAAATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATgaga                                                                                                                                            <  1:1794057/150‑1 (MQ=255)
         cGGTAATAATAAATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTa                                                                                                                                         <  1:3942087/150‑1 (MQ=255)
            taataataaATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGAt                                                                                                                                      >  1:547022/1‑150 (MQ=255)
             aataataaATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGAtt                                                                                                                                     >  1:3107693/1‑150 (MQ=255)
              ataataaATGATGCGCGAGTTAATGGATTTAACCCTTGTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGAttt                                                                                                                                    <  2:489565/150‑1 (MQ=255)
                aataaATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGAttttt                                                                                                                                  <  1:2751992/150‑1 (MQ=255)
                    aaTGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTtata                                                                                                                              >  1:1150064/1‑150 (MQ=255)
                     atgatgCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATaa                                                                                                                             <  2:2488455/150‑1 (MQ=255)
                     atgatgCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATaa                                                                                                                             <  1:2082865/150‑1 (MQ=255)
                           cgcgAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCtt                                                                                                                       <  2:1444009/150‑1 (MQ=255)
                           cgcgAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCtt                                                                                                                       <  2:2876892/150‑1 (MQ=255)
                             cgAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCtttt                                                                                                                     <  2:792549/150‑1 (MQ=255)
                             cgAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCAAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCTTCCGCAATGAGAGTAGATTTTTTATAATCTCtttt                                                                                                                     <  1:3843438/150‑1 (MQ=255)
                               aGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTGCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTaa                                                                                                                   >  2:4563591/1‑150 (MQ=255)
                                   aaTGGATTTAACCCTTTTTCGCCAAACCCACTGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATgtc                                                                                                               <  1:3885755/150‑1 (MQ=255)
                                   aaTGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATgtc                                                                                                               <  1:2986962/150‑1 (MQ=255)
                                    aTGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATgtcg                                                                                                              >  2:1490960/1‑150 (MQ=255)
                                       gATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTcc                                                                                                           >  1:892369/1‑150 (MQ=255)
                                       gATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTcc                                                                                                           >  1:3271442/1‑150 (MQ=255)
                                          ttAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCt                                                                                                        >  1:4053168/1‑150 (MQ=255)
                                           tAACCCTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTtg                                                                                                      >  2:2307576/1‑150 (MQ=255)
                                            aaCCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTtg                                                                                                      <  2:1214596/150‑1 (MQ=255)
                                             aCCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTtgt                                                                                                     <  1:1388661/150‑1 (MQ=255)
                                               ccTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGc                                                                                                   >  2:2342412/1‑150 (MQ=255)
                                                cTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGcc                                                                                                  <  1:3557498/150‑1 (MQ=255)
                                                 tttttCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCa                                                                                                 <  1:55307/150‑1 (MQ=255)
                                                 tttttCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCa                                                                                                 <  1:113069/150‑1 (MQ=255)
                                                 tttttCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCa                                                                                                 <  2:2209433/150‑1 (MQ=255)
                                                 tttttCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCCGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCa                                                                                                 >  2:1272187/1‑150 (MQ=255)
                                                    ttCGCCAAACCCAAGGCCCTTCATTGTATCATCGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCAttt                                                                                              <  1:1230647/150‑1 (MQ=255)
                                                        ccAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGa                                                                                          <  1:1838083/150‑1 (MQ=255)
                                                        ccAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGa                                                                                          >  1:2819035/1‑150 (MQ=255)
                                                         cAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGaa                                                                                         >  2:4035459/1‑150 (MQ=255)
                                                          aaaCCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGaaa                                                                                        >  1:3707681/1‑150 (MQ=255)
                                                          aaaCCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGaaa                                                                                        <  2:4575353/150‑1 (MQ=255)
                                                           aaCCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGaaaa                                                                                       >  2:1234651/1‑150 (MQ=255)
                                                           aaCCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGaaaa                                                                                       <  1:432880/150‑1 (MQ=255)
                                                            aCCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAAt                                                                                      <  2:3931826/150‑1 (MQ=255)
                                                              ccAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAg                                                                                    <  2:4204291/150‑1 (MQ=255)
                                                              ccAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAg                                                                                    <  1:142074/150‑1 (MQ=255)
                                                               cAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGc                                                                                   >  2:3031826/1‑150 (MQ=255)
                                                               cAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGc                                                                                   <  2:1064813/150‑1 (MQ=255)
                                                                 aTGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCtt                                                                                 >  1:4471547/1‑150 (MQ=255)
                                                                  tGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCttt                                                                                <  2:4225014/150‑1 (MQ=255)
                                                                  tGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAgaatt                                                                                >  2:261875/1‑146 (MQ=255)
                                                                     ccTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTcc                                                                             >  1:4038682/1‑150 (MQ=255)
                                                                      cTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTccc                                                                            <  2:2965628/150‑1 (MQ=255)
                                                                      cTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTccc                                                                            >  1:1526890/1‑150 (MQ=255)
                                                                        tCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCgg                                                                          <  2:2923329/150‑1 (MQ=255)
                                                                          aTTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGAc                                                                        >  1:1897880/1‑150 (MQ=255)
                                                                          aTTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGAc                                                                        <  2:3976816/150‑1 (MQ=255)
                                                                          aTTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGAc                                                                        <  2:930834/150‑1 (MQ=255)
                                                                            tGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAg                                                                      >  1:3968147/1‑150 (MQ=255)
                                                                            tGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAg                                                                      >  1:2628485/1‑150 (MQ=255)
                                                                             gTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGt                                                                     <  2:1907892/150‑1 (MQ=255)
                                                                               aTCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTc                                                                   >  2:645487/1‑150 (MQ=255)
                                                                                   aCGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAg                                                               >  2:1027315/1‑150 (MQ=255)
                                                                                   aCGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAg                                                               >  2:3359969/1‑150 (MQ=255)
                                                                                         aCATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAAt                                                         <  1:1703568/150‑1 (MQ=255)
                                                                                         aCATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAAAAGCTTTTCCCGGACAGTTCAAAGCGTAAt                                                         >  2:3439157/1‑150 (MQ=255)
                                                                                          cATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATg                                                        <  2:226245/150‑1 (MQ=255)
                                                                                             cACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTac                                                     <  2:3521332/150‑1 (MQ=255)
                                                                                                tCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTacgac                                                  <  1:99668/150‑1 (MQ=255)
                                                                                                 cAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTacgacg                                                 <  2:3791861/150‑1 (MQ=255)
                                                                                                        tGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTcc                                          >  1:177831/1‑150 (MQ=255)
                                                                                                         gAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCa                                         <  1:4198599/150‑1 (MQ=255)
                                                                                                          aaTAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAg                                        >  1:3328422/1‑150 (MQ=255)
                                                                                                          aaTAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAg                                        <  2:1397582/150‑1 (MQ=255)
                                                                                                          aaTAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAg                                        <  2:3290002/150‑1 (MQ=255)
                                                                                                           aTAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGt                                       >  2:1568918/1‑150 (MQ=255)
                                                                                                             aGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTc                                     >  1:1843276/1‑150 (MQ=255)
                                                                                                               cGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCAc                                   >  1:1656258/1‑150 (MQ=255)
                                                                                                                gTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACg                                  <  2:970435/150‑1 (MQ=255)
                                                                                                                gTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACg                                  >  2:890164/1‑150 (MQ=255)
                                                                                                                gTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACg                                  <  1:2804713/150‑1 (MQ=255)
                                                                                                                gTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACg                                  <  1:3538345/150‑1 (MQ=255)
                                                                                                                  gCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGaa                                >  2:2975354/1‑150 (MQ=255)
                                                                                                                        ctctGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCt                          >  2:1542430/1‑150 (MQ=255)
                                                                                                                          ctgtctgCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTGTCCCGGACAGTTCTAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCg                        <  2:1776189/150‑1 (MQ=255)
                                                                                                                               tgCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTACTGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGtt                   <  2:2658736/150‑1 (MQ=255)
                                                                                                                               tgCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGtt                   >  1:2535997/1‑150 (MQ=255)
                                                                                                                                gcgcCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTc                  <  1:2255660/150‑1 (MQ=255)
                                                                                                                                 cgcCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCa                 <  2:4443301/150‑1 (MQ=255)
                                                                                                                                  gcCCAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAg                >  1:1338726/1‑150 (MQ=255)
                                                                                                                                   cccAAAACGGGCATCAGCTATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCACAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGt               <  2:1579361/150‑1 (MQ=255)
                                                                                                                                   cccAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGt               >  2:3939416/1‑150 (MQ=255)
                                                                                                                                    ccAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAATACCTCGCCGTTCAGtt              >  1:4356002/1‑150 (MQ=255)
                                                                                                                                     cAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGttt             <  2:1861293/150‑1 (MQ=255)
                                                                                                                                     cAAAACGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGttt             <  1:1109776/150‑1 (MQ=255)
                                                                                                                                      aaaaCGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGtttt            <  2:1642582/150‑1 (MQ=255)
                                                                                                                                      aaaaCGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGtttt            <  2:2140584/150‑1 (MQ=255)
                                                                                                                                      aaaaCGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGtttt            <  1:3621520/150‑1 (MQ=255)
                                                                                                                                       aaaCGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGttttt           <  1:837119/150‑1 (MQ=255)
                                                                                                                                       aaaCGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGttttt           >  2:2461912/1‑150 (MQ=255)
                                                                                                                                        aaCGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGtttttt          <  2:4442055/150‑1 (MQ=255)
                                                                                                                                          cGGGCATCAGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGa        <  2:1872589/150‑1 (MQ=255)
                                                                                                                                           gggCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAATACCTCGCCGTTCAGTTTTTTGAc       >  1:2208969/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGAGATCTTCGGTAATAATAAATGATGCGCGAGTTAATGGATTTAACCCTTTTTCGCCAAACCCAATGCCCTTCATTGTATCAACGCCAACATCACTCAGGCGATGAATAGCGTGCCACTCTCTGTCTGCGCCCAAAACGGGCATCCGCAATGAGAGTAGATTTTTTATAATCTCTTTTAATGTCGTCCCCTTGTGCCATTTAAGAAAATAGCTTTTCCCGGACAGTTCAAAGCGTAATGTACGACGAGTCTCCAGTTCACGAAAGACCTCGCCGTTCAGTTTTTTGACCTCTG  >  CP009273/3798962‑3799255

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: