Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,151,901 T→C 5.8% T400A (ACG→GCG)  ppc ← phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,151,9010TC 5.8% 210.1 / 3.0 70T400A (ACG→GCG) ppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (37/28);  new base C (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

TTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCGGTTCCCACAGTT  >  NC_000913/4151757‑4152047
                                                                                                                                                |                                                                                                                                                  
ttGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGctc                                                                                                                                                >  8:107413/1‑149 (MQ=255)
ttGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGctc                                                                                                                                                >  2:342859/1‑149 (MQ=255)
  ggTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGCGCtctc                                                                                                                                              >  8:57160/1‑149 (MQ=255)
  ggTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATAAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGCGCtctc                                                                                                                                              >  4:104094/1‑149 (MQ=255)
   gTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCtctc                                                                                                                                              <  1:43589/148‑1 (MQ=255)
       ccAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGa                                                                                                                                          <  8:303376/148‑1 (MQ=255)
       ccAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGAc                                                                                                                                         >  2:348440/1‑149 (MQ=255)
          gTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGa                                                                                                                                      >  2:325088/1‑149 (MQ=255)
              cgcgGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGatat                                                                                                                                   <  1:342859/148‑1 (MQ=255)
                cgGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGCGCTCTCCTGACGGATATCaa                                                                                                                                <  7:57160/149‑1 (MQ=255)
                cgGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATAAGGAACGCCTGTTTGTCGGCCTCTGACCAGCGTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGCGCTCTCCTGACGGATATCaa                                                                                                                                <  3:104094/149‑1 (MQ=255)
                 ggCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCaa                                                                                                                                >  8:19976/1‑148 (MQ=255)
                 ggCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGGCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAAt                                                                                                                               <  3:217355/149‑1 (MQ=255)
                  gCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATa                                                                                                                              <  5:48035/149‑1 (MQ=255)
                           gACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGc                                                                                                                     <  8:312905/149‑1 (MQ=255)
                                ttGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTAc                                                                                                                >  8:38572/1‑149 (MQ=255)
                                      ttcagttcGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGCCCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGCTATCAATACGGCCCAGCGGTACGCCGa                                                                                                           >  4:271799/1‑148 (MQ=255)
                                        cagttcGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGaaa                                                                                                         >  6:99641/1‑148 (MQ=255)
                                         agttcGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACa                                                                                                       >  6:354512/1‑149 (MQ=255)
                                          gttcGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACAt                                                                                                      >  7:340007/1‑149 (MQ=255)
                                                   tCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACgcgg                                                                                             >  6:97389/1‑149 (MQ=255)
                                                   tCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACgcgg                                                                                             <  5:99641/149‑1 (MQ=255)
                                                   tCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAAATAGCTGCTCATATCGAACATAGCTTAGGTATGACGCGTGCTCTCCTGACCGAAATCAATACCGACCAGCGGTACGCCTAAACATTTAAagcgc                                                                                             >  2:60787/1‑148 (MQ=255)
                                                      ggAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACgcggcg                                                                                           >  2:361004/1‑148 (MQ=255)
                                                      ggAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACgcggcg                                                                                           >  7:161900/1‑148 (MQ=255)
                                                       gcaCGCCTGTTTGTCGGCCTCTGCCCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACgcggcg                                                                                           >  4:299295/3‑147 (MQ=255)
                                                       gAACGCCAGTTTGTCGGCCTCTGACCAGCTTTCGTACTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACgcggcg                                                                                           >  1:149874/1‑147 (MQ=255)
                                                         aCGCCTGTTTGTCGGCCTCTGCCCAGCTTTCGTAGGCGCGGATACCGAGGGAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAgg                                                                                       <  7:181819/149‑1 (MQ=255)
                                                         aCGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAgg                                                                                       <  1:348440/149‑1 (MQ=255)
                                                         aCGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAgg                                                                                       <  8:67693/149‑1 (MQ=255)
                                                         aCGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAgg                                                                                       <  1:361004/149‑1 (MQ=255)
                                                         aCGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAgg                                                                                       <  3:299295/149‑1 (MQ=255)
                                                           gCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGgt                                                                                     >  4:156729/1‑149 (MQ=255)
                                                             cTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGgtgt                                                                                   >  5:349915/1‑149 (MQ=255)
                                                             cTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGgtgt                                                                                   <  7:12953/149‑1 (MQ=255)
                                                                    tCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAg                                                                            <  1:325088/149‑1 (MQ=255)
                                                                      ggCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAg                                                                               >  4:270816/1‑144 (MQ=255)
                                                                      ggCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAg                                                                               >  4:91173/1‑144 (MQ=255)
                                                                      ggCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGAGCTTCGGTAAGACGCCTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCAAGCCGCGCAGGGTGTCGAg                                                                               >  7:231090/1‑144 (MQ=255)
                                                                        ccTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCg                                                                        <  4:137107/149‑1 (MQ=255)
                                                                         ctctGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGc                                                                       >  8:338978/1‑149 (MQ=255)
                                                                          tctGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGa                                                                           <  7:250630/144‑1 (MQ=255)
                                                                          tctGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCACCGGTACGCCGAAACAATTCACGCGGCGCAGGCTGTCGAGCAGa                                                                           >  8:250630/1‑144 (MQ=255)
                                                                            tGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAgg                                                                                       <  5:26340/130‑1 (MQ=255)
                                                                            tGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAgg                                                                                       >  6:26340/1‑130 (MQ=255)
                                                                              aCCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGCCGGATATCAATACGGACCAGCGGTACGCCGAACCATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTg                                                                  >  7:297147/1‑149 (MQ=255)
                                                                              aCCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGt                                                                    <  7:355692/147‑1 (MQ=255)
                                                                              aCCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGt                                                                    >  8:355692/1‑147 (MQ=255)
                                                                                  gCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGa                                                              <  8:180767/149‑1 (MQ=255)
                                                                                      tCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGataa                                                           <  7:247587/148‑1 (MQ=255)
                                                                                            tCGCCGATACCGAGGTAGCGTGTCAGCTCGCCCAGCGCTTCGGTATGACGCGGGCTCTGCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCAt                                                    <  1:69445/149‑1 (MQ=255)
                                                                                              gCCGATACCGAGGTAGCGGGTCAGCTCGAACAGCGCTTCGGTATCACGCTTGCTCTCCTGACGGAAATCAATACGGACCAGCGGTACGCCGAAACATTTCACCCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACACATGc                                                  >  8:280056/1‑149 (MQ=255)
                                                                                                          gTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTg                                        >  2:312379/1‑147 (MQ=255)
                                                                                                               gggTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGAc                                 <  7:19976/149‑1 (MQ=255)
                                                                                                                     gCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAg                           <  5:354512/149‑1 (MQ=255)
                                                                                                                       tCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACCCGCCTGAAGGGACTGGTAGc                          >  5:272386/1‑148 (MQ=255)
                                                                                                                           ccAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCg                     >  1:19089/1‑149 (MQ=255)
                                                                                                                           ccAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCg                     >  2:50297/1‑149 (MQ=255)
                                                                                                                           ccAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGACACACGCCTGAAGTGACTGGAAGACAGCg                     >  6:88222/1‑149 (MQ=255)
                                                                                                                            cAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGt                    <  8:161900/149‑1 (MQ=255)
                                                                                                                             aGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCg                     >  7:95574/1‑147 (MQ=255)
                                                                                                                             aGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTa                   <  7:338978/149‑1 (MQ=255)
                                                                                                                                  ttCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCg              >  5:48398/1‑149 (MQ=255)
                                                                                                                                     ggTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCGGtt           <  7:38572/149‑1 (MQ=255)
                                                                                                                                         tGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCGGTTCCca       <  3:156729/149‑1 (MQ=255)
                                                                                                                                          gACGCGTGCTCTCCTGCCGGATATCAATACGGACCAGCGGTACGCCGCAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCGGTTCCca       >  7:373133/1‑148 (MQ=255)
                                                                                                                                             gcgTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCGGTTCCCACAGt   >  3:249844/1‑149 (MQ=255)
                                                                                                                                              cgTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCGGTTCCCACAGtt  <  7:107413/149‑1 (MQ=255)
                                                                                                                                               gTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCGGTTCCCACAg    <  6:193465/146‑1 (MQ=255)
                                                                                                                                               gTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCGGTTCCCACAg    >  5:193465/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TTGGTTGCCAGTTGCGCGGCAGAAGCGGACGTTTGGAGTTCAGTTCGCGGATCAGGAACGCCTGTTTGTCGGCCTCTGACCAGCTTTCGTAGTCGCCGATACCGAGGTAGCGGGTCAGCTCGCCCAGCGCTTCGGTATGACGCGTGCTCTCCTGACGGATATCAATACGGACCAGCGGTACGCCGAAACATTTCACGCGGCGCAGGGTGTCGAGCAGATCGCCGTTGGCGATAATACCCATGCCACACGCCTGAAGTGACTGGTAGCAAGCGTAGAGCGGTTCCCACAGTT  >  NC_000913/4151757‑4152047

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