Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,259,126 C→A 5.3% G175V (GGT→GTT)  tdcG ← L‑serine dehydratase 3, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,259,1260CA 5.3% 377.4 / 5.9 122G175V (GGT→GTT) tdcGL‑serine dehydratase 3, anaerobic
Reads supporting (aligned to +/- strand):  ref base C (52/63);  new base A (3/4);  total (55/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCG  >  NC_000913/3258980‑3259275
                                                                                                                                                  |                                                                                                                                                     
tCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACctg                                                                                                                                                    >  2:346781/1‑150 (MQ=255)
 cccACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACctgc                                                                                                                                                   >  1:579108/1‑150 (MQ=255)
   cACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACctgctg                                                                                                                                                 <  2:890566/150‑1 (MQ=255)
   cACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACctgctg                                                                                                                                                 >  1:160697/1‑150 (MQ=255)
   cACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGGTGAActgctg                                                                                                                                                 <  2:890666/150‑1 (MQ=255)
     cGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGGTGAACTGCTGaa                                                                                                                                               <  2:828527/150‑2 (MQ=255)
      gTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGt                                                                                                                                              >  2:758258/1‑150 (MQ=255)
       ttCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTg                                                                                                                                             <  2:508056/150‑1 (MQ=255)
        tCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTgg                                                                                                                                            >  2:1206596/1‑150 (MQ=255)
        tCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGGTGAACTGCTGaaaat                                                                                                                                            <  1:798353/150‑5 (MQ=255)
           aTACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGaaa                                                                                                                                         >  1:1613238/1‑150 (MQ=255)
           aTACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGGTGAACTGCTGaaaatgtg                                                                                                                                         <  1:572557/150‑8 (MQ=255)
              ccGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCg                                                                                                                                      >  2:308172/1‑150 (MQ=255)
              ccGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCg                                                                                                                                      >  1:656498/1‑150 (MQ=255)
              ccGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCg                                                                                                                                      <  1:1509725/150‑1 (MQ=255)
                gTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTa                                                                                                                                    >  1:5078/1‑150 (MQ=255)
                  cATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAg                                                                                                                                  <  1:1221446/150‑1 (MQ=255)
                  cATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATAAGACCAGATAAAGACAGGCCCTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAg                                                                                                                                  >  2:1570662/1‑150 (MQ=255)
                    tGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGt                                                                                                                                >  2:1333877/1‑150 (MQ=255)
                    tGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAAGCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGt                                                                                                                                >  1:147157/1‑150 (MQ=255)
                     gCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTa                                                                                                                               >  2:1397961/1‑150 (MQ=255)
                      cATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTAc                                                                                                                              <  2:1578430/150‑1 (MQ=255)
                             tGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGtt                                                                                                                       >  2:491270/1‑150 (MQ=255)
                              gCCAGATACGGGCAAAAGCTGCGTCATTTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGGATCTTCAGACCAGATATTGACAGGCCGTTGTAATCACAGATTTTCTGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGttt                                                                                                                      <  2:793225/150‑1 (MQ=255)
                                ccgctACGGGCAAAACCGGCTGTAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTTAGGTACGGACGTTTCg                                                                                                                    <  2:1059579/146‑1 (MQ=255)
                                cAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCg                                                                                                                    <  2:130847/150‑1 (MQ=255)
                                cAGATACGGGCAAAACCGGCGGCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCg                                                                                                                    <  2:152625/150‑1 (MQ=255)
                                 aGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACTGACGTTTCGa                                                                                                                   >  2:638995/1‑150 (MQ=255)
                                 aGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGa                                                                                                                   <  1:1584980/150‑1 (MQ=255)
                                  gATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGAc                                                                                                                  >  1:912955/1‑150 (MQ=255)
                                  gATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGAc                                                                                                                  <  2:649622/150‑1 (MQ=255)
                                        ggCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTg                                                                                                            >  1:47495/1‑150 (MQ=255)
                                            aaaCCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGAc                                                                                                        <  1:9753/150‑1 (MQ=255)
                                             aaCCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACa                                                                                                       <  1:935803/150‑1 (MQ=255)
                                                    gTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGaa                                                                                                >  2:1622493/1‑150 (MQ=255)
                                                     tgaaTTTCCGCTTTGCTGCGCAGCGCCAGCTCGTTGTGCATCATCAGCCCAGACATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAg                                                                                               <  2:1555977/148‑1 (MQ=255)
                                                          ttCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATAATCACACCAGATATATACAGGACGATGTAATCACACATTTTCAGCAGATCACCTCCTCAGTGGAAATCGTAAGGTCCGGCCGTTACGACAACGTGCGACAGGCCGAAGAGttc                                                                                          >  1:735849/1‑150 (MQ=255)
                                                           tCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttct                                                                                         <  2:1074886/150‑1 (MQ=255)
                                                           tCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttct                                                                                         <  2:1710025/150‑1 (MQ=255)
                                                           tCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttct                                                                                         <  1:1405285/150‑1 (MQ=255)
                                                              gCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttcttct                                                                                      <  2:1638900/150‑1 (MQ=255)
                                                              gCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttcttct                                                                                      <  1:247372/150‑1 (MQ=255)
                                                                tttGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttcttcttc                                                                                    <  2:1017732/150‑1 (MQ=255)
                                                                 ttGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCg                                                                                   <  1:1163534/150‑1 (MQ=255)
                                                                 ttGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCg                                                                                   <  2:421091/150‑1 (MQ=255)
                                                                 ttGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCATAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTACg                                                                                   >  2:1536009/1‑150 (MQ=255)
                                                                  tgctgcGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGa                                                                                  >  2:734690/1‑150 (MQ=255)
                                                                     tgcGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCTGAAGTGTTCTTCTTCGACaa                                                                               <  2:463997/150‑1 (MQ=255)
                                                                      gcgcAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAAt                                                                              <  1:631598/150‑1 (MQ=255)
                                                                         cAGCGCTAGCTCGTTGTGCATCATCAGCCCAGATATAGACAGGCCGTTGTAATCACAGATTTTCAGCAGGTCACCGGCTGAGTGGAACTCGTAAGGTGCGGCCGTTTCGCCATGGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATaac                                                                           >  2:270255/1‑149 (MQ=255)
                                                                              cTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACccgc                                                                      <  2:848621/150‑1 (MQ=255)
                                                                                 gCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACccgccgc                                                                   >  1:150573/1‑150 (MQ=255)
                                                                                   tCGTTGTGCATCATCAGACCAGATATAGACAGGCCGGTGTAATCACACATTTTCCGCAGATCACCTGCTGAGTGGAAATCGTCAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACccgccgccg                                                                 <  1:107736/150‑1 (MQ=255)
                                                                                    cGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGa                                                                <  1:1120822/150‑1 (MQ=255)
                                                                                          gCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGt                                                          <  2:324147/150‑1 (MQ=255)
                                                                                              catcaGCCCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTAACCTGCTGAGTGGAAATAGTAAGGTACGGACGTTTCGACATCGTGCGACACGCCGAAGTGTTGTTCTTCGACAATACACCCGCCGCCGACAGAGTAAtt                                                      >  2:822640/1‑149 (MQ=255)
                                                                                              catcaGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttcttcttc                                                                                    <  1:1468941/120‑1 (MQ=255)
                                                                                              catcaGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttcttcttc                                                                                    >  2:1468941/1‑120 (MQ=255)
                                                                                               atcaGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAg                                                     <  1:592286/150‑1 (MQ=255)
                                                                                                tcaGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAgg                                                    >  2:709944/1‑150 (MQ=255)
                                                                                                 caGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGt                                                   >  2:1084911/1‑150 (MQ=255)
                                                                                                  aGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGtt                                                  <  2:719921/150‑1 (MQ=255)
                                                                                                   gACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACTGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCTACAGAGTAATAGGttt                                                 >  1:1628503/1‑150 (MQ=255)
                                                                                                   gACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGttt                                                 >  1:1681417/1‑150 (MQ=255)
                                                                                                     ccAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTa                                               >  2:1137730/1‑150 (MQ=255)
                                                                                                      cAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTAc                                              <  1:1397961/150‑1 (MQ=255)
                                                                                                        gATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACtt                                            <  2:579108/150‑1 (MQ=255)
                                                                                                           ataGACAGGCCGTTGTAATCACAAATTTTCAGCAGTTCACCTGCTGAGTGTAAATCGTAAGGTAAGGACGTTGCGACATAGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAAt                                         <  1:823577/150‑1 (MQ=255)
                                                                                                              gACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGc                                      <  1:758258/150‑1 (MQ=255)
                                                                                                                    ccGTTGTAATCACACATTTTCAGCAGTTCGCCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTcc                                <  1:338293/150‑1 (MQ=255)
                                                                                                                     cGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCt                               <  2:902484/150‑1 (MQ=255)
                                                                                                                      gTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTg                              >  2:698773/1‑150 (MQ=255)
                                                                                                                        tGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTc                            <  1:346781/150‑1 (MQ=255)
                                                                                                                           aaTCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCt                         <  2:642779/150‑1 (MQ=255)
                                                                                                                            aTCACACCTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGGGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCtt                        <  2:369535/150‑1 (MQ=255)
                                                                                                                            aTCACACAGTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGATGCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCtt                        <  2:1133132/150‑1 (MQ=255)
                                                                                                                               acacaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCa                     <  1:1137730/150‑1 (MQ=255)
                                                                                                                               acacaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCa                     <  1:1439716/150‑1 (MQ=255)
                                                                                                                               acacaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTATTCCTGTCCCTTCCa                     >  2:442650/1‑150 (MQ=255)
                                                                                                                                 ttccactCAGCAGGTGAACTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCAGCCGCCGACAGAGTAATAGGTGTTACTTAATAGCTCTTCCTGTCCCTTCCAgg                   >  1:1062015/7‑150 (MQ=255)
                                                                                                                                 acaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAgg                   <  2:419950/150‑1 (MQ=255)
                                                                                                                                  caTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGc                  >  2:1598761/1‑150 (MQ=255)
                                                                                                                                  caTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGc                  >  2:1569980/1‑150 (MQ=255)
                                                                                                                                   ccactCAGCAGGTGAACTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCa                 >  1:860077/5‑150 (MQ=255)
                                                                                                                                    ttttCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAATCCCGCCGCCGACTGAGTAATAGGTTTTAATTAATAGCTCTTCCTGTCCCTTCCAGGCAg                >  1:252597/1‑150 (MQ=255)
                                                                                                                                    ttttCAGCAGTTCACCTGCTGAGTGGAAATAGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAAGAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAg                >  1:256258/1‑150 (MQ=255)
                                                                                                                                     tttCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGt               <  1:67450/150‑1 (MQ=255)
                                                                                                                                      ttCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTg              >  2:1380756/1‑150 (MQ=255)
                                                                                                                                      ttCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTg              >  1:1375148/1‑150 (MQ=255)
                                                                                                                                        cagcagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGAt            >  1:623933/1‑150 (MQ=255)
                                                                                                                                           cagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATcca         <  1:1206596/150‑2 (MQ=255)
                                                                                                                                           cagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCg         <  2:108525/150‑1 (MQ=255)
                                                                                                                                           cagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCg         >  1:866640/1‑150 (MQ=255)
                                                                                                                                           cagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCg         <  1:601863/150‑1 (MQ=255)
                                                                                                                                            cgTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACCTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACATTAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGc        <  2:331677/149‑1 (MQ=255)
                                                                                                                                            agTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGc        >  1:1299047/1‑150 (MQ=255)
                                                                                                                                              ttCACCTGCTGAGTGGAAATCGTCAGGGACGGACGTTTCGACATCGTGCGACTGGCCGAAGTGTACTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTGTACTTAATAGCTCTACCTGGCCCTTCCAGGCAGTGAGCAGCAt      <  2:1313397/150‑1 (MQ=255)
                                                                                                                                                  ccTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGCCAATAAACCCGCCGCCGACAGAGTAAGAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCCTTCCg  >  2:1431155/1‑150 (MQ=255)
                                                                                                                                                  ccTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCg  >  2:1378115/1‑150 (MQ=255)
                                                                                                                                                  ccTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCg  >  2:1659644/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCG  >  NC_000913/3258980‑3259275

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