Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,331,082 C→A 5.3% A12E (GCA→GAA)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,331,0820CA 5.3% 294.0 / 2.7 89A12E (GCA→GAA) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base C (42/41);  new base A (2/3);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAA  >  NC_000913/1330934‑1331228
                                                                                                                                                    |                                                                                                                                                  
aCGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGc                                                                                                                                                    <  2:204437/149‑1 (MQ=255)
  gACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCa                                                                                                                                                   >  2:192077/1‑148 (MQ=255)
   aCAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCaa                                                                                                                                                  >  5:323280/1‑148 (MQ=255)
    cAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTTCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCaaaa                                                                                                                                                <  1:308621/149‑1 (MQ=255)
    cAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCaaaa                                                                                                                                                <  8:173219/149‑1 (MQ=255)
          gTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCaa                                                                                                                                           >  4:18720/1‑148 (MQ=255)
              ttatCCGGAGAGAGCGAGTCCATATCGTTAACTCTTTGACAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCGAAAACg                                                                                                                                       <  4:246404/147‑1 (MQ=255)
               tatCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGAt                                                                                                                                     >  7:209901/1‑149 (MQ=255)
                atCCGCAGAGAGCGAGTCCATATCGGTAAATCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATAAGGGTAAAGCTCTTGTCATCGATGAGTCCCCGGCAAAAGCCAAAACGAt                                                                                                                                     >  2:80561/1‑148 (MQ=255)
                 tCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATc                                                                                                                                    >  6:142479/1‑148 (MQ=255)
                 tCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATc                                                                                                                                    <  8:253019/148‑1 (MQ=255)
                 tCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATc                                                                                                                                    <  8:246658/148‑1 (MQ=255)
                          agCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATc                                                                                                                          <  4:159416/149‑1 (MQ=255)
                           gCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCt                                                                                                                         >  7:183559/1‑149 (MQ=255)
                            cGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATc                                                                                                                                    >  5:200511/1‑137 (MQ=255)
                            cGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATc                                                                                                                                    <  6:200511/137‑1 (MQ=255)
                                    tatCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAGAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCGAAAACGATCAACAAGTATCTGGGTAGTGa                                                                                                                >  6:233097/1‑149 (MQ=255)
                                    tatCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGa                                                                                                                <  1:192077/149‑1 (MQ=255)
                                    tatCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGa                                                                                                                <  8:87697/149‑1 (MQ=255)
                                             cTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtggt                                                                                                       <  8:209901/149‑1 (MQ=255)
                                             cTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtggt                                                                                                       <  8:183559/149‑1 (MQ=255)
                                              tCGTTGCCAGTGGAAGGTTTATCAACGTGCGCCGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCACAACGATCAACAAGTATCTGGGTAGTGACTACgtggt                                                                                                       >  2:108394/1‑148 (MQ=255)
                                              tCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtggt                                                                                                       >  5:408926/1‑148 (MQ=255)
                                                   gCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTGGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAAt                                                                                                  <  4:64630/148‑1 (MQ=255)
                                                    ccAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATc                                                                                                 <  6:408926/148‑1 (MQ=255)
                                                    ccAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGAAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATc                                                                                                 >  2:201748/1‑148 (MQ=255)
                                                    ccAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGAAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATc                                                                                                 >  2:207740/1‑148 (MQ=255)
                                                     cAGTGGAAGGTTTATCAACGTTCGACGGATTCCTGGAAGAATCAAATTATGTAAGGTTAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATcc                                                                                                <  4:369952/148‑1 (MQ=255)
                                                     cAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATcc                                                                                                >  3:369952/1‑148 (MQ=255)
                                                     cAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATc                                                                                                 >  6:235672/1‑147 (MQ=255)
                                                      aGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAg                                                                                              >  6:60949/1‑149 (MQ=255)
                                                       gTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGCAACAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGCGTCCCCGGCAAAAGCCAAAACGAt                                                                                                                                     >  1:297079/1‑109 (MQ=255)
                                                        tGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCg                                                                                            <  6:323280/149‑1 (MQ=255)
                                                        tGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCg                                                                                            >  8:1367/1‑149 (MQ=255)
                                                          gAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTggg                                                                                                                      <  7:374526/120‑1 (MQ=255)
                                                          gAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTggg                                                                                                                      >  8:374526/1‑120 (MQ=255)
                                                            aGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCgg                                                                                        >  3:367617/1‑149 (MQ=255)
                                                                ttATCAACGTGCGACGCATTCCTGGAAGAATAAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGGTGAGTCCCCGGCAAAAGCCCAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTca                                                                                     >  5:33451/1‑148 (MQ=255)
                                                                     aaCGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGTGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCcg                                                                               <  6:258784/149‑1 (MQ=255)
                                                                     aaCGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTAATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATAcg                                                                               >  5:119005/1‑149 (MQ=255)
                                                                     aaCGTCCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGAGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCcg                                                                               >  6:151639/1‑149 (MQ=255)
                                                                      aCGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCcgc                                                                              <  8:193622/149‑1 (MQ=255)
                                                                        gTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGAAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGa                                                                            <  1:201748/149‑1 (MQ=255)
                                                                        gTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGAAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGa                                                                            <  1:207740/149‑1 (MQ=255)
                                                                         tGCGACGCATTCCTGGTAGAATCAAATTAGGTAAGGTGAATATTGGTAAAGCTCTTGGCTTCGTTGAGTCCCCGGGAAATGCCAAAACGATCAACAAGTATCTGGGTAGTGACTAGGTGGTGAAATCCAGCGTCGGTCACATCCGCGAt                                                                           <  5:219969/149‑1 (MQ=255)
                                                                         tGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTAt                                                                                                                           <  6:387822/101‑1 (MQ=255)
                                                                         tGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTAt                                                                                                                           >  5:387822/1‑101 (MQ=255)
                                                                           cGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGAtt                                                                          <  8:79117/148‑1 (MQ=255)
                                                                            gACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGt                                                                                       <  5:298565/134‑1 (MQ=255)
                                                                            gACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGGCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGt                                                                                       >  6:298565/1‑134 (MQ=255)
                                                                             aCGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGc                                                                       <  4:367617/149‑1 (MQ=255)
                                                                                  ttCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGa                                                                                                                >  4:33647/1‑103 (MQ=255)
                                                                                  ttCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGa                                                                                                                <  3:33647/103‑1 (MQ=255)
                                                                                    ccTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCa                                                                 <  5:235672/148‑1 (MQ=255)
                                                                                     cTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAg                                                                <  1:80561/148‑1 (MQ=255)
                                                                                     cTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGt                                                               >  7:246256/1‑149 (MQ=255)
                                                                                        gaagaaTCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGc                                                            >  3:385698/1‑149 (MQ=255)
                                                                                            aaTCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCa                                                                                               <  2:410371/110‑1 (MQ=255)
                                                                                            aaTCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCa                                                                                               >  1:410371/1‑110 (MQ=255)
                                                                                              tCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGt                                                               >  2:427903/1‑140 (MQ=255)
                                                                                              tCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGt                                                               <  1:427903/140‑1 (MQ=255)
                                                                                              tCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                                      >  3:228422/1‑149 (MQ=255)
                                                                                                   ttAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCa                                                  >  3:322078/1‑148 (MQ=255)
                                                                                                      ggTAAGGTGAATATGGGTAAAGCTCTTGTCATCTTTGAGTCCCCGGCAAAAGCCAAAAGGATCATCCAGTATCTGGGTAGTGACTACGTGGTGAAATCGAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaaa                                               <  4:257265/148‑1 (MQ=255)
                                                                                                       gTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGa                                                                                                                >  6:108926/1‑82 (MQ=255)
                                                                                                       gTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGa                                                                                                                <  5:108926/82‑1 (MQ=255)
                                                                                                       gTAAGGCGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAag                                             <  5:390557/149‑1 (MQ=255)
                                                                                                        tAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGAGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAaga                                            <  3:18720/149‑1 (MQ=255)
                                                                                                                    gggTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                                >  8:63347/1‑149 (MQ=255)
                                                                                                                    gggTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTa                                >  5:271261/1‑149 (MQ=255)
                                                                                                                     ggTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAc                               <  7:1367/149‑1 (MQ=255)
                                                                                                                       tAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTATTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTTCCAAAAAGAGTTCCGACTCTAACt                             >  7:131169/1‑149 (MQ=255)
                                                                                                                               ttGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccacc                        >  4:351369/1‑146 (MQ=255)
                                                                                                                                 gTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTggg                                                                                                                      >  6:66345/1‑50 (MQ=255)
                                                                                                                                 gTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTggg                                                                                                                      <  5:66345/50‑1 (MQ=255)
                                                                                                                                 gTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCaa                      <  4:322078/146‑1 (MQ=255)
                                                                                                                                 ctCATCTTTGCGTCCCCGGAATATGCAAAAAGTATCAACAAGTATTTGTGTAGTGACACCGTGGGTAAATCAAGCGTCGGCCACATCCGCTATTTGCCGACCTGTGGCGCATCGGCCAAAAAGAGGGCCGACTCTACCTCCACCAAGa                    <  8:86634/147‑1 (MQ=255)
                                                                                                                                  tCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGc                                                       >  3:20835/1‑112 (MQ=255)
                                                                                                                                  tCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGc                                                       <  4:20835/112‑1 (MQ=255)
                                                                                                                                   cATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACgg                 >  6:235898/1‑149 (MQ=255)
                                                                                                                                    aTCGTTTAGTCCCCGGCAAAAGCCAAAACTATCAACATGTATCTGGTTAGTTACTACGTGTTGAAATCCAGCGTCGGTCAAATCCGCGATTTGCCTACCAGTGGCTCAGCTGCCAAAATGAGTGCCGACTCTACCTCCACCAAGACGGc                <  6:119005/149‑1 (MQ=255)
                                                                                                                                     tCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc                           >  2:144213/1‑137 (MQ=255)
                                                                                                                                     tCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTcc                           <  1:144213/137‑1 (MQ=255)
                                                                                                                                           aGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTaaaaa          <  5:142479/148‑1 (MQ=255)
                                                                                                                                               cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGa                                                                            <  2:355423/78‑1 (MQ=255)
                                                                                                                                               cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGa                                                                            >  1:355423/1‑78 (MQ=255)
                                                                                                                                               cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                                           >  2:147910/1‑111 (MQ=255)
                                                                                                                                               cccGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAagag                                           <  1:147910/111‑1 (MQ=255)
                                                                                                                                                 cGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCaaa    >  3:85045/1‑148 (MQ=255)
                                                                                                                                                  ggCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTaaaa  <  5:60949/149‑1 (MQ=255)
                                                                                                                                                    cAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTAc                               >  2:83372/1‑118 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAA  >  NC_000913/1330934‑1331228

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