Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,221,214 C→A 93.4% 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,2140CA93.4% 817.9 / 40.5 293intergenic (+60/+73)rluF/yjbD23S rRNA pseudouridine(2604) synthase/DUF3811 family protein
Reads supporting (aligned to +/- strand):  ref base C (10/10);  new base A (128/145);  total (138/155)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 294 total aligned reads displayed.

AAGATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGCT  >  CP009273/4221065‑4221360
                                                                                                                                                     |                                                                                                                                                  
    tGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGa                                                                                                                                                >  2:2202082/1‑150 (MQ=255)
     ggAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGaa                                                                                                                                               <  1:3733405/150‑1 (MQ=255)
          aaaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAggg                                                                                                                                          >  2:4950715/1‑150 (MQ=255)
           aaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGt                                                                                                                                         >  1:389651/1‑150 (MQ=255)
           aaCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGt                                                                                                                                         >  1:2848473/1‑150 (MQ=255)
            aCGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGtt                                                                                                                                        >  1:2481699/1‑150 (MQ=255)
             cggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTa                                                                                                                                       <  2:3033756/150‑1 (MQ=255)
              ggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAg                                                                                                                                      <  1:874439/150‑1 (MQ=255)
              ggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAg                                                                                                                                      <  2:2814724/150‑1 (MQ=255)
              ggcggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAg                                                                                                                                      <  2:2929551/150‑1 (MQ=255)
                cggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAggg                                                                                                                                    >  1:2378019/1‑150 (MQ=255)
                 ggAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGt                                                                                                                                   >  1:1066514/1‑150 (MQ=255)
                    aaaaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCGGAGGGAACAGggtgagggtgag                                                                                                                                >  2:910821/1‑150 (MQ=255)
                      aaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAggg                                                                                                                              <  2:4250314/150‑1 (MQ=255)
                      aaGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAggg                                                                                                                              >  1:1296129/1‑150 (MQ=255)
                       aGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAgggg                                                                                                                             <  1:1935929/150‑1 (MQ=255)
                          cGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGaaa                                                                                                                          >  2:4052617/1‑150 (MQ=255)
                               gCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGACAATtgt                                                                                                                     >  1:3348456/1‑150 (MQ=255)
                               gCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATtgt                                                                                                                     >  2:3568344/1‑150 (MQ=255)
                               gCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATtgt                                                                                                                     >  2:1931590/1‑150 (MQ=255)
                               gCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATtgt                                                                                                                     >  1:396903/1‑150 (MQ=255)
                                cccGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATtgtg                                                                                                                    <  1:3834739/150‑1 (MQ=255)
                                 ccGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTgtgt                                                                                                                   <  2:439165/150‑1 (MQ=255)
                                 ccGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTgtgt                                                                                                                   >  2:4788592/1‑150 (MQ=255)
                                    gCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTca                                                                                                                <  2:3530851/150‑1 (MQ=255)
                                     cGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTcac                                                                                                               >  2:2697635/1‑150 (MQ=255)
                                      gTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTcaca                                                                                                              >  2:457472/1‑150 (MQ=255)
                                      gTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTcaca                                                                                                              <  2:3971617/150‑1 (MQ=255)
                                         cAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTg                                                                                                           >  2:2966053/1‑150 (MQ=255)
                                         cAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTg                                                                                                           <  1:4209728/150‑1 (MQ=255)
                                         cAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTg                                                                                                           <  1:4570648/150‑1 (MQ=255)
                                          aaCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGa                                                                                                          <  2:4176944/150‑1 (MQ=255)
                                          aaCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGa                                                                                                          <  2:1598920/150‑1 (MQ=255)
                                          aaCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGa                                                                                                          >  1:723736/1‑150 (MQ=255)
                                           aCGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGAc                                                                                                         >  2:4723156/1‑150 (MQ=255)
                                              gTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAg                                                                                                      >  2:1380565/1‑150 (MQ=255)
                                               tAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGt                                                                                                     <  2:2368879/150‑1 (MQ=255)
                                                     ttttACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGACAAGGGTTAGGGTGAGGGGAAAAGTGTGTTCACACTGACGAGTGAAAAt                                                                                               >  2:288216/1‑150 (MQ=255)
                                                       ttACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtct                                                                                             <  1:2506181/150‑1 (MQ=255)
                                                        tACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAAtctc                                                                                            >  1:1047302/1‑150 (MQ=255)
                                                          ccTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAg                                                                                          <  2:529506/150‑1 (MQ=255)
                                                            tCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGc                                                                                        <  1:210714/150‑1 (MQ=255)
                                                                cGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGACAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTg                                                                                    >  1:3400834/1‑150 (MQ=255)
                                                                cGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTg                                                                                    >  2:3727155/1‑150 (MQ=255)
                                                                     cGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGTGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtgt                                                                               >  1:4431223/1‑150 (MQ=255)
                                                                     cGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATtgt                                                                               <  2:2999264/150‑1 (MQ=255)
                                                                        aaaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGa                                                                            <  2:41794/150‑1 (MQ=255)
                                                                         aaaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGAt                                                                           >  1:1630548/1‑150 (MQ=255)
                                                                          aaaGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGAtt                                                                          <  1:1503497/150‑1 (MQ=255)
                                                                           aagaagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTa                                                                         <  2:4088996/150‑1 (MQ=255)
                                                                              aagGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAgcg                                                                      >  1:1001274/1‑150 (MQ=255)
                                                                                   gcgcTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGc                                                                 <  1:1183428/150‑1 (MQ=255)
                                                                                     gcTGATTAACGCCTCAGACTGATGAAAACGGCCCAANAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcg                                                               >  1:2674378/1‑150 (MQ=255)
                                                                                       tGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcg                                                             >  1:4826528/1‑150 (MQ=255)
                                                                                         aTTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCcgcgcg                                                           >  2:1997300/1‑150 (MQ=255)
                                                                                            aaCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCg                                                        >  2:1876681/1‑150 (MQ=255)
                                                                                            aaCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCg                                                        <  1:4654135/150‑1 (MQ=255)
                                                                                            aaCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCg                                                        >  2:1092460/1‑150 (MQ=255)
                                                                                            aaCGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCg                                                        >  2:3672735/1‑150 (MQ=255)
                                                                                              cGCCTCAGACTAATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGAc                                                      <  2:2328993/150‑1 (MQ=255)
                                                                                                ccTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGt                                                    >  1:2996268/1‑150 (MQ=255)
                                                                                                   cAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAtt                                                 >  1:403517/1‑150 (MQ=255)
                                                                                                   cAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTAtt                                                 >  2:1476664/1‑150 (MQ=255)
                                                                                                      aCTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGc                                              <  2:116696/150‑1 (MQ=255)
                                                                                                       cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCt                                             <  2:2569187/150‑1 (MQ=255)
                                                                                                       cTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCt                                             <  2:3123473/150‑1 (MQ=255)
                                                                                                        tgatgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                            >  2:4881409/1‑150 (MQ=255)
                                                                                                        tgatgaAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTg                                            >  1:17499/1‑150 (MQ=255)
                                                                                                            gaAAACGGCCCAACAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCa                                        <  2:4431221/150‑1 (MQ=255)
                                                                                                             aaaaCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAg                                       <  2:4194642/150‑1 (MQ=255)
                                                                                                              aaaCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGc                                      <  2:4972563/150‑1 (MQ=255)
                                                                                                                aCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTa                                    <  1:2057484/150‑1 (MQ=255)
                                                                                                                aCGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTa                                    <  1:5021010/150‑1 (MQ=255)
                                                                                                                  ggCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTaaa                                  >  2:4915370/1‑150 (MQ=255)
                                                                                                                  ggCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGGTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTaaa                                  >  2:3906291/1‑150 (MQ=255)
                                                                                                                       aaaaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACg                             >  2:760169/1‑150 (MQ=255)
                                                                                                                        aaaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGc                            >  2:640380/1‑150 (MQ=255)
                                                                                                                         aaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCt                           <  2:1260201/150‑1 (MQ=255)
                                                                                                                         aaGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAAATGTGTTCACACAGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGACGACGTATTTGCTGACCAGCTAAATGACGCt                           >  1:4718529/1‑150 (MQ=255)
                                                                                                                          aGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCtt                          >  1:3566219/1‑150 (MQ=255)
                                                                                                                          aGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCtt                          <  2:2481699/150‑1 (MQ=255)
                                                                                                                          aGTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCtt                          <  1:4121991/150‑1 (MQ=255)
                                                                                                                           gTGCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTc                         <  1:344399/150‑1 (MQ=255)
                                                                                                                           gTGCCCGGACAGTCCCCTCGCACCTTAGTGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACTAGTGAAAATCTCAGGAGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCGGACCAGCTAAATGACGCTTc                         <  2:3078977/150‑1 (MQ=255)
                                                                                                                             gCCCGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAg                       <  1:559273/150‑1 (MQ=255)
                                                                                                                              cccGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAgg                      <  2:3009574/150‑1 (MQ=255)
                                                                                                                               ccGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGa                     >  2:2147166/1‑150 (MQ=255)
                                                                                                                                cGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGAt                    >  2:3430654/1‑150 (MQ=255)
                                                                                                                                cGGACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGAt                    >  1:4357754/1‑150 (MQ=255)
                                                                                                                                 ggACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATc                   <  2:1066488/150‑1 (MQ=255)
                                                                                                                                 ggACAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATc                   >  2:2454367/1‑150 (MQ=255)
                                                                                                                                   aCAGTCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCg                 >  2:33893/1‑150 (MQ=255)
                                                                                                                                       tCCCCTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGttt             <  1:4723162/150‑1 (MQ=255)
                                                                                                                                        ccccTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTta            >  1:3254284/1‑150 (MQ=255)
                                                                                                                                        ccccTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTta            >  1:3143178/1‑150 (MQ=255)
                                                                                                                                         cccTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTtat           <  2:3275074/150‑1 (MQ=255)
                                                                                                                                           cTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATaa         >  1:1713246/1‑150 (MQ=255)
                                                                                                                                           cTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATaa         <  1:2221336/150‑1 (MQ=255)
                                                                                                                                           cTCGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATaa         >  2:1689986/1‑150 (MQ=255)
                                                                                                                                             cGCCCCTTAGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGc       >  1:900407/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAGATGGAAAAAACGGCGGAAAAAGGCGGTCGCCCGGCGTCCAACGGTAAGCGTTTTACCTCGCCGGGGCGTAAAAAGAAGGGGCGCTGATTAACGCCTCAGACTGATGAAAACGGCCCAAAAGTGCCCGGACAGTCCCCTCGCCCCTTCGGGAAAAGGGTTAGGGTGAGGGGAAAATTGTGTTCACACTGACGAGTGAAAATCTCAGGCGTTGATTGTGGATTAGCGCCTGCCGCGCGTCGACGTATTTGCTGACCAGCTAAATGACGCTTCAGGATCCGGTTTATAAGCCTGCT  >  CP009273/4221065‑4221360

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: