Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,257 G→C 8.9% L640L (CTG→CTC fhlA → formate hydrogenlyase transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,2570GC 8.9% 377.6 / 22.1 111L640L (CTG→CTCfhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (54/47);  new base C (8/2);  total (62/49)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 116 total aligned reads displayed.

ACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCG  >  NC_000913/2856109‑2856404
                                                                                                                                                    |                                                                                                                                                   
aCATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTg                                                                                                                                                     >  1:761867/1‑149 (MQ=255)
aCATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTg                                                                                                                                                     >  1:386626/1‑149 (MQ=255)
 cATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTgg                                                                                                                                                    >  3:1309867/1‑149 (MQ=255)
           ccGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGa                                                                                                                                          <  1:784255/149‑1 (MQ=255)
            cGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGAt                                                                                                                                         >  3:933292/1‑149 (MQ=255)
            cGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGAt                                                                                                                                         <  2:465844/149‑1 (MQ=255)
             gggTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATg                                                                                                                                        <  3:852814/149‑1 (MQ=255)
             gggTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATg                                                                                                                                        >  4:668726/1‑149 (MQ=255)
               gTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGaa                                                                                                                                      <  4:846511/149‑1 (MQ=255)
               gTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGaa                                                                                                                                      <  3:357524/149‑1 (MQ=255)
               gTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGaa                                                                                                                                      <  3:1384018/149‑1 (MQ=255)
                  aCGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAatat                                                                                                                                   <  2:988548/149‑1 (MQ=255)
                    gTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCa                                                                                                                                 <  4:1063698/149‑1 (MQ=255)
                    gTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCa                                                                                                                                 <  4:804269/149‑1 (MQ=255)
                    gTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCa                                                                                                                                 >  1:1421227/1‑149 (MQ=255)
                     tACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAg                                                                                                                                <  2:1075099/149‑1 (MQ=255)
                       cgcgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATAAATATCAGtt                                                                                                                              <  2:45278/149‑1 (MQ=255)
                        gcgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTg                                                                                                                             >  2:157665/1‑149 (MQ=255)
                              tGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATtgtg                                                                                                                       >  1:1335480/1‑149 (MQ=255)
                              tGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATtgtg                                                                                                                       >  3:491407/1‑149 (MQ=255)
                               ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgc                                                                                                                      >  2:461465/1‑149 (MQ=255)
                                   aaCGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTCCAGGGCGacaac                                                                                                                                         >  1:1220194/1‑122 (MQ=255)
                                   aaCGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTCCAGGGCGacaac                                                                                                                                         <  2:1220194/126‑5 (MQ=255)
                                    aCGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGGTGATTGTGCGCGTg                                                                                                                 >  4:827579/1‑149 (MQ=255)
                                     cGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGc                                                                                                                >  4:435364/1‑149 (MQ=255)
                                     cGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGc                                                                                                                >  3:233924/1‑149 (MQ=255)
                                      gTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCt                                                                                                               <  1:280395/149‑1 (MQ=255)
                                      gTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCt                                                                                                               >  4:1237396/1‑149 (MQ=255)
                                         aTTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaa                                                                                                            <  1:177495/149‑1 (MQ=255)
                                         aTTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaa                                                                                                            <  3:1049868/149‑1 (MQ=255)
                                          ttGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                           <  2:26031/149‑1 (MQ=255)
                                            gAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaaga                                                                                                         <  3:566619/149‑1 (MQ=255)
                                              gcgcgcGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaaa                                                                                                       <  1:1276341/149‑1 (MQ=255)
                                              gcgcgcGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTCATTGTGCGCGTGCTgaaagaaa                                                                                                       <  1:917294/149‑1 (MQ=255)
                                              gcgcgcGGTATTGCTAACACGCGGTAACGGGCTGCAGCTGTCATTGCCAGATATTGTTTTACGGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaaa                                                                                                       <  4:631171/149‑1 (MQ=255)
                                                        ttGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGaa                                                                                                                                      >  4:162784/1‑108 (MQ=255)
                                                         tGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTg                                                                                            <  1:492583/149‑1 (MQ=255)
                                                            tAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGtt                                                                                         >  1:832136/1‑149 (MQ=255)
                                                                acGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCg                                                                                     >  4:538363/1‑149 (MQ=255)
                                                                        aCGTTCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAGCCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTaaa                                                                             <  3:1315008/149‑1 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGTCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAgg                                                                           >  1:897754/1‑149 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAgg                                                                           >  3:187531/1‑149 (MQ=255)
                                                                          gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAgg                                                                           >  2:527259/1‑149 (MQ=255)
                                                                             ctgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGgcgc                                                                        >  3:250364/1‑149 (MQ=255)
                                                                             ctgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGgcgc                                                                        <  1:358247/149‑1 (MQ=255)
                                                                              tgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCt                                                                       >  4:1203300/1‑149 (MQ=255)
                                                                                    tGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCaa                                                                 >  4:401204/1‑149 (MQ=255)
                                                                                    tGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCaa                                                                 >  4:584044/1‑149 (MQ=255)
                                                                                     gTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGa                                                                                                                                                   >  2:164208/1‑66 (MQ=255)
                                                                                        aTTGCCGGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTc                                                             >  2:1374343/1‑149 (MQ=255)
                                                                                        aTTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTc                                                             >  3:335662/1‑149 (MQ=255)
                                                                                            ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAAc                                                                                                  >  3:1201888/1‑108 (MQ=255)
                                                                                            ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAAc                                                                                                  <  4:1201888/108‑1 (MQ=255)
                                                                                            ccAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTggg                                                         <  3:435364/149‑1 (MQ=255)
                                                                                             cAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGa                                                                                                                                                   >  2:883978/1‑58 (MQ=255)
                                                                                             cAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGa                                                                                                                                                   <  1:883978/58‑1 (MQ=255)
                                                                                             cAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgggg                                                        >  1:682634/1‑149 (MQ=255)
                                                                                             cAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgggg                                                        >  4:141463/1‑149 (MQ=255)
                                                                                               gATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCt                                                      <  1:1139897/149‑1 (MQ=255)
                                                                                                  atTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGaa                                                   >  4:403369/1‑149 (MQ=255)
                                                                                                      ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGTGCTGAAACGc                                               >  4:340108/1‑149 (MQ=255)
                                                                                                      ttttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                               >  2:713266/1‑149 (MQ=255)
                                                                                                        ttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCAc                                             <  1:600398/149‑1 (MQ=255)
                                                                                                           ccGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAc                                          >  3:386529/1‑149 (MQ=255)
                                                                                                             ggAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAccc                                        <  3:183621/149‑1 (MQ=255)
                                                                                                              gAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCct                                       <  4:1045121/149‑1 (MQ=255)
                                                                                                              gAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCct                                       >  2:457353/1‑149 (MQ=255)
                                                                                                               aaCCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctg                                      >  3:280704/1‑149 (MQ=255)
                                                                                                                  cTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgctg                                   >  1:1395231/1‑149 (MQ=255)
                                                                                                                  cTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgctg                                   >  4:1370630/1‑149 (MQ=255)
                                                                                                                   tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                                  >  1:163415/1‑149 (MQ=255)
                                                                                                                   tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                                  >  2:1091839/1‑149 (MQ=255)
                                                                                                                   tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                                  >  1:688620/1‑149 (MQ=255)
                                                                                                                   tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                                  >  1:679251/1‑149 (MQ=255)
                                                                                                                    gAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgggg                                                        >  3:273192/1‑126 (MQ=255)
                                                                                                                    gAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTgggg                                                        <  4:273192/126‑1 (MQ=255)
                                                                                                                    gAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTc                                 >  3:102583/1‑149 (MQ=255)
                                                                                                                         gccgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGa                            <  1:1091839/149‑1 (MQ=255)
                                                                                                                              cTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAAc                                                 <  1:1147400/123‑1 (MQ=255)
                                                                                                                              cTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAAc                                                 >  2:1147400/1‑123 (MQ=255)
                                                                                                                                   gCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCt                  >  1:645169/1‑149 (MQ=255)
                                                                                                                                   gCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCt                  >  3:1176433/1‑149 (MQ=255)
                                                                                                                                   gCAACGGTTGTCGCCCTGGAGGGCGAAGAAGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCt                  >  3:1206283/1‑149 (MQ=255)
                                                                                                                                        ggTTGTCGCCCTGGCGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGaa             >  4:438631/1‑149 (MQ=255)
                                                                                                                                        catcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGaa             >  2:487500/6‑149 (MQ=255)
                                                                                                                                         gTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAt            <  1:1053057/149‑1 (MQ=255)
                                                                                                                                          tcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAtt           >  2:1132542/4‑149 (MQ=255)
                                                                                                                                          tcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAtt           >  2:847743/4‑149 (MQ=255)
                                                                                                                                           tGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTg          <  1:116442/149‑1 (MQ=255)
                                                                                                                                           cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTg          >  4:708669/3‑149 (MQ=255)
                                                                                                                                            ttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAtt           <  1:827778/146‑1 (MQ=255)
                                                                                                                                            ttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAtt           >  2:827778/2‑147 (MQ=255)
                                                                                                                                            ttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa         >  2:1143021/2‑149 (MQ=255)
                                                                                                                                            ttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa         >  1:1453947/2‑149 (MQ=255)
                                                                                                                                            gTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa         <  1:174063/149‑1 (MQ=255)
                                                                                                                                                cccTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATaaa     <  2:1488832/149‑1 (MQ=255)
                                                                                                                                                  cTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATc   >  1:1010440/1‑149 (MQ=255)
                                                                                                                                                   tGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCg  <  4:233924/149‑1 (MQ=255)
                                                                                                                                                    ggAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGt                                                                                             <  1:309295/57‑1 (MQ=255)
                                                                                                                                                    ggAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGt                                                                                             >  2:309295/1‑57 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCG  >  NC_000913/2856109‑2856404

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: