Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,432,469 C→A 5.1% Q196H (CAG→CAT) ‡ 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,4690CA 5.1% 280.3 / 4.4 78Q196H (CAG→CAT) ‡ytfFDMT transporter family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base A (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

TGAGCAAAGCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCGGGGAGT  >  NC_000913/4432322‑4432616
                                                                                                                                                   |                                                                                                                                                   
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                       >  1:15971/1‑139 (MQ=255)
                                                                                                                                       tgtgACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGt                       <  2:15971/139‑1 (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      <  3:333981/149‑1 (MQ=255)
                                                                                                                                              aGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCggg      <  1:688829/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)
                                                                                                                                                   |                                                                                                                                                   
TGAGCAAAGCACGGCTATCGCAACCATCAGACTAATAAACACCAGCGGACGGGGGCCAAAGGGTAGGGAGAAGTCCGGCGTTTGCGTATTCAGCCAGTAACAAGCGACGAGATAGCCGATAAGTGAAACCGGCAGTGTGACCAGCGCCTGCGCTGTTGCCCACATCATCGGATGTTTGTCGGGATTTTCCCGCAGCCAGCGGGCGTTGCGCAGGGCATACCACGCCCAGCAGACCACGGAAACTAACGCCAGCACGATGCCAGAGGTATAACGTGCCCAGTCAAAATCGGGGAGT  >  NC_000913/4432322‑4432616

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: