Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,645,196 G→A 59.0% intergenic (‑14/‑189) prlC ← / → rlmJ oligopeptidase A/23S rRNA m(6)A2030 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,645,1960GA59.0% 31.1 / 99.5 88intergenic (‑14/‑189)prlC/rlmJoligopeptidase A/23S rRNA m(6)A2030 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (14/22);  new base A (24/28);  total (38/50)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

GCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTACGCTTTCGCG  >  NC_000913/3645097‑3645295
                                                                                                   |                                                                                                   
gCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTgg                                                                                                     <  2:2276748/100‑1 (MQ=255)
  cTTAGTCACGGCTGGAACGACATGTTCCAGGAGAATTTCAGAAAACGGAGGTAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGAtt                                                                                                   <  1:1790290/100‑1 (MQ=255)
    tAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTaa                                                                                                 <  1:2287679/100‑1 (MQ=255)
      gTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAag                                                                                               <  2:876657/100‑1 (MQ=255)
      gTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAag                                                                                               >  2:824255/1‑100 (MQ=255)
        cACGGCTGGAACGACATGTTCCGGTACAACTTTAGTAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCGCGCGCGCAGCTCTGATTaaa                                                                                                >  1:1058338/1‑97 (MQ=255)
         aCGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAgg                                                                                            <  1:1395116/100‑1 (MQ=255)
          cGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGt                                                                                           <  2:47285/100‑1 (MQ=255)
          cGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGt                                                                                           <  1:1572070/100‑1 (MQ=255)
              tGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAg                                                                                       <  2:1512092/100‑1 (MQ=255)
               ggAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAgg                                                                                      <  1:1765157/100‑1 (MQ=255)
                 aaCGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGaa                                                                                    >  2:987630/1‑100 (MQ=255)
                 aaCGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGaa                                                                                    >  2:1087140/1‑100 (MQ=255)
                 aaCGACATGTACCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGaa                                                                                    >  1:544429/1‑100 (MQ=255)
                 aaCCACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGa                                                                                     <  2:213950/99‑1 (MQ=255)
                   cGACAGGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAgc                                                                                  <  1:502631/100‑1 (MQ=255)
                     aCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAgcgc                                                                                >  1:670812/1‑100 (MQ=255)
                      cATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCt                                                                               >  2:507609/1‑100 (MQ=255)
                        tGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCa                                                                             <  1:1911839/100‑1 (MQ=255)
                         gTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCaa                                                                            >  2:37675/1‑100 (MQ=255)
                           tCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACa                                                                          >  1:1954973/1‑100 (MQ=255)
                             cGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAgg                                                                        >  1:2587919/1‑100 (MQ=255)
                              gggAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGgc                                                                       <  1:124923/100‑1 (MQ=255)
                               ggAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGgcg                                                                      <  2:1176158/100‑1 (MQ=255)
                               ggAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGgcg                                                                      >  2:1560279/1‑100 (MQ=255)
                               ggAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGgcg                                                                      >  1:195450/1‑100 (MQ=255)
                                 agaATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCt                                                                    <  2:2418822/100‑1 (MQ=255)
                                    aTTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTa                                                                 >  2:57060/1‑100 (MQ=255)
                                    aTTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTa                                                                 >  2:53715/1‑100 (MQ=255)
                                    aTTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTa                                                                 >  2:2250520/1‑100 (MQ=255)
                                     ttttAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTAc                                                                >  1:1525134/1‑100 (MQ=255)
                                     ttttAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTACTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTAc                                                                >  2:2394677/1‑100 (MQ=255)
                                      tttAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACa                                                               <  1:2540365/100‑1 (MQ=255)
                                      tttAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGACCTGGTTAAAGAGGTTTAGGAAGCGCTCAACAGGAGCTTTACa                                                               <  1:116472/100‑1 (MQ=255)
                                       ttAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACAt                                                              <  2:1609346/100‑1 (MQ=255)
                                       ttAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACAt                                                              <  1:1404094/100‑1 (MQ=255)
                                         aGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATaa                                                            <  2:1133262/100‑1 (MQ=255)
                                         aGAAAACGGAGGCAAGTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATaa                                                            <  1:2469607/100‑1 (MQ=255)
                                          gAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAt                                                           >  2:1990230/1‑100 (MQ=255)
                                          gAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAt                                                           >  2:697707/1‑100 (MQ=255)
                                               cGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAtgtgtg                                                      <  1:55185/100‑1 (MQ=255)
                                               cGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAtgtgtg                                                      <  1:55239/100‑1 (MQ=255)
                                               cGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAtgtgtg                                                      <  1:606682/100‑1 (MQ=255)
                                                   ggCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGc                                                  >  1:326241/1‑100 (MQ=255)
                                                   ggCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGc                                                  <  2:1430055/100‑1 (MQ=255)
                                                   ggCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGc                                                  <  2:1616788/100‑1 (MQ=255)
                                                     cAATTCAAAGGGATTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAGAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCAt                                                >  2:1524656/1‑100 (MQ=255)
                                                     cAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCAt                                                >  1:1828743/1‑100 (MQ=255)
                                                     cAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCAt                                                >  2:792919/1‑100 (MQ=255)
                                                        ttCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATggg                                             >  1:235217/1‑100 (MQ=255)
                                                        ttCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATggg                                             <  1:140790/100‑1 (MQ=255)
                                                          cAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGt                                           <  2:954191/100‑1 (MQ=255)
                                                          cAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGt                                           >  1:1188921/1‑100 (MQ=255)
                                                           aaaGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATggggat                                          <  2:1759670/100‑3 (MQ=255)
                                                           aaaGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGtt                                          <  2:670815/100‑1 (MQ=255)
                                                           aaaGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGtt                                          <  1:1535007/100‑1 (MQ=255)
                                                                gAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAgtg                                     <  2:1595498/100‑1 (MQ=255)
                                                                 aGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAgtgt                                    <  1:126256/100‑1 (MQ=255)
                                                                 aGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAgtgt                                    >  1:1968158/1‑100 (MQ=255)
                                                                 aGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAgtgt                                    <  2:1129778/100‑1 (MQ=255)
                                                                  gtaaGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAgtgt                                    <  1:2101025/96‑1 (MQ=255)
                                                                   tCAGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAg                                  <  2:54137/100‑1 (MQ=255)
                                                                     aGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTg                                <  2:2593651/100‑1 (MQ=255)
                                                                     aGTAACGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTg                                >  2:2536107/1‑100 (MQ=255)
                                                                     aGTAACGGATTCGTCATTCCCGCCGTCCTGATTAAAGACGTTCAGGCACCGCTCAACACGCGCTTCACATAATCTGTGTAGCATGCGGCTAAGTGTAGtc                                <  1:1546058/100‑2 (MQ=255)
                                                                      gTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGa                               >  2:611145/1‑100 (MQ=255)
                                                                        aaCGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAAt                             >  1:142795/1‑100 (MQ=255)
                                                                         aCGGATTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAAtt                            <  2:1633324/100‑1 (MQ=255)
                                                                             aTTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCaa                        <  1:1665939/100‑1 (MQ=255)
                                                                             aTTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCaa                        <  2:2472200/100‑1 (MQ=255)
                                                                             aTTCGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCaa                        <  1:1296096/100‑1 (MQ=255)
                                                                              ttCGTCATTCGCGCAGTACTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAAt                       <  1:2394662/100‑1 (MQ=255)
                                                                               tCGTTATGCGCGCAGTCTTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGTTTTACATAATGTGTGTGGCATGGGGTTAAGTGTAGTGAATTTCAAt                       <  1:1932483/99‑1 (MQ=255)
                                                                               tCGTCATTGGCGCAGTCCTGATTACAGAGGGTAAGTAAGCGTTCGACAGGCGGTTTACATTATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCaa                        <  1:130073/98‑1 (MQ=255)
                                                                               tCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATg                      <  1:2249426/100‑1 (MQ=255)
                                                                               tCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATg                      <  2:326243/100‑1 (MQ=255)
                                                                               tCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATg                      <  1:2542884/100‑1 (MQ=255)
                                                                                cGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATga                     >  2:630031/1‑100 (MQ=255)
                                                                                cGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATga                     >  2:2152980/1‑100 (MQ=255)
                                                                                cGTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATga                     >  2:195378/1‑100 (MQ=255)
                                                                                 gTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATgag                    >  2:75061/1‑100 (MQ=255)
                                                                                 gTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATgag                    >  2:247221/1‑100 (MQ=255)
                                                                                 gTCATTCGCGCAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATgag                    <  1:2467997/100‑1 (MQ=255)
                                                                                   gatctgcgcAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGTTTAAGTGTAGTGAATTTCAATGAGaa                  <  2:2281848/95‑1 (MQ=255)
                                                                                   cATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGaa                  >  1:1059099/1‑100 (MQ=255)
                                                                                        gcgcAGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGTTTAAGTGTAGTGAATTTCAATGAGaaagatc             >  1:2281833/1‑96 (MQ=255)
                                                                                            aGTCCTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTACGc         <  2:1537558/100‑1 (MQ=255)
                                                                                               ccTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTACGCttt      >  1:1393493/1‑100 (MQ=255)
                                                                                               ccTGATTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTACGCttt      >  2:2221173/1‑100 (MQ=255)
                                                                                                   gTTAAAGAGGTTAAGTAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATTAGAAACGTTACGCTTTcgcg  <  1:2273980/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAGTGTAGTGAATTTCAATGAGAAACGTTACGCTTTCGCG  >  NC_000913/3645097‑3645295

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