Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,225,841 A→G T672A (ACC→GCC)  metH → cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,225,8410AG98.1% 382.7 / NA 106T672A (ACC→GCC) metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (43/61);  total (45/61)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

CACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCAC  >  NC_000913/4225708‑4225976
                                                                                                                                     |                                                                                                                                       
cACCGAGCGTTTACTCGAGCTTCCCAAGACATATCGCCGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGt                                                                                                                                   <  1:525270/140‑1 (MQ=255)
   cGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTa                                                                                                                                <  1:1835804/140‑1 (MQ=255)
   agAGCGTTTACTGGAGCTTGCCTAGGAATGTCGCGGCTGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTa                                                                                                                                <  1:849422/139‑1 (MQ=255)
        gTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAg                                                                                                                           <  1:1840724/140‑1 (MQ=255)
        gTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAg                                                                                                                           <  1:1850387/140‑1 (MQ=255)
        gTTTACTGGAGCTTGCCGAGAAATATCGCCGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAg                                                                                                                           <  1:898651/140‑1 (MQ=255)
          ttACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCa                                                                                                                         <  2:3748/140‑1 (MQ=255)
           tgcgggAGCTAGCCGTGAAATATCGCGGCAGCAAAACCGACGACGCCGCCAACGCCCAGCAGGCGGAGTGGCGGTCGTGGGAGGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAg                                                                                                                        <  2:1377989/136‑1 (MQ=255)
           tACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAg                                                                                                                        <  1:1585213/140‑1 (MQ=255)
               ggAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATa                                                                                                                    <  1:359609/140‑1 (MQ=255)
               ggAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATa                                                                                                                    <  2:144316/140‑1 (MQ=255)
               ggAGCTTGCCGAGAAATATCGCCGCAGCAAAACCGACTACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATa                                                                                                                    <  1:826607/140‑1 (MQ=255)
                gAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGGTGGTCAAGGGCATTACCGAGTTTATCGAGCAGGATAc                                                                                                                   >  1:528525/1‑140 (MQ=255)
                gAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCAGGAATACTCGCTCGGCAAAGGCATTACCGACCTTATCGAGCAGGAtct                                                                                                                   >  2:849422/1‑138 (MQ=255)
                 aGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTa                                                                                                                                >  1:1568124/1‑126 (MQ=255)
                 aGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTa                                                                                                                                <  2:1568124/126‑1 (MQ=255)
                 aGCTTGCCGAGAAATATCGCGGCAACAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATAcc                                                                                                                  >  2:1798475/1‑140 (MQ=255)
                   cTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCa                                                                                                                         >  1:12480/1‑131 (MQ=255)
                   cTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCa                                                                                                                         <  2:12480/131‑1 (MQ=255)
                   cTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACAga                                                                                                                <  2:1904154/140‑1 (MQ=255)
                     tGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAg                                                                                                                           <  1:1877892/127‑1 (MQ=255)
                     tGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAg                                                                                                                           >  2:1877892/1‑127 (MQ=255)
                       ccGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCgaagaa                                                                                                            >  1:1882035/1‑140 (MQ=255)
                        gccgaAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGaagaag                                                                                                           <  1:109486/137‑1 (MQ=255)
                          agaAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCGGAGTTTATCGAGCAGGATACCGAAGAAGcc                                                                                                         <  2:1883215/140‑1 (MQ=255)
                          agaAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGcc                                                                                                         >  2:1364922/1‑140 (MQ=255)
                          agaAATATCGCGGCAGAGAAACTGACGACGCTGCCTACGTCCAGCAGGAGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGcc                                                                                                         <  2:528525/140‑1 (MQ=255)
                             aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGATTACTCGCTGGTCAAAGGCATTGCCGAGGTTAACGGGCAGGAAACCGAAGAAGCCCGc                                                                                                      >  1:661565/1‑140 (MQ=255)
                                  cgcgGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCcagca                                                                                                 <  1:1864185/140‑1 (MQ=255)
                                  cgcgGCAACAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCcagca                                                                                                 <  1:1798475/140‑1 (MQ=255)
                                   gcgGCAGCAAAACCGACAACACAGCCAACGCTCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCcagcag                                                                                                <  2:784237/140‑1 (MQ=255)
                                      gcagcaAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGGAGGCTACCGAAGATGCCCGCCAGCAGGCt                                                                                             >  1:784237/1‑140 (MQ=255)
                                      gcagcaAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCt                                                                                             >  1:1576540/1‑140 (MQ=255)
                                      gcagcaAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCt                                                                                             >  1:1311648/1‑140 (MQ=255)
                                      gcagcaAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCt                                                                                             >  2:68621/1‑140 (MQ=255)
                                       cagcaAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGccc                                                                                                        <  1:1172817/128‑1 (MQ=255)
                                         gcaAAACCGACGACACCACCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTAcg                                                                                          >  2:1616406/1‑140 (MQ=255)
                                           aaaaCCGACTACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTAcgcg                                                                                        <  2:1073171/140‑1 (MQ=255)
                                           aaaaCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCGGGATACCGAAGAAGCCCGCCAGCAGGCTAcgcg                                                                                        <  2:58328/140‑1 (MQ=255)
                                           aaaaCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTAcgcg                                                                                        <  2:1576540/140‑1 (MQ=255)
                                           aaaaCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTAcgcg                                                                                        <  1:1364922/140‑1 (MQ=255)
                                           aaaaCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTAcgcg                                                                                        >  2:1188609/1‑140 (MQ=255)
                                           aaaaCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATGGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTAcgcg                                                                                        <  1:1616406/140‑1 (MQ=255)
                                             aaCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGcc                                                                                      <  1:1881340/140‑1 (MQ=255)
                                                cgacgaCACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGa                                                                                   <  2:1670358/140‑1 (MQ=255)
                                                cgacgaCACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGa                                                                                   <  1:1791001/140‑1 (MQ=255)
                                                   cgacACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATAGCCAAGAAGCCCGGCAGCAAGCTACGCACCCGATTg                                                                                >  2:923291/1‑140 (MQ=255)
                                                   cgacACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTg                                                                                <  1:68621/140‑1 (MQ=255)
                                                     acacCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGaa                                                                              <  2:1618833/140‑1 (MQ=255)
                                                         cGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGa                                                                          >  2:1515689/1‑140 (MQ=255)
                                                         cGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGa                                                                          >  2:10218/1‑140 (MQ=255)
                                                            cAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTg                                                                       <  1:1088228/140‑1 (MQ=255)
                                                             aaCGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGa                                                                      <  2:1311648/140‑1 (MQ=255)
                                                               cGCCCAGCAGGCGGATTGGCGCACGTGGGCAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAg                                                                    <  1:571584/140‑1 (MQ=255)
                                                                 cccAGCAGGCGGAGTGGTGCTCATGGGAAGTCAATAAACTTCTGGAATTCTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGc                                                                  <  2:763821/140‑1 (MQ=255)
                                                                  ccAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGcc                                                                 <  1:1188609/140‑1 (MQ=255)
                                                                     gcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGt                                                              >  2:1881495/1‑140 (MQ=255)
                                                                      cagGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGtt                                                             >  1:7478/1‑140 (MQ=255)
                                                                        ggCGGAGTGGCGCTCATGGGAAGTGTATACACGTCTGTAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGa                                                           <  1:592628/140‑1 (MQ=255)
                                                                          cGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCGCGCCAGCAGGCTGCGCACCCGATTGAAGTGATTGAAGGCCCGTTGATg                                                         >  2:571584/1‑140 (MQ=255)
                                                                           ggAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGa                                                                      <  1:365939/126‑1 (MQ=255)
                                                                           ggAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGa                                                                      >  2:365939/1‑126 (MQ=255)
                                                                            gAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGa                                                       <  2:1029818/140‑1 (MQ=255)
                                                                            gAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGa                                                       <  2:1825106/140‑1 (MQ=255)
                                                                            gAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGa                                                       <  1:1383910/140‑1 (MQ=255)
                                                                            gAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGa                                                       <  2:7478/140‑1 (MQ=255)
                                                                                ggCGCTCGTGGGAAGCGAATAAGCGACTGCAATACTCACTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGc                                                   <  1:793599/140‑1 (MQ=255)
                                                                                      cGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAAt                                             <  1:1758459/140‑1 (MQ=255)
                                                                                              gTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCgg                                     <  1:1881495/140‑1 (MQ=255)
                                                                                               tGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGAGTGAAGGCGCGTAGATGGACGGCATGAATGTGGTCGGc                                    >  1:508652/1‑140 (MQ=255)
                                                                                                 aaTAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGa                                  >  1:1443516/1‑140 (MQ=255)
                                                                                                 aaTAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGa                                  >  2:1649210/1‑140 (MQ=255)
                                                                                                 aaTAAACGTCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGa                                  >  1:152065/1‑140 (MQ=255)
                                                                                                 aaTAAACGTCTGGAATACTCGCCGGTCAAAGGCATGGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGa                                  <  2:784105/140‑1 (MQ=255)
                                                                                                  aTAAACGTCCGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATCCCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGAc                                 <  2:1853439/140‑1 (MQ=255)
                                                                                                         tCTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTg                          >  2:124552/1‑140 (MQ=255)
                                                                                                          cTGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTgg                         >  2:1851880/1‑140 (MQ=255)
                                                                                                           tGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGc                        >  1:759118/1‑140 (MQ=255)
                                                                                                           tGGAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGc                        >  2:1734992/1‑140 (MQ=255)
                                                                                                             gAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGa                      <  2:1443516/140‑1 (MQ=255)
                                                                                                             gAATACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGa                      <  1:124157/140‑1 (MQ=255)
                                                                                                              aaTACTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGaa                     <  1:1851880/140‑1 (MQ=255)
                                                                                                                  cTCGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGa                 >  2:1333057/1‑140 (MQ=255)
                                                                                                                   tCGCTGGTCAAAGGCATTGCCGAGTCTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGaa                >  2:1287007/1‑140 (MQ=255)
                                                                                                                    cGCTGGTCAAAGGCATTGCCGATTTTATCTAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGaaa               <  1:923291/140‑1 (MQ=255)
                                                                                                                    cGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGaaa               <  2:152065/140‑1 (MQ=255)
                                                                                                                    cGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGACATGAATGTGGTCGGCGACCTGTTTGGCg                       >  1:1221167/1‑132 (MQ=255)
                                                                                                                    cGCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGACATGAATGTGGTCGGCGACCTGTTTGGCg                       <  2:1221167/132‑1 (MQ=255)
                                                                                                                     gCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGaaaa              >  2:1917812/1‑140 (MQ=255)
                                                                                                                     gCTGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGAATGAATGTGGTGGGAGACCTGTTTGGCGAAGGGaaaa              >  1:488405/1‑140 (MQ=255)
                                                                                                                       tGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCg                       >  2:1111342/1‑129 (MQ=255)
                                                                                                                       tGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCg                       <  1:1111342/129‑1 (MQ=255)
                                                                                                                       tGGTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATg            >  2:128745/1‑140 (MQ=255)
                                                                                                                       tGGTCAAAGGCATTCCCGAGTTTATCGAGCAGGATCCCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATg            <  1:1484114/140‑1 (MQ=255)
                                                                                                                        ggTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCTTGAATATGGTCGGCGAACTGTAGGGCGAAGAGAAAATGt           >  2:920704/1‑140 (MQ=255)
                                                                                                                        ggTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGt           >  1:1078641/1‑140 (MQ=255)
                                                                                                                         gtaaaaGGAATTGCCAAGTTTACGGGGCACGATCCCGCAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGAGGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGtt          <  1:652443/137‑1 (MQ=255)
                                                                                                                         gTCAAAGGCATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGtt          <  1:1649210/140‑1 (MQ=255)
                                                                                                                          tCAAAGGCTTTGCCGAGTTTACCGAGCAGTATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTc         <  1:920704/140‑1 (MQ=255)
                                                                                                                                 cATTGCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCcac  >  2:1254721/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCCCGATTGAAGTGATTGAAGGCCCGTTGATGGACGGCATGAATGTGGTCGGCGACCTGTTTGGCGAAGGGAAAATGTTCCTGCCAC  >  NC_000913/4225708‑4225976

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