Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,584,756 A→T 100% *305K (TAA→AAA)  yjjM ← yjjMN operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,584,7560AT99.1% 1065.1 / ‑1.4 321*305K (TAA→AAA) yjjMyjjMN operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (156/162);  total (158/163)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 321 total aligned reads displayed.

CAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGA  >  CP009273/4584609‑4584902
                                                                                                                                                   |                                                                                                                                                  
cAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTAt                                                                                                                                                  >  2:3901961/1‑150 (MQ=255)
cAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTAATCTTTTGATAAAACCGTGCACCCGCGACTAGTGTGCTGTTTTTAATCCGTACTTTAt                                                                                                                                                  >  1:516442/1‑150 (MQ=255)
  ggTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGa                                                                                                                                                <  2:2809252/150‑1 (MQ=255)
   gTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAg                                                                                                                                               <  2:3521613/150‑1 (MQ=255)
    ttAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGc                                                                                                                                              <  1:4463445/150‑1 (MQ=255)
       aaGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATa                                                                                                                                           <  2:5245599/150‑1 (MQ=255)
        aGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATaa                                                                                                                                          >  2:2666562/1‑150 (MQ=255)
         gTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAAc                                                                                                                                         <  1:3321827/150‑1 (MQ=255)
          ttttAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAAcg                                                                                                                                        <  1:887230/150‑1 (MQ=255)
              aGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCgtg                                                                                                                                    >  2:2885472/1‑150 (MQ=255)
              aGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCgtg                                                                                                                                    >  1:3948795/1‑150 (MQ=255)
               ggACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCgtgt                                                                                                                                   >  1:1395696/1‑150 (MQ=255)
                 aCTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGtta                                                                                                                                 >  2:5167894/1‑149 (MQ=255)
                  cTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGttat                                                                                                                                <  1:1809377/150‑3 (MQ=255)
                       ctacCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTa                                                                                                                           >  2:3886024/1‑150 (MQ=255)
                       ctacCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTa                                                                                                                           <  1:3978308/150‑1 (MQ=255)
                       ctacCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTa                                                                                                                           >  1:3390746/1‑150 (MQ=255)
                          ccGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATg                                                                                                                        >  2:2266429/1‑150 (MQ=255)
                           cGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATgg                                                                                                                       >  2:917352/1‑150 (MQ=255)
                             cTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGAt                                                                                                                     <  2:3992650/150‑1 (MQ=255)
                             cTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGAt                                                                                                                     <  1:3817062/150‑1 (MQ=255)
                              tCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTATTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATc                                                                                                                    <  1:3495507/150‑1 (MQ=255)
                                acaAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGa                                                                                                                  >  1:4991200/1‑150 (MQ=255)
                                 caAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGaa                                                                                                                 >  1:4692819/1‑150 (MQ=255)
                                    aaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCa                                                                                                              >  1:2712586/1‑150 (MQ=255)
                                       cATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCattg                                                                                                           <  2:36692/150‑1 (MQ=255)
                                        aTTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCattga                                                                                                          <  2:1853075/150‑1 (MQ=255)
                                          ttttAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCattgatt                                                                                                        <  2:2116321/150‑1 (MQ=255)
                                           tttAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCattgattg                                                                                                       <  1:459421/150‑1 (MQ=255)
                                             tAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGtt                                                                                                     >  1:727237/1‑150 (MQ=255)
                                              aGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGttt                                                                                                    <  2:4682296/150‑1 (MQ=255)
                                                    cTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGa                                                                                              >  2:4990066/1‑150 (MQ=255)
                                                    cTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGa                                                                                              >  1:546049/1‑150 (MQ=255)
                                                      gATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAAt                                                                                            >  1:663312/1‑150 (MQ=255)
                                                       aTTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAAtt                                                                                           <  1:3169434/150‑1 (MQ=255)
                                                            gAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAAt                                                                                      >  1:5322471/1‑150 (MQ=255)
                                                             aGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATg                                                                                     >  2:2387460/1‑150 (MQ=255)
                                                             aGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATg                                                                                     >  1:3042524/1‑150 (MQ=255)
                                                             aGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATg                                                                                     <  1:1744215/150‑1 (MQ=255)
                                                               cATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGt                                                                                   >  1:346643/1‑150 (MQ=255)
                                                               cATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGt                                                                                   >  1:1597321/1‑150 (MQ=255)
                                                               cATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGt                                                                                   <  2:4812544/150‑1 (MQ=255)
                                                                aTAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGtt                                                                                  >  2:3105371/1‑150 (MQ=255)
                                                                 tAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTa                                                                                 <  1:3718272/150‑1 (MQ=255)
                                                                  tgttAAGCGGGAATTTCCCGGTATTTATCTTTTGATATAACCGTGCACCCGCTAATAGTTTGCTGTTTTTTATCCTTACTTTATGATCATAAAGCGTGTTATCATTTATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATTGTTAc                                                                                <  2:1917181/147‑1 (MQ=255)
                                                                   ggTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACg                                                                               <  1:3989821/150‑1 (MQ=255)
                                                                   ggTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACg                                                                               >  2:4586565/1‑150 (MQ=255)
                                                                   ggTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACg                                                                               <  1:3573820/150‑1 (MQ=255)
                                                                      aaGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAg                                                                            <  1:2214352/150‑1 (MQ=255)
                                                                      aaGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAg                                                                            <  1:2972258/150‑1 (MQ=255)
                                                                      aaGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAg                                                                            >  1:2356623/1‑150 (MQ=255)
                                                                      aaGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAg                                                                            <  1:113811/150‑1 (MQ=255)
                                                                        gCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAgcg                                                                          <  1:1258044/150‑1 (MQ=255)
                                                                         cGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAgcgc                                                                         >  1:1886490/1‑150 (MQ=255)
                                                                             aaTTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAgg                                                                     >  2:619524/1‑150 (MQ=255)
                                                                              aTTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAggg                                                                    <  1:2437138/150‑1 (MQ=255)
                                                                              aTTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAggg                                                                    >  1:2599431/1‑150 (MQ=255)
                                                                                   ccGGTATTTATCTTTTGATAAAACCGTGCACCCGCTCATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTCGCTGACTTCAAATGGGTACGGAGCGCACGGGTGGCa                                                               >  1:3912315/1‑150 (MQ=255)
                                                                                     ggTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATc                                                             >  1:492297/1‑150 (MQ=255)
                                                                                      gTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCa                                                            <  2:116449/150‑1 (MQ=255)
                                                                                      gTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCa                                                            >  2:1941745/1‑150 (MQ=255)
                                                                                       tATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCaa                                                           >  2:3814845/1‑150 (MQ=255)
                                                                                        aTTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAg                                                          <  2:2801446/150‑1 (MQ=255)
                                                                                            aTCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCa                                                      >  1:2020690/1‑150 (MQ=255)
                                                                                              ctattGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAATGGTGGCATCAATGTCACt                                                    <  1:3136772/147‑1 (MQ=255)
                                                                                                tttGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGc                                                  <  2:3605316/150‑1 (MQ=255)
                                                                                                tttGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGc                                                  >  2:2531321/1‑150 (MQ=255)
                                                                                                  tGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATTGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACTGAGCGCAAGGGTGGCATCAAGGTCACTGCgg                                                >  2:3233321/1‑150 (MQ=255)
                                                                                                       aaaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGATCGCAAGGGTGGCATCAAGGTCACTGCGGCAAAt                                           <  2:1442611/150‑1 (MQ=255)
                                                                                                       aaaCCGTGCACCCGCTAAGAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGATCGCAAGGGTGGCATCAAGGTCACTGCGGCAAAt                                           <  2:1442758/150‑1 (MQ=255)
                                                                                                        aaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATc                                          <  1:3245629/150‑1 (MQ=255)
                                                                                                         aCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCa                                         <  2:4123260/150‑1 (MQ=255)
                                                                                                         aCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCa                                         <  2:433199/150‑1 (MQ=255)
                                                                                                           cGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCaga                                       >  1:5203504/1‑150 (MQ=255)
                                                                                                           cGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCaga                                       >  1:2786840/1‑150 (MQ=255)
                                                                                                            gTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGTGGCAAATCagag                                      >  1:4187421/1‑150 (MQ=255)
                                                                                                             tGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAgc                                     <  1:104698/150‑1 (MQ=255)
                                                                                                              gCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAgcn                                    <  1:4886828/150‑2 (MQ=255)
                                                                                                                aCCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCt                                  >  1:4749917/1‑150 (MQ=255)
                                                                                                                      tAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAAt                            >  2:693490/1‑150 (MQ=255)
                                                                                                                          aGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGACTGGTc                        >  2:5076757/1‑150 (MQ=255)
                                                                                                                           gTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCa                       <  2:4846717/150‑1 (MQ=255)
                                                                                                                              tGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCatat                    >  2:2201679/1‑150 (MQ=255)
                                                                                                                              tGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCatat                    <  2:1678548/150‑1 (MQ=255)
                                                                                                                                 tGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCt                 >  2:3095338/1‑150 (MQ=255)
                                                                                                                                 tGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCt                 <  1:3686218/150‑1 (MQ=255)
                                                                                                                                  gTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTc                >  1:2844944/1‑150 (MQ=255)
                                                                                                                                   tttttAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCa               <  1:1969533/150‑1 (MQ=255)
                                                                                                                                   tttttAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCa               <  2:3030879/150‑1 (MQ=255)
                                                                                                                                     tttAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATc             <  2:1258985/150‑1 (MQ=255)
                                                                                                                                      ttAATCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCg            <  1:964646/150‑1 (MQ=255)
                                                                                                                                        aaTCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGAc          <  2:5012383/150‑1 (MQ=255)
                                                                                                                                        aaTCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGAc          <  2:839712/150‑1 (MQ=255)
                                                                                                                                          tCCGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACgg        >  1:4269647/1‑150 (MQ=255)
                                                                                                                                            cGTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCa      >  2:1343888/1‑150 (MQ=255)
                                                                                                                                             gTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCaa     >  1:4846669/1‑150 (MQ=255)
                                                                                                                                             gTACTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCAGAATTCAAAAGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGACTGGTCATATGCTCATCGACGGCaa     >  2:5316569/1‑150 (MQ=255)
                                                                                                                                               aCTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATg   >  2:234901/1‑150 (MQ=255)
                                                                                                                                                cTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGa  >  2:1808121/1‑150 (MQ=255)
                                                                                                                                                cTTTATGAGCATAACGCGTGTTATCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGa  <  2:3727941/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGA  >  CP009273/4584609‑4584902

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: