Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,522,187 C→T A189T (GCG→ACG)  hofM ← DNA utilization protein HofM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,522,1870CT94.1% 178.7 / ‑5.9 68A189T (GCG→ACG) hofMDNA utilization protein HofM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/39);  minor base G (3/1);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATAT  >  NC_000913/3521946‑3522420
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
gggCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCa                                                                                                                                                                                                                                   <  1:178660/251‑1 (MQ=255)
         attaatgcgCGGTTTCACCAGTACGCCTTCCCAAGCGCCAGCCCCGGCGCGATGGCAAAGTCGCCGCCTGGGGTCGGTAGCGGCGGCTGAGGAACAGTAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATGTCGCGACGCTTGCGGGATCAACGGAAATCGCGGCCGCCAGCGCCTACGCGCTAGTCATCCCCCCCGCCAGTGTGCGCCCCCAGGTAGAGCTCGTCGTCCTCAGCCACTGTTCGtt                                                                                                                                                                                                                         <  2:412741/243‑1 (MQ=255)
               tGGCGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCAc                                                                                                                                                                                                                    <  1:151830/250‑1 (MQ=255)
                  cGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGCCCCCCCCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACgcca                                                                                                                                                                                                                <  1:147558/251‑1 (MQ=255)
                     gTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCCCCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTAGAGCGCGTCTTCCACAGCCACTGTTCGTTATCACgccaggc                                                                                                                                                                                                             <  2:254990/251‑1 (MQ=255)
                     gTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACCGCCTCCCAGGGATCAAATCCACCTTCGCCACCTATCGCGCCGCCTTCGGGATCAACGGATCACGCTGCCGCCCGCTCCTTCGCGCTAGTCCTCCCCCCCGCCAGTTTGCGCCCCCAGCTCTAGCGCGTCGTCCACAGCCACTGTTCGTTATCCCgccaggc                                                                                                                                                                                                             >  2:188942/1‑251 (MQ=255)
                       tcatcaGTCCGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACgccaggcc                                                                                                                                                                                                            <  1:120399/250‑1 (MQ=255)
                                cGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGATGAAGAACAGAAACGGCCTCACAGGGATCAAATCCGCCTTCCCCACAAATCACGACGCTTTCGGGATCAACGGATAACGATGCCGCCCGCTCCTTCGCGCTAGTCATCCCCACCGCCAAATTGCGCCCCCAGCTATAGAGCGTCGTCCACAACCACCGTTCCGTCACACCCCAGGCCAGaactttct                                                                                                                                                                                                  >  1:19014/1‑244 (MQ=255)
                                 gCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCCGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGGTTGCGCCCCCAGCTATAGCGCGTCGGCCACAGCCACTGTTCGTTATCCCTCCAGGCCAGACATTGCtg                                                                                                                                                                                                 >  1:101214/1‑251 (MQ=255)
                                  ccTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAGAAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCCCCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCtg                                                                                                                                                                                                 >  1:197129/1‑250 (MQ=255)
                                         aaGCGCAAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGCCGCACAGAAACGGCCTCCCAGGGATCCAATCCCCCTTCGCCACATATCGCGACGCTTTCGGGATCCACGGACAACGCTGCCGCCAGCTCCTTCGCGCTACTCATCCCCCCCCCCCGTTTGCGCCCCCAGCTACAGCCCCTCGTTCACAGCCACCGGTCGTTATCACccccggcccgaaaatgctggtgcggcg                                                                                                                                                                                         >  2:193974/1‑232 (MQ=255)
                                             gcCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGAGACGCTTTCGGGATCCACGGATAACGCTGCCGCCCGCTCCTTCGTGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATCGCGCGTCGTCCCCAGCCACCGTTCGTTATCAAGCCAGGCCAGACCTTGCTGATGAGAAGGTa                                                                                                                                                                                      >  1:53678/1‑250 (MQ=255)
                                                     cAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATACCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCa                                                                                                                                                                             <  1:277165/251‑1 (MQ=255)
                                                        cgcgATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGCCGAACAGAAACGGCCTCCCAGGGAGCAAATCAGCCTTCGCCACATATCGCGCCGCTTTCGGGGTCAACGGATAAAGCTGCCGCCCGCTCCTTCGCGCCAGTCATCACCACCCCCCGTTTGCGCCCCCACCTATAGCGCGTCGTCCCCAACCACCTTTCGTTATCACTCCAGGCCCGACATTGCTGGTGAGAAGGGAAAAAAGGCAggg                                                                                                                                                                          >  2:76188/1‑250 (MQ=255)
                                                          cgATGGCAAATTCTCCACCGGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAAt                                                                                                                                                                        <  1:251593/251‑1 (MQ=255)
                                                           gATGGCAAAGTCTCCACCGGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCACCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATc                                                                                                                                                                       <  1:188942/251‑1 (MQ=255)
                                                              ggCAAAGTCACCCCCTGGCGTCGGTAGCGGCGGCTGACGAACATAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCAGCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGc                                                                                                                                                                     <  1:363510/250‑1 (MQ=255)
                                                                cAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAAAAGAAACGGCCTCCCCGGGATCAAAACCGCCTTCGCCACATATCGCGACGCTTTCGGGGTCACCGTATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCAACGCCAGTTTGAGCCCCCAGCTATAGCCCGTCGGCCACAGCCCCCGTTCGTTTTCACCCCCGGCCAGCCATTGCTGATGAGAAGGTaaaaaaagaaagaacagctgt                                                                                                                                                                  >  2:315547/1‑236 (MQ=255)
                                                                cAAAGTCCCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCGCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGt                                                                                                                                                                  <  1:182810/251‑1 (MQ=255)
                                                                   aGTCTCCACCGGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGAACAAATCCGCCTCCGCCACATACCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAAt                                                                                                                                                               <  2:46766/251‑1 (MQ=255)
                                                                      ctcCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGGTATCACGCCAGGCCAGACATTGCTGATGTGAAGGTAAAAAAGGCAGGAATCGCTGTAATGcc                                                                                                                                                            >  1:4663/1‑250 (MQ=255)
                                                                               gcggctGTAGCGGCGGCTGACGAACAAAACCGCCCTCCCGGGGATCAAACCCGCCGTCGCCACATATCGCGCCGCTTTCGGGACCAACGGAGAACGCTGCCGCCAGCTCCTTCGCGCTAGTCACCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATcc                                                                                                                                                   <  2:293292/251‑1 (MQ=255)
                                                                               gcggcgGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATACCGCGACGCTTGCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCGCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATcc                                                                                                                                                   <  1:23128/251‑1 (MQ=255)
                                                                                cggcggTAGCGGCGGCTGACGAACAGAAACGGCCCCCCCGGGATCAAATCCGCCTTCGCCCCATATCTCGACGCTTTCGGGATCAAAGGAAAAAGCTGCCGCCAGCTCCCTCGCGCTAGGCATCCCCACCGCCAGTTTTCGCCCCCCACTATAGAGCGTCGGCCACAGCCACTGTTCGTTTTCCCGCCAGGACAGACATTGCTGATGGGAAGGGACAAAAAGCAGGAAACGCCGTAATGCACTGGCAtacc                                                                                                                                                  >  2:314474/1‑248 (MQ=255)
                                                                                  gcggTGGCGGCGGCTGCCGAACAGAAACGGCCTCCAAGGGATCAAACCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAAGGCTGCCGCCAGCTCCTTCGCGCTAGTCACCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCggg                                                                                                                                                <  2:380764/251‑1 (MQ=255)
                                                                                    ggTAGCGGCGGCTGCCGAACAGAAACGGCCCCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATACCGCGGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGTCAGTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGt                                                                                                                                              <  2:212455/251‑1 (MQ=255)
                                                                                    ggTAGCGGCGGCTGAAGAACAGAAAAAGCCTCCCCGGGAACAAATCCGCCTTCGCCACATATCGAGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAACTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAACCACTGTTCGTTATCACGCCAGGCCAGCCATTGCTGATGAGAAGGTAAAAAAAGCAGGAATAGCTGTAATGCACTGGCATCCggggg                                                                                                                                              >  2:342618/1‑250 (MQ=255)
                                                                                       aGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAAGCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGa                                                                                                                                            <  1:294482/250‑1 (MQ=255)
                                                                                         cggcggCTGACGAACAGAAACCGCCTCCCAGGGCTCACATCAGCCTTCGCCACATCTCGACACGCTTTCGGGATCCACGGCTAACCCTGCCGCCAGAACCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCCGACATTGCTGATGAGCAGGTAACAAAGGCAGGAATCGCTGTAATGCCCCGGCATCCGGGGTGATcg                                                                                                                                         >  2:411067/1‑251 (MQ=255)
                                                                                                   cGAACAGAAACGGCCCCCCAGGGAACAAATCCGCCTTCGCCCCAAATCGCGACGCTTTCGGGATCAACGGATATCGCTGCCGCCAGCTCCTTCGCGCATGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCGCTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATg                                                                                                                               <  2:31052/251‑1 (MQ=255)
                                                                                                      aCAGAAACGGCCTCCCAGGGATCAAACCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAac                                                                                                                            <  1:149177/251‑1 (MQ=255)
                                                                                                         gAAACGTCCTCCCAGGGATCAAATCCGCCGTCGCCACTAACCGCGACGCTTTCGGTATCAACGGAGAACGCTTCCGCCAGCTCCTTCGCGCTAGTCACCCCCACCGCCAGGTTGCGCCCCCAGCTATAGCGCGTCGGCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATTCACTGGCATCCGGGGTGATCGCGCTCACATGAACACg                                                                                                                         <  2:340245/251‑1 (MQ=255)
                                                                                                            aCGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCa                                                                                                                       >  1:260486/1‑250 (MQ=255)
                                                                                                              ggCCTCCCAGGGATCAAATCCGCCTTCGCCACATACCGCGACGCTTTCGGGATCAACGAATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCACCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAAcc                                                                                                                    <  1:205878/251‑1 (MQ=255)
                                                                                                                            caaacccactttcctcacaTACCGCGACGTTTTCGGAACACACGGTTAACTCTGCCGCCAGGTCCTTCGCTCTAAACACCCCCTCCGCCAGTTTGCGCCCCCAGATATAGCTCGTCGTCCACAGCCACTGTTCGGTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCTATGGCATCCGGGGTGATCGCGCTCACATGGACACGCAACCTTTCTGCCagagga                                                                                                      <  2:12970/236‑3 (MQ=255)
                                                                                                                             aaaccctcctgcGCCACATAGCGCGACGCTTTCGTGTTCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTaa                                                                                                     <  1:266329/240‑1 (MQ=255)
                                                                                                                                      ttCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTg                                                                                            >  2:359572/1‑251 (MQ=255)
                                                                                                                                       gtgccctctttcgtgactcTTTCGGGCTCACCGGATCACGCTGCAGCCCTCTCCTTCGCGCTAGGCACCCCAGCCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCGCTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCGGTAATGCACTGGCGTCCGGGGTGATCGCGCTCACAGGAACACGCAACCTTGCTTCCAGAGTAAGCAGCGtt                                                                                             <  2:27700/231‑1 (MQ=255)
                                                                                                                                         gccccaTATCGGGACGCTTGCGGGATCAACGGATAACGCGGCCGCCAGCTCCTTCGCGCTAGTCACCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCGTTCGGCCAGAGTAAGCGGCGGTGCCa                                                                                         <  2:179357/247‑1 (MQ=255)
                                                                                                                                                   gctaCGCTTTCGGGATCAACGGATACCTCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTctttgct                                                                               <  1:359249/248‑1 (MQ=255)
                                                                                                                                                   gcgACGCTTTCGGGAGCAACGGATAACGCGGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTGCTGCCAGAGTAAGCAGCGTTGCCAGCTctttgct                                                                               <  2:81174/251‑1 (MQ=255)
                                                                                                                                                           tgcGGGACCAACGGAAAACGCTGCCGCCAGCTCCTCCGCGCTAGTCACCCCCACCGCCAGGTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGTGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGGTGCCAGCTCTTTGCTTTgcgc                                                                         <  2:197648/247‑1 (MQ=255)
                                                                                                                                                             ggggtccAACGGATAACGCTGCCGCCACCTCCTTCGCTCTAGTCACCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTTTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTTCTCGGCa                                                                     <  2:247491/245‑1 (MQ=255)
                                                                                                                                                             cGGGATCCACGGACACCGCTGCCGCCAGCGCCTTCGCGCTAGTAACCCCCTCCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCa                                                                     <  2:39973/251‑1 (MQ=255)
                                                                                                                                                               ggATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGt                                                                   >  1:111451/1‑251 (MQ=255)
                                                                                                                                                                 aTCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTc                                                                  >  1:176353/1‑250 (MQ=255)
                                                                                                                                                                   cAACGTATAACGCTGCCGCCACCTCCTTCGCGCTAGTTATCCCCGCCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCTGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCGCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACg                                                               <  1:379825/251‑1 (MQ=255)
                                                                                                                                                                   cAACGGAAAAAGCTGCAGCCAACTCCTTCGCGCTCGTCATACCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTGATCACGCCAGGCCAGACATTGCTGATGAGACGGTACAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGACCACGCAACATTTCTGCCATCGTAAGCAGCGTTGCCTGCACTTCGCTTTGCGCGGCTGTCACg                                                               >  2:156856/1‑251 (MQ=255)
                                                                                                                                                                    aaCGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCACCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGt                                                              <  1:51145/251‑1 (MQ=255)
                                                                                                                                                                        ggATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTtat                                                           >  1:187630/1‑251 (MQ=255)
                                                                                                                                                                           tAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAgg                                                        >  1:325341/1‑251 (MQ=255)
                                                                                                                                                                            aaCGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGc                                                       >  1:277460/1‑251 (MQ=255)
                                                                                                                                                                                      ccAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGAGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACtgagg                                             >  1:295284/1‑250 (MQ=255)
                                                                                                                                                                                            cctccgggcTAGTCACCCCCTCCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTTTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCt                                       <  2:364629/244‑1 (MQ=255)
                                                                                                                                                                                            ccTTCGCGCTAGACATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCAGCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCt                                       <  2:197129/251‑1 (MQ=255)
                                                                                                                                                                                               tCGCGCTAGACATCCCCACCGCCAGTTGGCGCCCCCAGATATATCGCGTCGTCCACAGCCACTGTTCGGTATCACGCCAGGCCAGACAGTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTGTGCTTTGCGCGGCAGTCACGTTATAGGCGGTACTGAGTGAGTCTTc                                     <  2:57827/250‑1 (MQ=255)
                                                                                                                                                                                                cgcgCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTACGCAGCGTTGCCGGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACtgcggggagcttcgc                                   >  2:128933/1‑238 (MQ=255)
                                                                                                                                                                                                  cgcTAGTCATACCCACAGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGCCATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCCCATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGAGAGTCTTTGCTg                                 >  2:399062/1‑251 (MQ=255)
                                                                                                                                                                                                   gcTAGTCATCCCCACAGCCAGTTTGCGACCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCCGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCggggcatggtgcgtctttcctgt                                >  2:346918/1‑231 (MQ=255)
                                                                                                                                                                                                      aGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGCGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACtagaacgggagagccgcgggg                             >  2:63303/1‑231 (MQ=255)
                                                                                                                                                                                                      aGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACt                                                 <  1:63303/231‑1 (MQ=255)
                                                                                                                                                                                                         cATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCg                           <  1:141589/250‑1 (MQ=255)
                                                                                                                                                                                                             cccACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACCGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGCAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGCGGGAGTCTTTGCTGGAATcacaacg                      >  2:398837/1‑245 (MQ=255)
                                                                                                                                                                                                                       aGTTTGCGCCCCCAGCTATAGCGCGTCGTCCCCAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGATGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGACATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGAGCGGCAGTCACGTTATAGGCGGGACAGGGGTTGTCTTCGCTGTAATccgagctcaaggagtc             >  1:355519/1‑236 (MQ=255)
                                                                                                                                                                                                                        gtgtgctcccccAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGTCATTGCTGATGAGAAGGTAAAAAAGGCAGGAACCGCTTTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAATCGCAGGTAGTcc            <  2:399808/244‑1 (MQ=255)
                                                                                                                                                                                                                               cccccAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAAGCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCa     <  1:163611/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccccAGCTATAGCGCGTCGTCCACAGCCACTGTTCGGTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCat    <  1:13659/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cAGCTATAGCGCGTCGTCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGGATGGGTGAGGCTTCGCTGGAAACAAAGCGCACGGGGTCCGGGACCatat  >  2:236247/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
GGGCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCAACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCACCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATAT  >  NC_000913/3521946‑3522420

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