Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,003,597 G→T 100% D100Y (GAC→TAC)  rlmL → fused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,003,5970GT99.2% 849.9 / ‑2.8 258D100Y (GAC→TAC) rlmLfused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (131/125);  total (133/125)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-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 258 total aligned reads displayed.

GCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGCC  >  CP009273/1003451‑1003743
                                                                                                                                                  |                                                                                                                                                  
  cgccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACAcc                                                                                                                                               <  1:29353/150‑1 (MQ=255)
  cgccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACAcc                                                                                                                                               <  1:2300218/150‑1 (MQ=255)
     cTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATc                                                                                                                                            >  1:2609059/1‑150 (MQ=255)
      tGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATcc                                                                                                                                           <  1:3723677/150‑1 (MQ=255)
      tGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATcc                                                                                                                                           >  2:584285/1‑150 (MQ=255)
         ccTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCa                                                                                                                                        >  2:2104209/1‑150 (MQ=255)
            cgcgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGTGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGAAAca                                                                                                                                     >  2:5145043/1‑150 (MQ=255)
            cgcgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAca                                                                                                                                     <  2:1528254/150‑1 (MQ=255)
            cgcgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAca                                                                                                                                     <  2:78451/150‑1 (MQ=255)
             gcgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcag                                                                                                                                    <  1:4925876/150‑1 (MQ=255)
             gcgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcag                                                                                                                                    <  2:3630878/150‑1 (MQ=255)
             gcgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcag                                                                                                                                    >  2:1188138/1‑150 (MQ=255)
                tattatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtca                                                                                                                                 <  2:271827/150‑1 (MQ=255)
                tattatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtca                                                                                                                                 <  1:4854937/150‑1 (MQ=255)
                tattatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtca                                                                                                                                 <  2:4346622/150‑1 (MQ=255)
                 attatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtcag                                                                                                                                >  2:568291/1‑150 (MQ=255)
                 attatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtcag                                                                                                                                >  1:1571089/1‑150 (MQ=255)
                  ttatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtcagt                                                                                                                               <  1:2383742/150‑1 (MQ=255)
                       ttGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGgc                                                                                                                          >  1:5036181/1‑150 (MQ=255)
                          ccGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGcgcg                                                                                                                       <  2:1971280/150‑1 (MQ=255)
                          ccGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGcgcg                                                                                                                       <  1:104264/150‑1 (MQ=255)
                           agCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGa                                                                                                                      >  1:2057707/2‑150 (MQ=255)
                             cTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtg                                                                                                                    <  1:2976952/150‑1 (MQ=255)
                              tGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgc                                                                                                                   >  2:3430893/1‑149 (MQ=255)
                              tGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtga                                                                                                                   >  2:4319091/1‑150 (MQ=255)
                              tGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtga                                                                                                                   <  2:3875598/150‑1 (MQ=255)
                                  cGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagt                                                                                                               <  2:407001/150‑1 (MQ=255)
                                  cGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagt                                                                                                               <  2:406849/150‑1 (MQ=255)
                                   gAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtg                                                                                                              >  2:305604/1‑150 (MQ=255)
                                      tgtAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaa                                                                                                           <  2:1694747/150‑1 (MQ=255)
                                      tgtAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaa                                                                                                           <  2:1556352/150‑1 (MQ=255)
                                        tAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGa                                                                                                         >  2:1268531/1‑150 (MQ=255)
                                        tAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGa                                                                                                         <  2:2394206/150‑1 (MQ=255)
                                         aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAc                                                                                                        >  2:3329962/1‑150 (MQ=255)
                                          aGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAcg                                                                                                       <  1:1206018/150‑1 (MQ=255)
                                               tACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtc                                                                                                  <  2:284650/150‑1 (MQ=255)
                                               tACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtc                                                                                                  >  1:4310146/1‑150 (MQ=255)
                                               tACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtc                                                                                                  <  2:1372360/150‑1 (MQ=255)
                                                aCAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtcg                                                                                                 >  2:2414381/1‑150 (MQ=255)
                                                    cGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGa                                                                                             <  1:3134552/150‑1 (MQ=255)
                                                       tttAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGAGGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGACAGACGCGATCGTCGATGc                                                                                          >  2:2236367/1‑150 (MQ=255)
                                                          aGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCttt                                                                                       <  2:3157771/150‑1 (MQ=255)
                                                            aCCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCa                                                                                     <  1:77948/150‑1 (MQ=255)
                                                              ctctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAcg                                                                                   <  1:1990309/150‑1 (MQ=255)
                                                                 tATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCgg                                                                                <  1:3228902/150‑1 (MQ=255)
                                                                 tATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCgg                                                                                <  1:1924808/150‑1 (MQ=255)
                                                                   tctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaa                                                                              <  1:3057869/150‑1 (MQ=255)
                                                                    ctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaa                                                                             >  1:2088758/1‑150 (MQ=255)
                                                                    ctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaa                                                                             <  1:4810072/150‑1 (MQ=255)
                                                                    ctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaa                                                                             <  2:3569608/150‑1 (MQ=255)
                                                                        gtgtTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATc                                                                         <  1:542299/150‑1 (MQ=255)
                                                                             cAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcg                                                                    >  1:2734635/1‑150 (MQ=255)
                                                                              aGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcgc                                                                   >  2:1853630/1‑150 (MQ=255)
                                                                              aGGCGAACAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGAGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAAACTGAcgc                                                                   >  2:2695213/1‑150 (MQ=255)
                                                                               ggCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcgcg                                                                  <  2:4841317/150‑1 (MQ=255)
                                                                                  gATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTcc                                                               >  1:1182510/1‑150 (MQ=255)
                                                                                     cAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCaaa                                                            >  1:3568048/1‑150 (MQ=255)
                                                                                     cAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCaaa                                                            >  1:3393320/1‑150 (MQ=255)
                                                                                     cAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCaaa                                                            >  2:2348544/1‑150 (MQ=255)
                                                                                        cTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGt                                                         >  1:3929208/1‑150 (MQ=255)
                                                                                            aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGAt                                                     >  1:4482164/1‑150 (MQ=255)
                                                                                            aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGAt                                                     <  1:4919190/150‑1 (MQ=255)
                                                                                             cAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATc                                                    >  1:3779399/1‑150 (MQ=255)
                                                                                             cAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATc                                                    <  1:937532/150‑1 (MQ=255)
                                                                                              agagATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAATTCTGCCGCGTCCAAATGTTGATcg                                                   <  1:174646/150‑1 (MQ=255)
                                                                                               gagaTGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATcgc                                                  <  1:4129213/150‑1 (MQ=255)
                                                                                                   tGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATg                                              <  1:4589579/150‑1 (MQ=255)
                                                                                                     tttAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATgcg                                            <  2:2922716/150‑1 (MQ=255)
                                                                                                         aTCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCgg                                        <  2:1990968/150‑1 (MQ=255)
                                                                                                         aTCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCgg                                        <  1:2486419/150‑1 (MQ=255)
                                                                                                          tCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAACGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGa                                       <  2:928865/150‑1 (MQ=255)
                                                                                                           ccTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGat                                      >  1:4386195/1‑150 (MQ=255)
                                                                                                           ccTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATAGCGATGCGCCGGat                                      >  2:1071644/1‑150 (MQ=255)
                                                                                                            cTGGCGCGTCCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGata                                     >  1:1638378/1‑150 (MQ=255)
                                                                                                            cTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGata                                     >  2:5106566/1‑150 (MQ=255)
                                                                                                             tGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGatat                                    >  2:4934218/1‑150 (MQ=255)
                                                                                                              ggCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATc                                   >  2:1227825/1‑150 (MQ=255)
                                                                                                                   gACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGt                              >  2:4723232/1‑150 (MQ=255)
                                                                                                                      cTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTaa                           >  2:3347515/1‑150 (MQ=255)
                                                                                                                      cTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTaa                           >  2:2573903/1‑150 (MQ=255)
                                                                                                                       ttCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAAc                          <  2:4666613/150‑1 (MQ=255)
                                                                                                                       ttCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAAc                          >  2:1946123/1‑150 (MQ=255)
                                                                                                                         cGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGt                        >  2:944133/1‑150 (MQ=255)
                                                                                                                         cGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGt                        >  1:1360901/1‑150 (MQ=255)
                                                                                                                         cGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGt                        >  2:4387385/1‑150 (MQ=255)
                                                                                                                            tGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTg                     <  2:5049517/150‑1 (MQ=255)
                                                                                                                            tGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTg                     >  2:2811615/1‑150 (MQ=255)
                                                                                                                            tGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTg                     >  2:4224374/1‑150 (MQ=255)
                                                                                                                             gTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTgg                    <  2:4538700/150‑1 (MQ=255)
                                                                                                                                cACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTg                 <  2:488396/150‑1 (MQ=255)
                                                                                                                                    tCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATa             <  1:555278/150‑1 (MQ=255)
                                                                                                                                     cAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaa            <  1:3382916/150‑1 (MQ=255)
                                                                                                                                         ggTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa        <  2:1080407/150‑1 (MQ=255)
                                                                                                                                         ggTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa        <  2:1314227/150‑1 (MQ=255)
                                                                                                                                         ggTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa        <  2:1570407/150‑1 (MQ=255)
                                                                                                                                            ttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAAcc     <  1:3607008/150‑1 (MQ=255)
                                                                                                                                            ttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAAcc     <  2:3963430/150‑1 (MQ=255)
                                                                                                                                            ttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAAcc     >  2:2974482/1‑150 (MQ=255)
                                                                                                                                             tGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCg    <  2:4796849/150‑1 (MQ=255)
                                                                                                                                               aaTTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTACACGCGGAAAAATCTGACGCGTCCAAAAGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGACACCGcc  >  2:4411531/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGCC  >  CP009273/1003451‑1003743

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: