Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,012,830 A→G T730A (ACG→GCG) 
V297A (GTT→GCT) 		 metE →
ysgA ←		 tetrahydropteroyltriglutamate methyltransferase
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,012,8300AG84.0% 723.9 / 108.6 288T730A (ACG→GCG) 
V297A (GTT→GCT) 		metE
ysgA		tetrahydropteroyltriglutamate methyltransferase
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/23);  new base G (134/108);  total (157/131)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 304 total aligned reads displayed.

CGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTG  >  NC_000913/4012738‑4012919
                                                                                            |                                                                                         
  aGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGcg                                                                                         >  2:5003211/1‑93 (MQ=255)
  aGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAAcg                                                                                         <  1:167067/93‑1 (MQ=255)
   gCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAgcgc                                                                                        >  2:2289402/1‑93 (MQ=255)
   gCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAAcgc                                                                                        <  1:6642341/93‑1 (MQ=255)
     gTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAgcgcgc                                                                                      <  2:5361590/93‑1 (MQ=255)
       ggAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCgg                                                                                    >  1:6462510/1‑93 (MQ=255)
       ggAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCgg                                                                                    >  1:1740361/1‑93 (MQ=255)
         aaTGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCt                                                                                  >  2:869563/1‑93 (MQ=255)
         aaTGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCt                                                                                  >  1:4967005/1‑93 (MQ=255)
         aaTGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCt                                                                                  >  2:5611224/1‑93 (MQ=255)
          aTGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTg                                                                                 <  2:5239313/93‑1 (MQ=255)
          aTGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTg                                                                                 <  1:1999285/93‑1 (MQ=255)
             gATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGcc                                                                              <  1:2642038/93‑1 (MQ=255)
             gATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGcc                                                                              >  2:6303386/1‑93 (MQ=255)
              aTTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCTCCGCTGGCCa                                                                             >  1:4745119/1‑93 (MQ=255)
               ttGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAg                                                                            <  2:4835282/93‑1 (MQ=255)
               ttGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAg                                                                            <  2:5431519/93‑1 (MQ=255)
                 gAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGaa                                                                          >  2:3990695/1‑93 (MQ=255)
                  aaGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGaaa                                                                         >  1:2953941/1‑93 (MQ=255)
                  aaGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGAGCGCGGCTGGCCAGaaa                                                                         >  1:3342508/1‑93 (MQ=255)
                    gCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAAcc                                                                       >  2:3896999/1‑93 (MQ=255)
                     ccTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAAccc                                                                      >  1:4140671/1‑93 (MQ=255)
                     ccTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAAccc                                                                      >  2:2685862/1‑93 (MQ=255)
                     ccTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAAccc                                                                      >  1:5426076/1‑93 (MQ=255)
                      cTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCcg                                                                     <  1:4606698/93‑1 (MQ=255)
                       ttGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCCGAAAACcgc                                                                    >  2:3369528/1‑93 (MQ=255)
                       ttGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCcgc                                                                    <  1:6078705/93‑1 (MQ=255)
                        tGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCcgcg                                                                   >  1:4561567/1‑93 (MQ=255)
                        tGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCcgcg                                                                   >  1:3171190/1‑93 (MQ=255)
                          cTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCcgcgcg                                                                 <  1:566665/93‑1 (MQ=255)
                           tGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCgg                                                                >  1:5322644/1‑93 (MQ=255)
                            gaagaaAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGc                                                               >  2:1470229/1‑93 (MQ=255)
                             aagaaAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCa                                                              >  2:947025/1‑93 (MQ=255)
                              agaaAGCGGGAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCAc                                                             <  1:1178820/93‑1 (MQ=255)
                               gaaAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACt                                                            >  1:1391171/1‑93 (MQ=255)
                               gaaAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACt                                                            <  1:5859875/93‑1 (MQ=255)
                                 aaGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTgg                                                          >  2:2908121/1‑93 (MQ=255)
                                 aaGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTgg                                                          <  2:3562302/93‑1 (MQ=255)
                                 aaGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTgg                                                          >  2:1405174/1‑93 (MQ=255)
                                 aaGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTgg                                                          >  2:982327/1‑93 (MQ=255)
                                 aaGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTgg                                                          >  1:1915773/1‑93 (MQ=255)
                                   gCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCg                                                        >  2:3516503/1‑93 (MQ=255)
                                    cGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGa                                                       >  2:3191432/1‑93 (MQ=255)
                                     ggCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGaa                                                      <  2:3799543/93‑1 (MQ=255)
                                      gCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAAc                                                     >  2:1224892/1‑93 (MQ=255)
                                       cAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACa                                                    <  2:3684883/93‑1 (MQ=255)
                                        aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACAt                                                   >  1:389714/1‑93 (MQ=255)
                                        aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACAt                                                   >  1:2065750/1‑93 (MQ=255)
                                         aaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATg                                                  <  2:240897/93‑1 (MQ=255)
                                         aaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATg                                                  <  2:6118870/93‑1 (MQ=255)
                                           aCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGt                                                <  1:5795501/93‑1 (MQ=255)
                                            cGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTg                                               <  2:6079828/93‑1 (MQ=255)
                                            cGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTg                                               <  2:3622115/93‑1 (MQ=255)
                                            cGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTg                                               >  2:5541662/1‑93 (MQ=255)
                                            cGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCACGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTg                                               <  2:6647706/93‑1 (MQ=255)
                                            cGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGtc                                               >  2:5385415/1‑92 (MQ=255)
                                            cGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTg                                               <  2:1067764/93‑1 (MQ=255)
                                             gCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGc                                              >  2:6527832/1‑93 (MQ=255)
                                              cATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCa                                             <  2:2778085/93‑1 (MQ=255)
                                               aTTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAg                                            >  2:2550428/1‑93 (MQ=255)
                                                  ccGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAggcg                                         <  1:2897355/93‑1 (MQ=255)
                                                  ccGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAggcg                                         <  1:30390/93‑1 (MQ=255)
                                                    ggCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAggcggc                                       <  2:1343860/93‑1 (MQ=255)
                                                       agagCGCCTGTGGGTCAACCCTGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCa                                    <  1:5148675/93‑1 (MQ=255)
                                                         agCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGa                                  >  2:432933/1‑93 (MQ=255)
                                                         agCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGa                                  <  2:6365620/93‑1 (MQ=255)
                                                          gcgcCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGaa                                 >  2:184034/1‑93 (MQ=255)
                                                          gcgcCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGaa                                 <  1:6328954/93‑1 (MQ=255)
                                                            gcCTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACt                               >  1:3618905/1‑93 (MQ=255)
                                                             ccTGTGGGTCACCCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACtt                              <  1:2335312/93‑1 (MQ=255)
                                                             ccTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACtt                              <  1:5484618/93‑1 (MQ=255)
                                                             ccTGTGGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACtt                              >  2:4393399/1‑93 (MQ=255)
                                                             ccTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACtt                              >  1:343200/1‑93 (MQ=255)
                                                               tgtgGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGc                            >  1:5020530/1‑93 (MQ=255)
                                                                gtgGGTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCg                           >  1:3812862/1‑93 (MQ=255)
                                                                  gggTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTc                         >  1:2786755/1‑93 (MQ=255)
                                                                    gTCAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCgg                       <  2:5071080/93‑1 (MQ=255)
                                                                    gTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCgg                       <  1:4656288/93‑1 (MQ=255)
                                                                      cAACCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCgggg                     <  2:3125703/93‑1 (MQ=255)
                                                                      cAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCgggg                     >  2:5680294/1‑93 (MQ=255)
                                                                       aaCCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCggggg                    >  2:4629120/1‑93 (MQ=255)
                                                                        aCCCGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGt                   >  1:954891/1‑93 (MQ=255)
                                                                         cccGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTa                  <  2:2061740/93‑1 (MQ=255)
                                                                          ccGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTaa                 >  1:1061968/1‑93 (MQ=255)
                                                                           cGGACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCACAACTTGCGTCGGGGGTaaa                >  1:6583195/1‑93 (MQ=255)
                                                                           cGGACTGTGGCCTGAAAGCGCGCGGCTGGCAAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTaaa                <  1:1263396/93‑1 (MQ=255)
                                                                           cGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGGaaa                >  2:1723137/1‑93 (MQ=255)
                                                                            ggACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTaaaa               >  2:1036615/1‑93 (MQ=255)
                                                                             gACTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAAt              >  2:5937937/1‑93 (MQ=255)
                                                                              aCTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATc             <  2:1532525/93‑1 (MQ=255)
                                                                               cTGTGGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATcc            <  1:909195/93‑1 (MQ=255)
                                                                                tgtgGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCa           <  2:454226/93‑1 (MQ=255)
                                                                                  tgGCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCaaa         <  1:2114019/93‑1 (MQ=255)
                                                                                   ggCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAAc        <  2:4900418/93‑1 (MQ=255)
                                                                                    gCCTGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAAcc       <  2:5859210/93‑1 (MQ=255)
                                                                                    gCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAAcc       >  2:510937/1‑93 (MQ=255)
                                                                                       tGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCggg    >  2:6246331/1‑93 (MQ=255)
                                                                                       tGAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCggg    >  1:2398928/1‑93 (MQ=255)
                                                                                        gAAAGCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGggt   <  1:3124538/93‑1 (MQ=255)
                                                                                         aaagCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGggtg  >  1:5370883/1‑93 (MQ=255)
                                                                                            |                                                                                         
CGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTG  >  NC_000913/4012738‑4012919

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑