Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,258 G→C 5.4% E641Q (GAG→CAG)  fhlA → formate hydrogenlyase transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,2580GC 5.4% 389.9 / 8.3 121E641Q (GAG→CAG) fhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (55/59);  new base C (5/2);  total (60/61)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 123 total aligned reads displayed.

CATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCA  >  NC_000913/2856110‑2856407
                                                                                                                                                    |                                                                                                                                                     
cATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGa                                                                                                                                                      >  2:974468/1‑150 (MQ=255)
cATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAAAGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGAGATTGTTTTACCGGACCCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTggg                                                                                                                                                      >  1:582562/1‑149 (MQ=255)
   ggAGTGGCCGGGTAACGTACTCGAAGTGGAAAACGTCATTGAGCTCGCGGTATTGCTAACTCGCGGTAACGTGCTTCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAggg                                                                                                                                                   <  2:1127675/150‑1 (MQ=255)
   ggAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAggg                                                                                                                                                   <  1:1073010/150‑1 (MQ=255)
           cGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTCCAGGGCGacaacc                                                                                                                                           <  1:1890572/150‑6 (MQ=255)
              gTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAat                                                                                                                                        <  2:554056/150‑1 (MQ=255)
              gTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAat                                                                                                                                        >  2:13531/1‑150 (MQ=255)
               tAACGTACGCGCACTGCAAAACGTCCTTGAGCGCGCGGTATAGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAata                                                                                                                                       <  2:150842/150‑1 (MQ=255)
                aaCGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGATCCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAttat                                                                                                                                      >  2:524289/1‑150 (MQ=255)
                  cGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCa                                                                                                                                    >  2:1677754/1‑150 (MQ=255)
                  cGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCa                                                                                                                                    >  1:1304016/1‑150 (MQ=255)
                   gTACGCGCACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACCCGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTAACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAg                                                                                                                                   <  1:1884557/150‑1 (MQ=255)
                   gTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGTTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAg                                                                                                                                   <  2:1198963/150‑1 (MQ=255)
                    tACGCGAACTGGAACACGTCATTGAGCGCGCGGTATTGCAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGtt                                                                                                                                 <  2:1014189/150‑1 (MQ=255)
                    tACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGt                                                                                                                                  >  1:1721315/1‑150 (MQ=255)
                    tACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGt                                                                                                                                  >  1:1713479/1‑150 (MQ=255)
                    tACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGt                                                                                                                                  >  2:1350687/1‑150 (MQ=255)
                     ccgcgAACTGGAAAACGTCATAGACCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTCCCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGtt                                                                                                                                 <  2:486476/149‑1 (MQ=255)
                     aCGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAg                                                                                                                                              >  2:1016062/1‑137 (MQ=255)
                     aCGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAg                                                                                                                                              <  1:1016062/137‑1 (MQ=255)
                      ctcgcaatGGAAAACGTCATTGAGCGCGCGGTATTGCTAACTCGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTg                                                                                                                                <  2:941380/143‑1 (MQ=255)
                      cgcgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGATTATCAGTTg                                                                                                                                >  2:1167845/1‑150 (MQ=255)
                              ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcg                                                                                                                        >  1:54720/1‑150 (MQ=255)
                              ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcg                                                                                                                        <  1:1591460/150‑1 (MQ=255)
                              ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcg                                                                                                                        >  2:78250/1‑150 (MQ=255)
                               gAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcgc                                                                                                                       >  1:210690/1‑150 (MQ=255)
                                aaaaCGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCGGCCGAAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGcgcg                                                                                                                      <  1:42074/150‑1 (MQ=255)
                                  aaCGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTg                                                                                                                    <  1:1350687/150‑1 (MQ=255)
                                   aCGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGc                                                                                                                   <  2:726527/150‑1 (MQ=255)
                                       cATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGGGGGCGAAGATGAATATCAGTTGATTGTGCGCGGGCTgaa                                                                                                               <  2:1454480/150‑1 (MQ=255)
                                       cATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCACTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaa                                                                                                               >  1:352798/1‑150 (MQ=255)
                                        aTTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaa                                                                                                              <  1:1047468/150‑1 (MQ=255)
                                         ttGAGCGCGCGGTATTGCTAACACGCGGTACCGTGCGGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGGGCTgaaag                                                                                                             <  1:980670/150‑1 (MQ=255)
                                           gAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaaagaa                                                                                                           <  2:770456/150‑1 (MQ=255)
                                              cgcgcgGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCAGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAAcc                                                                                                        <  2:1563718/150‑1 (MQ=255)
                                                         gCTACCACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTGTTACCGGAACCTGAAACGCCGCCTCCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTggg                                                                                             <  2:1116148/150‑2 (MQ=255)
                                                         gCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGCCGTGGt                                                                                             <  1:1015352/150‑1 (MQ=255)
                                                          cTACCACGCGCTAACGTGCTGCCGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCAAGCCGCAACGGTTGTCGACCTGGAGAGAGAAGATGCATATCACTTAATTCTCCGCTTGCTGAAAGAAACCAACGTCGTGGtt                                                                                            >  1:79815/1‑150 (MQ=255)
                                                           tAACCCGCGGTAACGTGCTGCAGCTGTCATTGCCAGATCTTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTg                                                                                           <  1:92130/150‑1 (MQ=255)
                                                           tAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTg                                                                                           >  1:535817/1‑150 (MQ=255)
                                                            aaCACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGc                                                                                          <  2:77921/150‑1 (MQ=255)
                                                            aaCACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACAGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGGTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGc                                                                                          >  2:1106922/1‑150 (MQ=255)
                                                             acacGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGcc                                                                                         <  1:1126414/150‑1 (MQ=255)
                                                              caggcgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACAGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGTTGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCg                                                                                        <  1:947296/147‑1 (MQ=255)
                                                              cacGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCg                                                                                        >  2:146395/1‑150 (MQ=255)
                                                                cgcgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCggg                                                                                      <  1:751249/150‑1 (MQ=255)
                                                                cgcgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCggg                                                                                      >  1:962295/1‑150 (MQ=255)
                                                                  cgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCGCTGGAGGGCGAAGATGAATAGCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGcc                                                                                    <  1:65166/150‑1 (MQ=255)
                                                                  cgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGcc                                                                                    <  2:1583443/150‑1 (MQ=255)
                                                                  cgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGcc                                                                                    <  1:406400/150‑1 (MQ=255)
                                                                  cgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGGCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGcc                                                                                    <  1:1903604/150‑1 (MQ=255)
                                                                  cgGTAACGTGCTGCAGCTGGCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGcc                                                                                    <  2:1215477/150‑1 (MQ=255)
                                                                  cgGTAACGTGCGGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGcc                                                                                    <  2:1605238/150‑1 (MQ=255)
                                                                     tAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTaa                                                                                 >  2:1634343/1‑150 (MQ=255)
                                                                         gTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGgc                                                                             <  1:1026992/150‑1 (MQ=255)
                                                                             tgcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAACCCACCGGCGTGGTTGCCGGGCCTAAAGGCGCTg                                                                         >  1:393039/1‑150 (MQ=255)
                                                                               cAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTgcg                                                                       <  1:476650/150‑1 (MQ=255)
                                                                                aGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTgcgc                                                                      >  2:1721225/1‑150 (MQ=255)
                                                                                aGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTgcgc                                                                      >  1:781927/1‑150 (MQ=255)
                                                                                 gCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCa                                                                     <  1:955173/150‑1 (MQ=255)
                                                                                 gCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCa                                                                     >  1:1576050/1‑150 (MQ=255)
                                                                                 gCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCa                                                                     >  1:1879547/1‑150 (MQ=255)
                                                                                  cTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCaa                                                                    >  2:934874/1‑150 (MQ=255)
                                                                                       aTTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCt                                                               >  1:727263/1‑150 (MQ=255)
                                                                                        ttGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTg                                                              >  1:1880924/1‑150 (MQ=255)
                                                                                          gCCAGATATTGTTTTACCGGAACCTGTAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTggg                                                            <  1:593920/150‑1 (MQ=255)
                                                                                          gCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTggg                                                            <  2:575215/150‑1 (MQ=255)
                                                                                          gCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTggg                                                            <  2:1614801/150‑1 (MQ=255)
                                                                                                 atTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGaaa                                                     >  1:947446/1‑150 (MQ=255)
                                                                                                    gTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                                  >  2:1472086/1‑150 (MQ=255)
                                                                                                      tttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACTGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCAc                                                <  1:1167845/150‑1 (MQ=255)
                                                                                                      tttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGGCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTCGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCAc                                                <  2:1890097/150‑1 (MQ=255)
                                                                                                        tACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACgc                                              >  1:1432104/1‑149 (MQ=255)
                                                                                                         aCCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAc                                             <  2:1688921/150‑1 (MQ=255)
                                                                                                           cGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGCAACGCACGAccc                                           >  1:1867528/1‑150 (MQ=255)
                                                                                                                ccTGAACCGCCGCCTGCCGCAACGGTTGTCGCCCTGGTGGGCGAAGATTAATATCAGTTGATTGTGCGCGTGCGGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgctg                                      <  2:1576050/150‑1 (MQ=255)
                                                                                                                 cTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                                     >  2:1076124/1‑150 (MQ=255)
                                                                                                                  tGAACCGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATAGCAGTTGATTGTGCGCGTGCTGAAAGACACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAAAGTCTGGGGCTGAAACGCACGACCCTGCTGTc                                    <  1:693661/150‑1 (MQ=255)
                                                                                                                  tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGATGAGTGTGCGCGTGCTGAAAGAACCCACCGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCATCGTCTGGGGCTGAAACGCACGACCCTGCTGTc                                    >  1:107379/1‑150 (MQ=255)
                                                                                                                        gccgccTGCCGCAACGGATGTCGCCCTGGAGGCCGAAGATGAATAACAGTTGATTGTGCGCGAGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAATGCGCTGCGCAACGGCTGGGGCTGTAACGCACGACCCTGCTGTCACGGAt                              <  1:173178/150‑1 (MQ=255)
                                                                                                                         ccgccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGAGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATg                             >  2:22935/1‑150 (MQ=255)
                                                                                                                            ccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGCAACGCACGACCCTGCTGTCACGGATGAAg                          >  1:1363675/1‑150 (MQ=255)
                                                                                                                            ccTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAg                          >  2:232167/1‑150 (MQ=255)
                                                                                                                              tGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATTAATATCAGTTGATTGTGCGCGTGCTGAAAGCAACCAACGGCGTGGTTGCCGGGACTAAAGGCGATGCGCAACGTCTGGGGCTGAAACGCACGCCCCTGCTGTCACGGATGAAGCg                        >  1:1341426/1‑150 (MQ=255)
                                                                                                                              tGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCg                        <  2:54720/150‑1 (MQ=255)
                                                                                                                               gCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGACTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCgg                       >  2:259696/1‑150 (MQ=255)
                                                                                                                                 cGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTAGCCGGGCCGAAAGGCGCTGCGCACCGTCTGGGGCTGAAACGCACGACCCTGCAGTCACGGATGAGGCGGCt                     >  2:926198/1‑150 (MQ=255)
                                                                                                                                     ttcatcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGa                 >  1:1483385/8‑150 (MQ=255)
                                                                                                                                      tcatcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGAGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATCAACCGGCTGGGaa                >  2:348052/7‑150 (MQ=255)
                                                                                                                                      ggggtGTCGCCCTGCAGGGCGAAGTTGAATATCAGTTGATTGTGCGCGTGCTTAAAGTACCCAACGGCGGGGGTGCCGGGCCTAAAGGCGCTGCGCAACGGCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGaa                <  2:191452/146‑1 (MQ=255)
                                                                                                                                       ggTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAt               <  2:1528765/150‑1 (MQ=255)
                                                                                                                                        gTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTGGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAtt              <  1:192734/150‑1 (MQ=255)
                                                                                                                                         ttGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTg             <  1:908490/150‑1 (MQ=255)
                                                                                                                                         tcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTg             >  1:704216/4‑150 (MQ=255)
                                                                                                                                          tGTCTCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGAGTGCTGATCGACAGGAACGGCGGGGTGGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGACACGCACGACCCTGCTGTCACGGATGAAGCGGCGGGGAATTGa            <  1:289033/150‑1 (MQ=255)
                                                                                                                                          ggTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa            <  1:1041846/149‑1 (MQ=255)
                                                                                                                                           gTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGAt           >  1:1625161/1‑150 (MQ=255)
                                                                                                                                            tCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTTCTGAAAGAAAGCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATa          <  2:210690/150‑1 (MQ=255)
                                                                                                                                            tCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATa          >  1:642374/1‑150 (MQ=255)
                                                                                                                                                    gAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCa  >  1:615446/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCA  >  NC_000913/2856110‑2856407

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