Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,407,969 C→A 100% N249K (AAC→AAA acrE → cytoplasmic membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,407,9690CA99.0% 999.9 / 2.5 313coding (747/1158 nt)acrEcytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base A (154/155);  minor base . (1/2);  total (155/158)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 319 total aligned reads displayed.

TGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATTGATGA  >  CP009273/3407822‑3408115
                                                                                                                                                   |                                                                                                                                                  
tGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAgg                                                                                                                                                  <  2:393962/150‑1 (MQ=255)
      ggCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAAc                                                                                                                                            >  1:4645068/1‑150 (MQ=255)
      ggCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAAc                                                                                                                                            <  1:1506573/150‑1 (MQ=255)
      ggCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAAc                                                                                                                                            >  2:4135506/1‑150 (MQ=255)
        cGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCt                                                                                                                                          <  1:2978397/150‑1 (MQ=255)
         gACGGTCCAGCAGCTCGAGCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATACGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTa                                                                                                                                         <  2:610097/150‑1 (MQ=255)
            tGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATcc                                                                                                                                      >  1:1708161/1‑150 (MQ=255)
            tGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATcc                                                                                                                                      >  1:4365278/1‑150 (MQ=255)
               ccAGCAGCTCTATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAATTTTGCATAAGTAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCt                                                                                                                                   <  1:4940215/150‑1 (MQ=255)
               ccAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCt                                                                                                                                   <  1:4270419/150‑1 (MQ=255)
                ctgctgcTCGATCATATCTACGTTTATGTGACCCAATCAATAAACTATTTTATGTGGCTGAAGAAAACCGTAGATAAAGGAAATTTGAATAAGTAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTg                                                                                                                                  <  1:1038988/145‑1 (MQ=255)
                cagcagCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTg                                                                                                                                  <  2:2123901/150‑1 (MQ=255)
                 agcagcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGa                                                                                                                                 <  2:1039448/150‑1 (MQ=255)
                 agcagcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGa                                                                                                                                 >  1:3758886/1‑150 (MQ=255)
                 agcagcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGa                                                                                                                                 >  2:2821650/1‑150 (MQ=255)
                 agcagcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGa                                                                                                                                 >  2:3066044/1‑150 (MQ=255)
                    agcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAg                                                                                                                              >  1:589935/1‑150 (MQ=255)
                    agcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAg                                                                                                                              >  1:658864/1‑150 (MQ=255)
                      cTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAATCCTATCCCCTGAAAGGt                                                                                                                            >  2:4838295/1‑150 (MQ=255)
                      cTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGt                                                                                                                            <  2:5227830/150‑1 (MQ=255)
                        cGATCCTATCTACTTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTAc                                                                                                                          <  1:4052230/150‑1 (MQ=255)
                        cGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTAc                                                                                                                          <  1:5245397/150‑1 (MQ=255)
                        cGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTAc                                                                                                                          <  1:351470/150‑1 (MQ=255)
                         gATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACg                                                                                                                         >  2:2553615/1‑150 (MQ=255)
                          aTCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGc                                                                                                                        >  2:3821081/1‑150 (MQ=255)
                          aTCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGc                                                                                                                        >  2:3578644/1‑150 (MQ=255)
                           tCCTATCTCCGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGTGCAAGGAAATTGGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCt                                                                                                                       <  2:4659784/150‑1 (MQ=255)
                             cTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGc                                                                                                                     >  1:1965518/1‑150 (MQ=255)
                                   aCGTTGATGTGACCCAATCCCGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtct                                                                                                               <  2:732658/150‑1 (MQ=255)
                                   aCGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtct                                                                                                               <  1:4805666/150‑1 (MQ=255)
                                     gTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGTAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTcc                                                                                                             <  1:4554751/150‑1 (MQ=255)
                                      ttGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCg                                                                                                            <  1:2234214/150‑1 (MQ=255)
                                       tGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGa                                                                                                           <  2:4294088/150‑1 (MQ=255)
                                        gATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAt                                                                                                          >  2:1961573/1‑150 (MQ=255)
                                         aTGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtg                                                                                                         <  1:5084207/150‑1 (MQ=255)
                                              aCCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGAcc                                                                                                    <  2:4501633/150‑1 (MQ=255)
                                               cccAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCg                                                                                                   >  1:3966643/1‑150 (MQ=255)
                                                    tCCAGCAACTATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGTAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtgat                                                                                              <  1:3586659/150‑1 (MQ=255)
                                                     ccAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtgatg                                                                                             >  2:2321668/1‑150 (MQ=255)
                                                       aGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGaa                                                                                           >  2:2157737/1‑150 (MQ=255)
                                                        gCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGTAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGaaa                                                                                          <  2:3634331/150‑1 (MQ=255)
                                                        gCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGaaa                                                                                          >  2:2942755/1‑150 (MQ=255)
                                                           aCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCa                                                                                       <  1:2876950/150‑1 (MQ=255)
                                                                  tATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTc                                                                                <  1:757894/150‑1 (MQ=255)
                                                                   aTGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTcc                                                                               >  2:4163294/1‑150 (MQ=255)
                                                                   aTGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTcc                                                                               >  2:1380162/1‑150 (MQ=255)
                                                                      aGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATa                                                                            <  1:5183451/150‑1 (MQ=255)
                                                                          tGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAAccc                                                                        >  2:1719986/1‑150 (MQ=255)
                                                                          tGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAAccc                                                                        >  2:1252656/1‑150 (MQ=255)
                                                                          tGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAAccc                                                                        <  1:1104920/150‑1 (MQ=255)
                                                                             aGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTAc                                                                     <  1:2899479/150‑1 (MQ=255)
                                                                                aaTCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTg                                                                  >  2:2622884/1‑150 (MQ=255)
                                                                                   ccGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTg                                                               <  2:2475106/150‑1 (MQ=255)
                                                                                    cGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGt                                                              >  1:1203976/1‑150 (MQ=255)
                                                                                    cGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGt                                                              >  2:2012820/1‑150 (MQ=255)
                                                                                      tAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCt                                                            >  1:1797165/1‑150 (MQ=255)
                                                                                       agagCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCtt                                                           <  2:855122/150‑1 (MQ=255)
                                                                                       agagCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCtt                                                           >  1:2479499/1‑150 (MQ=255)
                                                                                          gCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTccc                                                        <  2:4290320/150‑1 (MQ=255)
                                                                                          gCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTccc                                                        <  2:3713478/150‑1 (MQ=255)
                                                                                           cAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCt                                                       >  2:2657027/1‑150 (MQ=255)
                                                                                            aaGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTa                                                      >  2:3357783/1‑150 (MQ=255)
                                                                                            aaGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTACGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGATGTCTTCCCTa                                                      >  1:838325/1‑150 (MQ=255)
                                                                                             aGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAAc                                                    <  1:2410214/150‑1 (MQ=255)
                                                                                               gAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAAcc                                                   >  1:3478056/1‑150 (MQ=255)
                                                                                                   tttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCa                                               <  1:3059720/150‑1 (MQ=255)
                                                                                                       cATCAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAAGTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACAt                                           <  2:1593621/150‑1 (MQ=255)
                                                                                                       cATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACAt                                           <  2:4254818/150‑1 (MQ=255)
                                                                                                        aTAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                          >  1:2780482/1‑150 (MQ=255)
                                                                                                        aTAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                          >  2:4596243/1‑150 (MQ=255)
                                                                                                        aTAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                          >  1:4230670/1‑150 (MQ=255)
                                                                                                         tAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATAc                                         >  1:3272485/1‑150 (MQ=255)
                                                                                                          aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACg                                        >  2:3308755/1‑150 (MQ=255)
                                                                                                              aaaaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCttt                                    <  1:3107713/150‑1 (MQ=255)
                                                                                                               aaaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCtttt                                   <  2:3117301/150‑1 (MQ=255)
                                                                                                               aaaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCtttt                                   >  2:1007881/1‑150 (MQ=255)
                                                                                                                aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTg                                  <  1:1375265/150‑1 (MQ=255)
                                                                                                                aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTg                                  <  1:2551952/150‑1 (MQ=255)
                                                                                                                 aCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGc                                 <  1:1421891/150‑1 (MQ=255)
                                                                                                                    ccaccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCgg                              >  2:2258468/1‑150 (MQ=255)
                                                                                                                    ccaccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCgg                              <  2:1867742/150‑1 (MQ=255)
                                                                                                                     caccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCggg                             >  2:528550/1‑150 (MQ=255)
                                                                                                                        caGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTAt                          <  1:2086615/150‑1 (MQ=255)
                                                                                                                            aaCGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGttt                      >  1:4006330/1‑150 (MQ=255)
                                                                                                                            aaCGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGttt                      >  2:1402257/1‑150 (MQ=255)
                                                                                                                               gTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTtgtg                   >  2:4004287/1‑150 (MQ=255)
                                                                                                                                  gagTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTgtgcgt                <  2:3444539/150‑1 (MQ=255)
                                                                                                                                   agTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTgtgcgtg               <  1:186150/150‑1 (MQ=255)
                                                                                                                                     ttGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCa             >  1:5124149/1‑150 (MQ=255)
                                                                                                                                      tGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCAc            <  1:3821081/150‑1 (MQ=255)
                                                                                                                                      tGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCAc            <  2:3534282/150‑1 (MQ=255)
                                                                                                                                        gTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACgg          >  1:4743886/1‑150 (MQ=255)
                                                                                                                                          cATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGAt        <  2:4349380/150‑1 (MQ=255)
                                                                                                                                           aTGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGAtt       <  2:2087230/150‑1 (MQ=255)
                                                                                                                                            tGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtg      <  2:1657655/150‑1 (MQ=255)
                                                                                                                                            tGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtg      >  1:225892/1‑150 (MQ=255)
                                                                                                                                            tGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtg      <  2:4867438/150‑1 (MQ=255)
                                                                                                                                              gAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtgat    <  2:2050538/150‑1 (MQ=255)
                                                                                                                                               aaaaAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtgatg   >  1:3367954/1‑150 (MQ=255)
                                                                                                                                                aaaaGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtgatga  >  2:2891069/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGAACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATTGATGA  >  CP009273/3407822‑3408115

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: