Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,221,214 C→A 70.2% intergenic (+60/+73) rluF → / ← yjbD 23S rRNA pseudouridine(2604) synthase/DUF3811 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,221,2140CA70.2% 343.5 / inf 343intergenic (+60/+73)rluF/yjbD23S rRNA pseudouridine(2604) synthase/DUF3811 family protein
Reads supporting (aligned to +/- strand):  ref base C (47/55);  new base A (118/123);  total (165/178)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 354 total aligned reads displayed.

AAGATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGCTTTT  >  CP009273/4221065‑4221363
                                                                                                                                                     |                                                                                                                                                     
aaGATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTc                                                                                                                                                       >  1:1314753/1‑150 (MQ=255)
aaGATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTc                                                                                                                                                       <  2:4992363/150‑1 (MQ=255)
     ggAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGac                                                                                                                                                  >  2:608906/1‑149 (MQ=255)
     ggAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGaa                                                                                                                                                  >  2:3940849/1‑150 (MQ=255)
          aaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAggg                                                                                                                                             >  2:4647376/1‑150 (MQ=255)
          aaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAggg                                                                                                                                             >  1:1739535/1‑150 (MQ=255)
              ggcggAAACAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAg                                                                                                                                         <  1:4217395/150‑1 (MQ=255)
               gcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAgg                                                                                                                                        <  2:5120157/150‑1 (MQ=255)
                cggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAggg                                                                                                                                       <  1:2963361/150‑1 (MQ=255)
                 ggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGt                                                                                                                                      >  1:5165476/1‑150 (MQ=255)
                 ggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGt                                                                                                                                      >  2:4162146/1‑150 (MQ=255)
                  gAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACTCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTg                                                                                                                                     >  2:3494831/1‑150 (MQ=255)
                   aaaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGa                                                                                                                                    >  2:4674036/1‑150 (MQ=255)
                   aaaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGAAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGa                                                                                                                                    <  2:4762401/150‑1 (MQ=255)
                            gTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAAt                                                                                                                           <  1:706786/150‑1 (MQ=255)
                             tCGCCCGGCGTCCATCTGTTAGAGTTTTACCTCGCCGGGGAGTTAAAAGAAGGGGCGCTGATTAACGCCTCAGACTTATGAAAACGTCACAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGTTTTAGGGTGAGGGGAAAAtt                                                                                                                          <  2:3956973/150‑1 (MQ=37)
                             tCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAAtt                                                                                                                          <  2:517992/150‑1 (MQ=255)
                              cGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATtg                                                                                                                         <  2:5145914/150‑1 (MQ=255)
                               gCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATtgt                                                                                                                        >  2:558255/1‑150 (MQ=255)
                                 ccGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTgtgt                                                                                                                      >  2:3411042/1‑150 (MQ=255)
                                     cGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTcac                                                                                                                  <  1:294203/150‑1 (MQ=255)
                                          aaCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGa                                                                                                             <  1:340451/150‑1 (MQ=255)
                                           aCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGAc                                                                                                            <  1:4893315/150‑1 (MQ=255)
                                           aCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGAc                                                                                                            <  2:303785/150‑1 (MQ=255)
                                           aCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGAc                                                                                                            <  2:1036964/150‑1 (MQ=255)
                                            cGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACg                                                                                                           <  2:2404417/150‑1 (MQ=255)
                                            cGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACg                                                                                                           <  2:4522697/150‑1 (MQ=255)
                                                   cGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGaaa                                                                                                    >  1:1126609/1‑150 (MQ=255)
                                                    gTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGaaaa                                                                                                   >  2:74066/1‑150 (MQ=255)
                                                        tACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtctc                                                                                               >  1:1915498/1‑150 (MQ=255)
                                                        tACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtctc                                                                                               <  1:2297649/150‑1 (MQ=255)
                                                        tACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtctc                                                                                               <  1:2011706/150‑1 (MQ=255)
                                                          ccTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAg                                                                                             <  1:4426697/150‑1 (MQ=255)
                                                          ccTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAg                                                                                             >  1:427994/1‑150 (MQ=255)
                                                           cTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAgg                                                                                            <  1:1176707/150‑1 (MQ=255)
                                                           cTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAgg                                                                                            >  1:2522310/1‑150 (MQ=255)
                                                                    gCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtg                                                                                   >  2:3653872/1‑150 (MQ=255)
                                                                     cGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGGTGATggt                                                                                  >  2:3529104/1‑150 (MQ=255)
                                                                     agTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtgt                                                                                  >  2:2786556/2‑150 (MQ=255)
                                                                      gttatatgaagGGGCGCTTATTCACGCCTCAGCCTGATGAACACGGACCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTATGGAAAATGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGAAGAGTGAAAATCTCAGGCGTTGATtgtg                                                                                 <  1:1442524/143‑1 (MQ=255)
                                                                       tAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTgg                                                                                >  2:172770/1‑150 (MQ=255)
                                                                       tAAAAAGAAGGGGAGCTGATTAACGCCTCAGACTTAAGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCATTAGGGAAAAGGGTTCGGGTGAGGGGAAAATTGTTTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTgg                                                                                <  1:4789604/150‑1 (MQ=255)
                                                                        aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAACACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAAGTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGa                                                                               <  2:4926199/150‑1 (MQ=255)
                                                                        aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGCAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGa                                                                               >  2:4205070/1‑150 (MQ=255)
                                                                        aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGa                                                                               >  1:521929/1‑150 (MQ=255)
                                                                         aaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGGTGATTGTGGAt                                                                              >  2:1828081/1‑150 (MQ=255)
                                                                          aaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGAtt                                                                             >  1:1196349/1‑150 (MQ=255)
                                                                           aagaagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTa                                                                            <  1:3495790/150‑1 (MQ=255)
                                                                           aagaagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTa                                                                            >  1:809119/1‑150 (MQ=255)
                                                                             gaagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAgc                                                                          >  1:4862002/1‑150 (MQ=255)
                                                                                 gggCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCt                                                                      <  1:3885732/150‑1 (MQ=255)
                                                                                 gggCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCt                                                                      <  2:3428948/150‑1 (MQ=255)
                                                                                  ggCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTg                                                                     >  1:4647191/1‑150 (MQ=255)
                                                                                    cgcTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGcc                                                                   <  1:1795035/150‑1 (MQ=255)
                                                                                           tAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTc                                                            <  1:4552373/150‑1 (MQ=255)
                                                                                           tAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTc                                                            >  2:3862935/1‑150 (MQ=255)
                                                                                           tAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTc                                                            <  1:354156/150‑1 (MQ=255)
                                                                                           tAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTc                                                            <  2:4338036/150‑1 (MQ=255)
                                                                                            aaCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCg                                                           <  1:3697475/150‑1 (MQ=255)
                                                                                             aCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAATTGTGTTCACACTGATGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGAc                                                         <  1:2106054/150‑1 (MQ=255)
                                                                                              cGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGAc                                                         <  1:3111055/150‑1 (MQ=255)
                                                                                               gCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACg                                                        >  2:1531938/1‑150 (MQ=255)
                                                                                                    aGACTGCTGAAAACGGCCCAACAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAttt                                                   <  2:2583570/150‑1 (MQ=255)
                                                                                                    aGACTGATGACAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAttt                                                   <  2:327554/150‑1 (MQ=255)
                                                                                                    aGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAttt                                                   <  1:3210673/150‑1 (MQ=255)
                                                                                                    aGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAttt                                                   <  1:1504939/150‑1 (MQ=255)
                                                                                                    aGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAttt                                                   <  1:1104637/150‑1 (MQ=255)
                                                                                                     gACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTg                                                  <  2:3735528/150‑1 (MQ=255)
                                                                                                      aCTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGc                                                 >  2:2560115/1‑150 (MQ=255)
                                                                                                      aCTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGc                                                 >  2:1152078/1‑150 (MQ=255)
                                                                                                       cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCt                                                <  1:2378487/150‑1 (MQ=255)
                                                                                                       cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCt                                                <  1:2420489/150‑1 (MQ=255)
                                                                                                        tgatgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                               >  1:8807/1‑150 (MQ=255)
                                                                                                             aaaaCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAg                                          <  2:809119/150‑1 (MQ=255)
                                                                                                               aaCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCt                                        >  2:4383285/1‑150 (MQ=255)
                                                                                                               aaCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCt                                        >  1:1006192/1‑150 (MQ=255)
                                                                                                                aCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTa                                       <  1:74066/150‑1 (MQ=255)
                                                                                                                    cccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATg                                   <  2:2875961/150‑1 (MQ=255)
                                                                                                                    cccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATg                                   >  2:2310293/1‑150 (MQ=255)
                                                                                                                    cccAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATg                                   <  1:3351805/150‑1 (MQ=255)
                                                                                                                       aaaaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACg                                <  2:917549/150‑1 (MQ=255)
                                                                                                                       aaaaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACg                                <  1:4017166/150‑1 (MQ=255)
                                                                                                                       aaaaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACg                                <  1:4162146/150‑1 (MQ=255)
                                                                                                                       aaaaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACg                                >  2:1825825/1‑150 (MQ=255)
                                                                                                                        aaaGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTCAATGACGc                               >  2:329524/1‑150 (MQ=255)
                                                                                                                           gTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTc                            >  1:3922257/1‑150 (MQ=255)
                                                                                                                           gTGCCCGGACAGTCCCCTCGCCCCTTAGGGATAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTc                            <  2:1640106/150‑1 (MQ=255)
                                                                                                                           gTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTc                            <  2:4103482/150‑1 (MQ=255)
                                                                                                                             gCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAg                          >  2:1748465/1‑150 (MQ=255)
                                                                                                                                 ggACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATc                      >  1:2779347/1‑150 (MQ=255)
                                                                                                                                      gTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGtt                 <  2:3207561/150‑1 (MQ=255)
                                                                                                                                          ccTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTtata             <  2:2196624/150‑1 (MQ=255)
                                                                                                                                            tCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAg           <  1:2656326/150‑1 (MQ=255)
                                                                                                                                             cGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGc          >  1:588222/1‑150 (MQ=255)
                                                                                                                                              gCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGcc         <  2:2700821/150‑1 (MQ=255)
                                                                                                                                               ccccTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCt        <  1:684867/150‑1 (MQ=255)
                                                                                                                                               ccccTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCt        >  2:2080216/1‑150 (MQ=255)
                                                                                                                                                cccTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTg       <  2:1366671/150‑1 (MQ=255)
                                                                                                                                                cccTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTg       <  2:5165557/150‑1 (MQ=255)
                                                                                                                                                  cTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGCt     <  1:4511770/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGCTTTT  >  CP009273/4221065‑4221363

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: