Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,239,109 A→C 5.1% intergenic (‑78/‑297) argF ← / → PP_1080 ornithine carbamoyltransferase/molybdopterin‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,239,1090AC 5.1% 223.8 / 3.8 79intergenic (‑78/‑297)argF/PP_1080ornithine carbamoyltransferase/molybdopterin‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (37/38);  new base C (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

TCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCTAC  >  NC_002947/1238971‑1239254
                                                                                                                                          |                                                                                                                                                 
tccccgcGAGCAATGCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACgggg                                                                                                                                        <  2:243200/144‑1 (MQ=255)
 cACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGc                                                                                                                                       <  1:847134/150‑1 (MQ=255)
   ccccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAg                                                                                                                                                  >  1:825595/1‑137 (MQ=255)
   ccccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAg                                                                                                                                                  <  2:825595/137‑1 (MQ=255)
   ccccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCGCGGGGCgg                                                                                                                                     >  2:637434/1‑150 (MQ=255)
   ccccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCgg                                                                                                                                     <  1:702409/150‑1 (MQ=255)
   ccccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCgg                                                                                                                                     >  1:1921892/1‑150 (MQ=255)
    cccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGa                                                                                                                                    >  2:809584/1‑150 (MQ=255)
      cGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACTGGGGTGGCCTGCGGCGTAAGCCGCAAGGGGCGGACg                                                                                                                                  >  1:877981/1‑150 (MQ=255)
       gAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGa                                                                                                                                 >  1:744561/1‑150 (MQ=255)
          cAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAAt                                                                                                                              >  1:1630871/1‑150 (MQ=255)
           aaTTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATa                                                                                                                             <  2:1559422/150‑1 (MQ=255)
            aTTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAg                                                                                                                            >  2:1082937/1‑150 (MQ=255)
            aTTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAg                                                                                                                            <  1:1082937/150‑1 (MQ=255)
             ttCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACGTGCGGCTTACGCCGCa                                                                                                                                             >  1:1602897/1‑132 (MQ=255)
             ttCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACGTGCGGCTTACGCCGCa                                                                                                                                             >  2:1897037/1‑132 (MQ=255)
             ttCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACGTGCGGCTTACGCCGCa                                                                                                                                             <  1:1897037/132‑1 (MQ=255)
             ttCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACGTGCGGCTTACGCCGCa                                                                                                                                             <  2:1602897/132‑1 (MQ=255)
                gtcgGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGAACTAGGCGGACGaaagtgtgac                                                                                                                        >  2:111003/1‑143 (MQ=255)
                        tgAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCaa                                                                                                                >  1:883185/1‑150 (MQ=255)
                          aaGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATa                                                                                                              <  1:762843/150‑1 (MQ=255)
                               cagcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAgg                                                                                                                           >  1:1836201/1‑132 (MQ=255)
                               cagcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAgg                                                                                                                           <  2:1836201/132‑1 (MQ=255)
                                agcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATaa                                                                                                        <  2:883185/150‑1 (MQ=255)
                                agcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATaa                                                                                                        <  1:809584/150‑1 (MQ=255)
                                    gggAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGAAACGGTATGCAGGTCGACCATATGGCGTGAGTGGACAACGGGGATGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTaa                                                                                                    >  1:1492919/1‑150 (MQ=255)
                                     ggAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTaaa                                                                                                   >  1:330432/1‑150 (MQ=255)
                                           aaaGCGCCTGGCGCGCATGAGTCGCTGCCTTAGCTGCAACGGTAGGCAGGGCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGc                                                                                             <  1:1567371/150‑1 (MQ=255)
                                                 ccTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAg                                                                                                                            >  2:1707135/1‑113 (MQ=255)
                                                 ccTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAg                                                                                                                            <  1:1707135/113‑1 (MQ=255)
                                                 ccTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAg                                                                                       >  2:1413330/1‑150 (MQ=255)
                                                   tAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCCGAGTACTATAAGTAAATGTCGCCTGAAGcc                                                                                     >  1:974789/1‑150 (MQ=255)
                                                         tCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAgggg                                                                               <  2:330432/150‑1 (MQ=255)
                                                         tCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAgggg                                                                               >  1:1821429/1‑150 (MQ=255)
                                                         tCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAgggg                                                                               >  2:1764891/1‑150 (MQ=255)
                                                          cATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAggggg                                                                              <  1:637434/150‑1 (MQ=255)
                                                             gAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGaa                                                                           <  2:1630871/150‑1 (MQ=255)
                                                                tGACCACCTTATGTGCCCCGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGCACTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAg                                                                        <  1:1529641/150‑1 (MQ=255)
                                                                    tACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAgcggc                                                                    <  2:744561/150‑1 (MQ=255)
                                                                    caCCTTCTCTGCAAGGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGAGGGGTAAGCCGCACGTGGCGGACGAAATAGGCAAAGGCGCAATACTATAAGTAAATGTCGCTTGAAGCCAAGGGGGGAAAAgcggc                                                                    <  1:867402/149‑1 (MQ=255)
                                                                      ccgtATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGt                                                                  >  1:1936353/4‑150 (MQ=255)
                                                                      ccTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATaa                                                                                                        >  1:504548/1‑112 (MQ=255)
                                                                      ccTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATaa                                                                                                        <  2:504548/112‑1 (MQ=255)
                                                                        ttATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAGAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGa                                                                >  2:1503850/1‑150 (MQ=255)
                                                                        ttATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGa                                                                >  2:1492158/1‑150 (MQ=255)
                                                                        ctATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATa                                                                                                                             >  1:1771280/2‑89 (MQ=255)
                                                                        ctATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATa                                                                                                                             <  2:1771280/88‑1 (MQ=255)
                                                                          aTCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAtt                                                              <  1:1764891/150‑1 (MQ=255)
                                                                          aTCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAtt                                                              <  2:865413/150‑1 (MQ=255)
                                                                             tGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgca                                                           >  2:1634131/1‑150 (MQ=255)
                                                                              gCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAAt                                                                                                                              <  1:1667848/82‑1 (MQ=255)
                                                                              gCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAAt                                                                                                                              >  2:1667848/1‑82 (MQ=255)
                                                                               cAACGGCATGGAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaa                                                         <  1:103872/150‑1 (MQ=255)
                                                                                aaCGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaag                                                        >  1:1588380/1‑150 (MQ=255)
                                                                                    ggaTGCAGGTCGCCCGGATGGTTTTTTTTGACAACTGGGGTGACCTGCGGCGTAAGCCGCCCGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGGCCCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtg                                                    <  2:877981/148‑1 (MQ=255)
                                                                                    gTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACtat                                                                                                          <  2:1549790/96‑1 (MQ=255)
                                                                                    gTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACtat                                                                                                          >  1:1549790/1‑96 (MQ=255)
                                                                                        gCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgc                                                <  2:1821429/150‑1 (MQ=255)
                                                                                          aGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAgggg                                                                               >  2:955908/1‑117 (MQ=255)
                                                                                          aGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAgggg                                                                               <  1:955908/117‑1 (MQ=255)
                                                                                             tCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgcaaagg                                           >  1:783946/1‑150 (MQ=255)
                                                                                                    tATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCtt                                    <  2:1492919/150‑1 (MQ=255)
                                                                                                       ggTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTAc                                 <  2:1921892/150‑1 (MQ=255)
                                                                                                       ggTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTAc                                 <  1:1712392/150‑1 (MQ=255)
                                                                                                         ttttttttGACAACTGTGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTAccc                               >  2:1090037/1‑150 (MQ=255)
                                                                                                           ttttttGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTc                                      <  1:152492/141‑1 (MQ=255)
                                                                                                           ttttttGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTc                                      >  2:152492/1‑141 (MQ=255)
                                                                                                              tttGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTgg                          >  1:62063/1‑150 (MQ=255)
                                                                                                                        ggggTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGt                <  1:1634131/150‑1 (MQ=255)
                                                                                                                            tGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAgggg                                                                               >  2:1504611/1‑83 (MQ=255)
                                                                                                                            tGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAgggg                                                                               <  1:1504611/83‑1 (MQ=255)
                                                                                                                              aCCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTa          <  1:1533344/150‑1 (MQ=255)
                                                                                                                               ccTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTAt         <  1:1503850/150‑1 (MQ=255)
                                                                                                                                cTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGt                                                                                                      >  1:113200/1‑56 (MQ=255)
                                                                                                                                cTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGt                                                                                                      <  2:113200/56‑1 (MQ=255)
                                                                                                                                 tGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTAttt       >  1:22152/1‑150 (MQ=255)
                                                                                                                                  gcggcgTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaa                                                         >  2:663205/1‑99 (MQ=255)
                                                                                                                                  gcggcgTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaa                                                         <  1:663205/99‑1 (MQ=255)
                                                                                                                                      cgTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCTAc  >  1:1045617/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
TCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTATGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTATTTGCTAC  >  NC_002947/1238971‑1239254

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: