Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,342,897 C→A 100% D19Y (GAC→TAC)  inaA ← acid‑inducible Kdo/WaaP family putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,342,8970CA99.5% 883.2 / ‑4.3 269D19Y (GAC→TAC) inaAacid‑inducible Kdo/WaaP family putative kinase
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (115/151);  total (115/154)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 269 total aligned reads displayed.

GGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTCGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCGTATTAATGAATTTCT  >  CP009273/2342750‑2343043
                                                                                                                                                   |                                                                                                                                                  
ggACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGc                                                                                                                                                  <  2:3260374/150‑1 (MQ=255)
 gACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGcc                                                                                                                                                 >  1:1043067/1‑150 (MQ=255)
  aCAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGccc                                                                                                                                                <  2:2148779/150‑1 (MQ=255)
  aCAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGccc                                                                                                                                                <  2:3189754/150‑1 (MQ=255)
   cAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCct                                                                                                                                               >  2:3127789/1‑150 (MQ=255)
   cAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCct                                                                                                                                               >  2:85938/1‑150 (MQ=255)
    aaTCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCctc                                                                                                                                              <  1:1519362/150‑1 (MQ=255)
      tCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTg                                                                                                                                            <  2:2500083/150‑1 (MQ=255)
        gTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGtt                                                                                                                                          <  2:3735672/150‑1 (MQ=255)
         ttGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTGTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTAATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTg                                                                                                                                         <  2:3027520/150‑1 (MQ=255)
            gTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGccc                                                                                                                                      >  2:3076821/1‑150 (MQ=255)
             tCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCcca                                                                                                                                     <  2:377853/150‑1 (MQ=255)
                gCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCccacca                                                                                                                                  <  2:4742999/150‑1 (MQ=255)
                  cGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGt                                                                                                                                >  2:4549279/1‑150 (MQ=255)
                  cGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGt                                                                                                                                >  1:2132041/1‑150 (MQ=255)
                    aaCGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGa                                                                                                                              >  1:4647518/1‑150 (MQ=255)
                    aaCGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGa                                                                                                                              >  1:3502872/1‑150 (MQ=255)
                       ggATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTa                                                                                                                           >  1:4926470/1‑150 (MQ=255)
                         aTAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTaaa                                                                                                                         >  2:2580401/1‑150 (MQ=255)
                             cGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATtcg                                                                                                                     <  2:1311723/150‑1 (MQ=255)
                              gTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATtcgt                                                                                                                    >  2:565823/1‑150 (MQ=255)
                              gTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATtcgt                                                                                                                    <  1:532821/150‑1 (MQ=255)
                                  ggAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTa                                                                                                                >  1:58912/1‑150 (MQ=255)
                                     aTGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACtt                                                                                                             <  1:1185714/150‑1 (MQ=255)
                                      tGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACttt                                                                                                            >  1:1947303/1‑150 (MQ=255)
                                      tGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACttt                                                                                                            <  2:3643845/150‑1 (MQ=255)
                                        aaaCAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGc                                                                                                          >  1:969987/1‑150 (MQ=255)
                                        aaaCAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGc                                                                                                          >  1:1291718/1‑150 (MQ=255)
                                         aaCAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCt                                                                                                         >  1:1076667/1‑150 (MQ=255)
                                           cAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGa                                                                                                       >  2:1721460/1‑150 (MQ=255)
                                            aGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGaa                                                                                                      >  1:2879454/1‑150 (MQ=255)
                                            aGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGaa                                                                                                      <  1:1660343/150‑1 (MQ=255)
                                               tgatgCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACt                                                                                                   <  1:890362/150‑1 (MQ=255)
                                               tgatgCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACt                                                                                                   <  1:305956/150‑1 (MQ=255)
                                                  tgCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGcc                                                                                                <  1:4267690/150‑1 (MQ=255)
                                                  tgCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGcc                                                                                                <  2:1057905/150‑1 (MQ=255)
                                                     gTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATa                                                                                             <  1:1685599/150‑1 (MQ=255)
                                                     gTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATa                                                                                             <  2:4146147/150‑1 (MQ=255)
                                                      tCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATAc                                                                                            >  2:1207810/1‑150 (MQ=255)
                                                       cATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATAcc                                                                                           <  2:2844140/150‑1 (MQ=255)
                                                        aTGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATAccc                                                                                          <  2:2244684/150‑1 (MQ=255)
                                                        aTGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATAccc                                                                                          >  2:2171458/1‑150 (MQ=255)
                                                        aTGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATAccc                                                                                          <  2:4340057/150‑1 (MQ=255)
                                                        aTGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATAccc                                                                                          >  1:3009146/1‑150 (MQ=255)
                                                         tGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCt                                                                                         >  2:4427896/1‑150 (MQ=255)
                                                         tGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCt                                                                                         >  1:488543/1‑150 (MQ=255)
                                                          gcgcTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTa                                                                                        >  1:1849808/1‑150 (MQ=255)
                                                           cgcTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTAt                                                                                       <  2:4863733/150‑1 (MQ=255)
                                                            gcTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATc                                                                                      <  1:2932713/150‑1 (MQ=255)
                                                                 acatacaGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTg                                                                                 <  1:2011350/150‑1 (MQ=255)
                                                                 acatacaGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTg                                                                                 <  2:2360333/150‑1 (MQ=255)
                                                                  catacaGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGt                                                                                <  2:2019180/150‑1 (MQ=255)
                                                                   atacaGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTc                                                                               <  2:2129079/150‑1 (MQ=255)
                                                                       aGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCAttt                                                                           >  2:2960668/1‑150 (MQ=255)
                                                                        gCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCAtttt                                                                          <  2:2751724/150‑1 (MQ=255)
                                                                         cTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCAttttt                                                                         <  2:4826405/150‑1 (MQ=255)
                                                                           ttttGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTttat                                                                       >  2:2827731/1‑150 (MQ=255)
                                                                             ttGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTttatta                                                                     <  1:544906/150‑1 (MQ=255)
                                                                              tGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTaa                                                                    >  1:49435/1‑150 (MQ=255)
                                                                                ccGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACg                                                                  >  2:3510677/1‑150 (MQ=255)
                                                                                 cGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGa                                                                 >  1:2970365/1‑150 (MQ=255)
                                                                                  gTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGAc                                                                <  2:140136/150‑1 (MQ=255)
                                                                                     gcgcTCGACTCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTTATTAAATTCGTCGTACTTTGCTGAAACTGCAATACCCTATCGCCTGTCATTTTTATTAACGACAAt                                                             <  2:3938870/150‑1 (MQ=255)
                                                                                        cTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGAc                                                          >  1:1451078/1‑150 (MQ=255)
                                                                                         ccGACGCACTGCACCCCGCTCCTACCATCACGGCGATAGTTAGGCTCTTCAACCCAGTCGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACt                                                         <  1:2685124/149‑1 (MQ=255)
                                                                                          cGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACta                                                        >  2:3298564/1‑150 (MQ=255)
                                                                                           gACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACtat                                                       >  1:1031792/1‑150 (MQ=255)
                                                                                              gCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAgg                                                    <  1:1268455/150‑1 (MQ=255)
                                                                                                aTTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAggtg                                                  >  2:1109373/1‑150 (MQ=255)
                                                                                                  tGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAggtggt                                                >  2:3341070/1‑150 (MQ=255)
                                                                                                   gCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGtt                                               <  2:3831888/150‑1 (MQ=255)
                                                                                                   gCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGtt                                               >  1:336689/1‑150 (MQ=255)
                                                                                                      ccccGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTAcc                                            >  2:1012894/1‑150 (MQ=255)
                                                                                                       cccGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCt                                           >  1:720422/1‑150 (MQ=255)
                                                                                                          gCTCATACCTTTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTTTTGCCCACCAGTGATTAATTTCGTCGTACTGTGCTGAAACTGCCATACCCTATCGCCTTTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAg                                        <  1:3127789/150‑1 (MQ=255)
                                                                                                           cTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAgg                                       <  2:2960830/150‑1 (MQ=255)
                                                                                                           cTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAgg                                       <  1:5110982/150‑1 (MQ=255)
                                                                                                              aTACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGaaa                                    >  2:1979304/1‑150 (MQ=255)
                                                                                                                 ccATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATc                                 <  2:944745/150‑1 (MQ=255)
                                                                                                                  cATTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCt                                <  1:841984/150‑1 (MQ=255)
                                                                                                                   aTTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCtt                               <  2:5014071/150‑1 (MQ=255)
                                                                                                                   aTTACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCtt                               >  1:2072533/1‑150 (MQ=255)
                                                                                                                    ttACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTa                              <  2:4926416/150‑1 (MQ=255)
                                                                                                                     tACGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTaa                             >  2:4705588/1‑150 (MQ=255)
                                                                                                                       cGGCGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATg                           >  2:2900783/1‑150 (MQ=255)
                                                                                                                        ggCGATAGTTAGGCACTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGa                          <  2:3487324/150‑1 (MQ=255)
                                                                                                                          cGATAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGaaa                        >  1:1418921/1‑150 (MQ=255)
                                                                                                                            aTAGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACg                      >  2:2651766/1‑150 (MQ=255)
                                                                                                                              aGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACgtg                    <  2:1576553/150‑1 (MQ=255)
                                                                                                                              aGTTAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCGGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACgtg                    <  2:1341794/150‑1 (MQ=255)
                                                                                                                                ttAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTc                  >  2:1191758/1‑150 (MQ=255)
                                                                                                                                ttAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTc                  >  1:295402/1‑150 (MQ=255)
                                                                                                                                 tAGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCg                 <  2:759329/150‑1 (MQ=255)
                                                                                                                                  aGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCGt                >  1:3945977/1‑150 (MQ=255)
                                                                                                                                  aGGCTCTTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCGt                >  2:2314231/1‑150 (MQ=255)
                                                                                                                                     ctctTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCGTAtt             <  1:4506897/150‑1 (MQ=255)
                                                                                                                                     ctctTCAACCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCGTAtt             <  1:3575643/150‑1 (MQ=255)
                                                                                                                                           aaCCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCGTATTAATGaa       >  1:4541986/1‑150 (MQ=255)
                                                                                                                                           aaCCCAGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCGTATTAATGaa       >  1:2503609/1‑150 (MQ=255)
                                                                                                                                                aGTAGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCGTATTAATGAATTTCt  >  1:402668/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGACAATCGTTGGTCGGCCGAACGGATAACGTACGGAATGAAACAGATGATGCGTCATGCGCTTTACATACAGCTTTTTGCCGTTGCGCTCGACGCATTGCACCCCGCTCATACCATTACGGCGATAGTTAGGCTCTTCAACCCAGTCGCCCTCTGTTGCCCACCAGTGATTAAATTCGTCGTACTTTGCTGAAACTGCCATACCCTATCGCCTGTCATTTTTATTAACGACAATGACTATAGGTGGTTACCTGAGGAAAATCTTAATGAAACGTGTCGTATTAATGAATTTCT  >  CP009273/2342750‑2343043

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: