Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,355,833 G→C 8.7% S85C (TCC→TGC) ‡ pnp ← polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,355,8330GC 8.7% 557.0 / 33.5 172S85C (TCC→TGC) ‡pnppolyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (75/82);  new base C (8/7);  total (83/89)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 176 total aligned reads displayed.

AGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTT  >  NC_002947/5355684‑5355981
                                                                                                                                                     |                                                                                                                                                    
aGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCgg                                                                                                                                                      >  1:1136853/1‑150 (MQ=255)
  gTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGaa                                                                                                                                                    <  2:212854/150‑1 (MQ=255)
  gTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGaa                                                                                                                                                    <  1:1017350/150‑1 (MQ=255)
   tACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAg                                                                                                                                                   >  2:1882390/1‑150 (MQ=255)
        atcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACgg                                                                                                                                              >  1:3091574/1‑150 (MQ=255)
        atcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACgg                                                                                                                                              >  1:3693730/1‑150 (MQ=255)
         tcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGc                                                                                                                                             >  1:1417326/1‑150 (MQ=255)
         tcatcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGc                                                                                                                                             <  2:3867582/150‑1 (MQ=255)
          catcGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGcc                                                                                                                                            <  1:1081854/150‑1 (MQ=255)
              gCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCa                                                                                                                                        >  1:3128301/1‑150 (MQ=255)
              gCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCa                                                                                                                                        >  2:1242491/1‑150 (MQ=255)
                      gTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaa                                                                                                                                >  1:1434722/1‑150 (MQ=255)
                      gTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaa                                                                                                                                <  1:1618250/150‑1 (MQ=255)
                      gTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaa                                                                                                                                <  1:3169673/150‑1 (MQ=255)
                         cGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaagaa                                                                                                                             <  2:1909683/150‑1 (MQ=255)
                         cGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTgaagaa                                                                                                                             <  2:1647934/150‑1 (MQ=255)
                                   ttcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTggg                                                                                                                   >  1:738126/1‑150 (MQ=255)
                                   ttcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTggg                                                                                                                   >  2:572996/1‑150 (MQ=255)
                                    tcttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGa                                                                                                                  <  2:3219146/150‑1 (MQ=255)
                                     cttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGAt                                                                                                                 <  1:3813570/150‑1 (MQ=255)
                                     cttGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGAt                                                                                                                 <  1:4205614/150‑1 (MQ=255)
                                        gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                                                                              <  2:799797/150‑1 (MQ=255)
                                        gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                                                                              <  1:2200851/150‑1 (MQ=255)
                                        gCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCtt                                                                                                              <  2:462357/150‑1 (MQ=255)
                                         cTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTg                                                                                                             >  1:4128899/1‑150 (MQ=255)
                                            gtggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCa                                                                                                          <  1:1961384/150‑1 (MQ=255)
                                              ggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgc                                                                                                        >  1:360802/1‑150 (MQ=255)
                                              ggAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgc                                                                                                        >  2:1056148/1‑150 (MQ=255)
                                               gAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgcg                                                                                                       >  1:4934414/1‑150 (MQ=255)
                                                 aaCCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAgcggc                                                                                                     >  1:3962242/1‑150 (MQ=255)
                                                    cACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTa                                                                                                  >  1:1322581/1‑150 (MQ=255)
                                                        gTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTc                                                                                              >  2:4497782/1‑150 (MQ=255)
                                                             gacgacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTc                                                                                         <  2:2200317/150‑1 (MQ=255)
                                                             gacgacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTc                                                                                         <  1:3521915/150‑1 (MQ=255)
                                                              acgacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTcc                                                                                        <  2:2062010/150‑1 (MQ=255)
                                                                gacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTg                                                                                      >  1:905530/1‑150 (MQ=255)
                                                                gacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTg                                                                                      <  1:4187523/150‑1 (MQ=255)
                                                                gacCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTg                                                                                      <  1:4178532/150‑1 (MQ=255)
                                                                 acCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTgg                                                                                     <  2:2271758/150‑1 (MQ=255)
                                                                 acCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTgg                                                                                     <  1:4610957/150‑1 (MQ=255)
                                                                   ctgcgtgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATa                                                                                   <  2:3807330/143‑1 (MQ=255)
                                                                   ctgcgtgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATa                                                                                   <  1:2971208/143‑1 (MQ=255)
                                                                   cTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTa                                                                                   <  1:3976441/150‑1 (MQ=255)
                                                                   cTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTa                                                                                   <  1:1730050/150‑1 (MQ=255)
                                                                   cTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTa                                                                                   <  1:481119/150‑1 (MQ=255)
                                                                    tGCACTTCGTTCATGAAGCCATCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGACAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGACACCTGGGATCTTGCCAGCGGCGAAGGTCTTTTCCTGGTAg                                                                                  >  2:4548324/1‑150 (MQ=255)
                                                                       gtgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGttg                                                                               <  1:498736/147‑3 (MQ=255)
                                                                        tgtCGTTCTTGTAGCCATCAACGAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCGGCGTAAGTCTTTTCCTGATAGttga                                                                              >  2:1582952/3‑147 (MQ=255)
                                                                        cTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGaa                                                                              <  2:3455844/150‑1 (MQ=255)
                                                                         ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGc                                                                                                                                             <  2:1148241/86‑1 (MQ=255)
                                                                         ttCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                             >  2:1301089/1‑150 (MQ=255)
                                                                         tgcGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAAc                                                                             >  1:2754749/3‑150 (MQ=255)
                                                                              tCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACa                                                                        >  1:1865252/1‑150 (MQ=255)
                                                                                aTGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGc                                                                      >  2:541292/1‑150 (MQ=255)
                                                                                  gAAGCCTTCAGGTAACAGCGGACGGATCTGACTGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCTGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCgg                                                                    <  2:3143473/150‑1 (MQ=255)
                                                                                     gCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCAc                                                                                                                                       <  1:3378489/80‑1 (MQ=255)
                                                                                        ttCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaagaa                                                              >  1:408538/1‑150 (MQ=255)
                                                                                        ttCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGgaagaa                                                              >  2:1425216/1‑150 (MQ=255)
                                                                                        atcaacgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGgaagaa                                                              >  1:2691825/7‑150 (MQ=255)
                                                                                          caacgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGaaat                                                            <  2:2932017/146‑2 (MQ=255)
                                                                                           cGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc                                                           <  2:3343106/150‑1 (MQ=255)
                                                                                           cGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc                                                           <  2:4019775/150‑1 (MQ=255)
                                                                                           cGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAAcc                                                           >  2:3310430/1‑150 (MQ=255)
                                                                                           aacgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGaaatc                                                           >  1:1860273/4‑148 (MQ=255)
                                                                                            acgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAAtcc                                                          >  2:5078149/3‑150 (MQ=255)
                                                                                            acgAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAAtcc                                                          >  2:1573333/3‑150 (MQ=255)
                                                                                              gAAAAGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAATCCtg                                                        >  1:1609826/1‑149 (MQ=255)
                                                                                                aCAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGc                                                      <  1:3740146/150‑1 (MQ=255)
                                                                                                 aaGCGGGCGGATCGGACGGTCGATCAGGCGCGAAACCAGAGTTTCGTTTTCGCTTGGACGACCTTCGCGCTTGAAGAAACCACCTGGGATCTTGCCAGCTGCGTAAGTCTTTTCCTGATAGTTGACGGTCAGCGGGAAGAAATCCtgacc                                                     >  2:1391405/2‑146 (MQ=255)
                                                                                                  aGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGCCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCGTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCt                                                    <  1:4497803/150‑1 (MQ=255)
                                                                                                    cggacggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTgg                                                  >  1:4676008/1‑150 (MQ=255)
                                                                                                     ggacggaTCGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGAt                                                >  1:966140/1‑150 (MQ=255)
                                                                                                      gacggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGACGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGAt                                                <  2:2055718/150‑1 (MQ=255)
                                                                                                       acggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATc                                               >  2:4543540/1‑150 (MQ=255)
                                                                                                       acggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATc                                               <  1:3518994/150‑1 (MQ=255)
                                                                                                         ggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAg                                             >  2:4388522/1‑150 (MQ=255)
                                                                                                         ggaTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAg                                             >  2:4388422/1‑150 (MQ=255)
                                                                                                           aTCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGcc                                           >  2:2426543/1‑150 (MQ=255)
                                                                                                               gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGtt                                       <  1:216626/150‑1 (MQ=255)
                                                                                                               gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGtt                                       >  1:3719564/1‑150 (MQ=255)
                                                                                                               gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGtt                                       >  2:3641957/1‑150 (MQ=255)
                                                                                                               gACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGttt                                      <  2:336283/150‑1 (MQ=255)
                                                                                                                    tCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGgcg                                  >  1:806686/1‑150 (MQ=255)
                                                                                                                      gATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGcc                                >  1:2613649/1‑150 (MQ=255)
                                                                                                                      gATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGcc                                >  1:170731/1‑150 (MQ=255)
                                                                                                                          aGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTAcc                            <  1:4235511/150‑1 (MQ=255)
                                                                                                                             cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                         <  1:2316541/150‑1 (MQ=255)
                                                                                                                             cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                         <  2:4022278/150‑1 (MQ=255)
                                                                                                                             cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                         <  1:630909/150‑1 (MQ=255)
                                                                                                                             cgcgAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                         <  1:3831709/150‑1 (MQ=255)
                                                                                                                             cgcgAGGTCAGCGTCTCGTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACcaca                         <  1:2737410/150‑1 (MQ=255)
                                                                                                                                gAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTc                      >  1:1024814/1‑150 (MQ=255)
                                                                                                                                      aGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGc                >  1:1985135/1‑150 (MQ=255)
                                                                                                                                       gCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCa               >  2:4326571/1‑150 (MQ=255)
                                                                                                                                           ctctTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGt           <  1:2584887/150‑1 (MQ=255)
                                                                                                                                              tttCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAc        >  2:4893823/1‑150 (MQ=255)
                                                                                                                                                tctcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAcgt      >  1:1239912/1‑150 (MQ=255)
                                                                                                                                                 ctcGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGAcgtc     <  2:4934414/150‑1 (MQ=255)
                                                                                                                                                    ggAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGtt  <  2:2797584/150‑1 (MQ=255)
                                                                                                                                                    ggAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTGTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGtt  <  2:3091574/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGTACCGATCATCGCAGCGATGTCCGGGTCGGTCTTCTTGCTGGTGGAAACCACGGTGCAGACGACCTGCACTTCGTTCATGAAGCCTTCCGGGAACAGCGGACGGATCGGACGGTCGATCAGGCGCGAGGTCAGCGTCTCTTTCTCGGAAGGACGGCCTTCACGCTTGAAGAAGCCACCTGGGATCTTGCCAGCGGCGTAGGTCTTTTCCTGGTAGTGAACCGACAGCGGGAAGAAACCCTTGCCTGGATCAGCCTGTTTGGCGCCTACCACAGTCACCAGCACGGTGACGTCGTT  >  NC_002947/5355684‑5355981

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: