Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,432,470 T→G 5.1% Q196P (CAG→CCG) ‡ ytfF ← DMT transporter family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,432,4700TG 5.1% 281.9 / 6.2 78Q196P (CAG→CCG) ‡ytfFDMT transporter family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base G (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TGAGCAAAGCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCGGGGAGTCC  >  NC_000913/4432322‑4432618
                                                                                                                                                    |                                                                                                                                                    
tGAGCAAAGCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCt                                                                                                                                                      >  4:53958/1‑149 (MQ=255)
  aGCAAAGCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTgc                                                                                                                                                    >  3:668264/1‑149 (MQ=255)
    cAAAGCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTgcgc                                                                                                                                                  >  2:427439/1‑149 (MQ=255)
    cAAAGCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTgcgc                                                                                                                                                  <  1:204567/149‑1 (MQ=255)
        gCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGtt                                                                                                                                              <  4:525122/149‑1 (MQ=255)
        gCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCAGGCGCTggt                                                                                                                                              <  1:218406/149‑3 (MQ=255)
          aCGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGc                                                                                                                                            <  1:205791/149‑1 (MQ=255)
           cGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGcc                                                                                                                                           >  4:604404/1‑149 (MQ=255)
           cGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGcc                                                                                                                                           <  3:147103/149‑1 (MQ=255)
             gCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCca                                                                                                                                         >  3:464436/1‑149 (MQ=255)
                 tCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCAcatc                                                                                                                                     >  1:321496/1‑149 (MQ=255)
                  cGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCAcatca                                                                                                                                    >  2:92195/1‑149 (MQ=255)
                   gCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCAcatcat                                                                                                                                   >  1:496209/1‑149 (MQ=255)
                    cAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACatcatc                                                                                                                                  <  1:92195/149‑1 (MQ=255)
                    cAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACatcatc                                                                                                                                  <  4:293947/149‑1 (MQ=255)
                    cAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACatcatc                                                                                                                                  <  2:146363/149‑1 (MQ=255)
                         aTCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATg                                                                                                                             >  2:21736/1‑149 (MQ=255)
                            aGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGttt                                                                                                                          >  4:333981/1‑149 (MQ=255)
                                 aataaACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCgg                                                                                                                     <  1:427439/149‑1 (MQ=255)
                                   taaACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGa                                                                                                                   <  4:161655/149‑1 (MQ=255)
                                   taaACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGa                                                                                                                   <  1:446730/149‑1 (MQ=255)
                                      acacCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTc                                                                                                                       <  1:535712/142‑1 (MQ=255)
                                      acacCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTc                                                                                                                       >  2:535712/1‑142 (MQ=255)
                                         ccAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTcc                                                                                                             <  1:151392/149‑1 (MQ=255)
                                           aGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCg                                                                                                           >  1:218695/1‑149 (MQ=255)
                                                  gggggCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGcagccag                                                                                                    >  2:613103/1‑149 (MQ=255)
                                                     ggCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCgg                                                                                                 >  1:515352/1‑149 (MQ=255)
                                                       ccAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCAGGCGCTggt                                                                                                                                              <  3:680540/102‑3 (MQ=255)
                                                       ccAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCAGGCGCTggt                                                                                                                                              >  4:680540/1‑100 (MQ=255)
                                                        cAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCg                                                                                              >  1:683596/1‑149 (MQ=255)
                                                        cAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCg                                                                                              >  4:494244/1‑149 (MQ=255)
                                                         aaaGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGcagc                                                                                                       <  1:620950/139‑1 (MQ=255)
                                                         aaaGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGcagc                                                                                                       >  2:620950/1‑139 (MQ=255)
                                                         aaaGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGt                                                                                             >  1:68423/1‑149 (MQ=255)
                                                         aaaGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGt                                                                                             >  1:84658/1‑149 (MQ=255)
                                                          aaGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGtt                                                                                            <  1:357059/149‑1 (MQ=255)
                                                          aaGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGtt                                                                                            <  1:83034/149‑1 (MQ=255)
                                                            gggTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTgc                                                                                          <  3:604404/149‑1 (MQ=255)
                                                                 gggAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAgg                                                                                     <  3:53958/149‑1 (MQ=255)
                                                                     gaAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCAt                                                                                 <  3:578369/149‑1 (MQ=255)
                                                                      aaGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATa                                                                                >  2:466867/1‑149 (MQ=255)
                                                                       aGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATAc                                                                               <  3:189321/149‑1 (MQ=255)
                                                                        gTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATAcc                                                                              >  2:252956/1‑149 (MQ=255)
                                                                          ccGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCAc                                                                            <  2:321496/149‑1 (MQ=255)
                                                                           cGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACg                                                                           >  3:524260/1‑149 (MQ=255)
                                                                           cGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACg                                                                           >  3:397236/1‑149 (MQ=255)
                                                                            ggCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGc                                                                          >  3:371716/1‑149 (MQ=255)
                                                                                tttGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCcag                                                                      <  1:252956/149‑1 (MQ=255)
                                                                                   gCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAgg                                                                                     >  2:500255/1‑131 (MQ=255)
                                                                                   gCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAgg                                                                                     <  1:500255/131‑1 (MQ=255)
                                                                                   gCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCcagcag                                                                   >  3:479015/1‑149 (MQ=255)
                                                                                     gTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGAc                                                                 <  2:683596/149‑1 (MQ=255)
                                                                                       aTTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCa                                                               <  1:613103/149‑1 (MQ=255)
                                                                                            gCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGaa                                                          <  2:68423/149‑1 (MQ=255)
                                                                                            gCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGaa                                                          >  3:376159/1‑149 (MQ=255)
                                                                                                 tAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTaa                                                     >  2:169056/1‑149 (MQ=255)
                                                                                                    cAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGc                                                  <  2:496209/149‑1 (MQ=255)
                                                                                                     aaGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGcc                                                 >  3:619515/1‑149 (MQ=255)
                                                                                                                          gTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATaa                            >  3:567572/1‑149 (MQ=255)
                                                                                                                           tGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAATGCCAGCACGATGCCAGAGGt                                <  2:641685/144‑1 (MQ=255)
                                                                                                                           tGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGt                                >  1:641685/1‑144 (MQ=255)
                                                                                                                             aaaCCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGt                         >  3:607947/1‑149 (MQ=255)
                                                                                                                             aaaCCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGt                         >  1:139782/1‑149 (MQ=255)
                                                                                                                               aCCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGc                       >  3:255572/1‑149 (MQ=255)
                                                                                                                                ccGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGcc                      <  4:524260/149‑1 (MQ=255)
                                                                                                                                       tgtgACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGt                         <  2:15971/139‑1 (MQ=255)
                                                                                                                                       tgtgACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGt                         >  1:15971/1‑139 (MQ=255)
                                                                                                                                       tgtgACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCaa               >  1:627886/1‑149 (MQ=255)
                                                                                                                                           aacAGCGCAGGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGATCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATc           >  3:596371/3‑149 (MQ=255)
                                                                                                                                           aCCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCCGTCAAAATc           >  2:395322/1‑149 (MQ=255)
                                                                                                                                           aCCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATc           >  2:219865/1‑149 (MQ=255)
                                                                                                                                             cAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCgg         <  2:515352/149‑1 (MQ=255)
                                                                                                                                             cAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCACAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCgg         >  3:481249/1‑149 (MQ=255)
                                                                                                                                              aGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCggg        <  1:688829/149‑1 (MQ=255)
                                                                                                                                              aGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCggg        <  3:333981/149‑1 (MQ=255)
                                                                                                                                              aGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCggg        <  4:619515/149‑1 (MQ=255)
                                                                                                                                                cgcCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCGGGGa      <  2:218695/149‑1 (MQ=255)
                                                                                                                                                cgcCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCGGGGa      <  4:668264/149‑1 (MQ=255)
                                                                                                                                                  ccTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCGGGGAGt    <  1:169056/149‑1 (MQ=255)
                                                                                                                                                    tGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCGGGGAGTcc  >  4:588126/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGAGCAAAGCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCGGGGAGTCC  >  NC_000913/4432322‑4432618

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑