Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,873,695 A→G 5.1% E41G (GAA→GGA) ‡ yeaK → mischarged aminoacyl‑tRNA deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,873,6950AG 5.1% 243.6 / 2.7 78E41G (GAA→GGA) ‡yeaKmischarged aminoacyl‑tRNA deacylase
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base G (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

TGTTTGCTTGTAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCGCCAGCCATATAGGCGGATTACGC  >  NC_000913/1873546‑1873843
                                                                                                                                                     |                                                                                                                                                    
tgttTGCTTGTAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGa                                                                                                                                                      <  1:208417/150‑1 (MQ=255)
 gttTGCTTGTAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGaa                                                                                                                                                     >  2:685088/1‑150 (MQ=255)
        tGTAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCTTGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTg                                                                                                                                              <  2:1083570/150‑1 (MQ=255)
         gTAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTgg                                                                                                                                             >  2:468658/1‑150 (MQ=255)
           aaTCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCa                                                                                                                                           >  1:262803/1‑150 (MQ=255)
           aaTCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCa                                                                                                                                           >  2:343814/1‑150 (MQ=255)
            aTCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCAc                                                                                                                                          >  2:553763/1‑150 (MQ=255)
             tctcAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCAcc                                                                                                                                         >  1:273301/1‑150 (MQ=255)
                 aGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGCCACGAATTTCTGACACCGctt                                                                                                                                     <  1:505487/150‑2 (MQ=255)
                 aGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGCCACGAATTTCTGACACCGctt                                                                                                                                     <  1:345064/150‑2 (MQ=255)
                              gACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGc                                                                                                                        >  2:460544/1‑150 (MQ=255)
                              gACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGc                                                                                                                        <  2:181849/150‑1 (MQ=255)
                                cTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTa                                                                                                                      >  1:1355322/1‑150 (MQ=255)
                                cTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTa                                                                                                                      >  2:627892/1‑150 (MQ=255)
                                  gAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAg                                                                                                                    <  2:281152/150‑1 (MQ=255)
                                     aTGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCa                                                                                                                 <  1:623672/150‑1 (MQ=255)
                                       ggCTAAAGGAAGGGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGATGCGTTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACt                                                                                                               <  2:598472/150‑1 (MQ=255)
                                       ggCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACt                                                                                                               >  1:1319883/1‑150 (MQ=255)
                                       ggCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACt                                                                                                               >  2:205267/1‑150 (MQ=255)
                                           aaaGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTa                                                                                                           <  1:358332/150‑1 (MQ=255)
                                           aaaGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTa                                                                                                           <  1:39798/150‑1 (MQ=255)
                                             aGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATg                                                                                                         >  1:1021736/1‑150 (MQ=255)
                                             aGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATg                                                                                                         <  2:306785/150‑1 (MQ=255)
                                              ggAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCt                                                                                                                       <  2:1008530/135‑1 (MQ=255)
                                              ggAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCt                                                                                                                       >  1:1008530/1‑135 (MQ=255)
                                              ggAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGt                                                                                                        <  1:252314/150‑1 (MQ=255)
                                              ggAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGt                                                                                                        >  2:102883/1‑150 (MQ=255)
                                                  gCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAg                                                                                                    >  1:1104148/1‑150 (MQ=255)
                                                     tGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCa                                                                                                 <  1:102883/150‑1 (MQ=255)
                                                          cATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGc                                                                                            <  2:118895/150‑1 (MQ=255)
                                                            tCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCaa                                                                                          >  1:768642/1‑150 (MQ=255)
                                                             aaGCGGTTAATCGTTTTATTATCACATTAGGTTCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATg                                                                                        <  1:1123222/149‑1 (MQ=255)
                                                               gCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATgg                                                                                       >  2:1212508/1‑150 (MQ=255)
                                                               gCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATgg                                                                                       >  2:1005409/1‑150 (MQ=255)
                                                                  gTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGt                                                                                    >  1:1422090/1‑150 (MQ=255)
                                                                   ttAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTa                                                                                   <  1:685088/150‑1 (MQ=255)
                                                                    tAATCGCTTTATTATCACAAGAAGGTGATGACTTCCGAGTCGTTACCCATGAAGCCGAAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACAGCTCTAGGCCAGGGTGCTAATGCACTGGTATGTAAAGTCAAAGGCAATGGCGTaa                                                                                  >  2:1406342/1‑150 (MQ=255)
                                                                         gCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGc                                                                                                                        <  1:792631/107‑1 (MQ=255)
                                                                         gCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGc                                                                                                                        >  2:792631/1‑107 (MQ=255)
                                                                         gCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCaa                                                                             <  2:177641/150‑1 (MQ=255)
                                                                          cTTTATTATCACAAGAATGTGCTTACTTCAGCGTCGTTACCCATTCAGCTGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGAACCGCTATAGGCCAGTGTGCTAAGTCACTGGTATGTAAAGTCAAAGGCAATGGCTTAAATCAAc                                                                            <  2:1345839/150‑1 (MQ=255)
                                                                          cTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAAc                                                                            <  2:1372532/150‑1 (MQ=255)
                                                                              attatCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGt                                                                        <  1:1111334/150‑1 (MQ=255)
                                                                                 atCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCt                                                                     >  2:170231/1‑150 (MQ=255)
                                                                                   cacaAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTgg                                                                   <  2:615409/150‑1 (MQ=255)
                                                                                    acaAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGc                                                                  <  2:1355322/150‑1 (MQ=255)
                                                                                      aagaagGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGa                                                                <  1:1063061/150‑1 (MQ=255)
                                                                                       agaagGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGAt                                                               >  1:1198447/1‑150 (MQ=255)
                                                                                       agaagGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGAt                                                               >  1:821536/1‑150 (MQ=255)
                                                                                          agGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCt                                                            >  1:975150/1‑150 (MQ=255)
                                                                                            gTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCg                                                          <  1:1376076/150‑1 (MQ=255)
                                                                                             tGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGc                                                         >  1:129865/1‑150 (MQ=255)
                                                                                                 gACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGcc                                                     <  1:1212508/150‑1 (MQ=255)
                                                                                                 gACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGcc                                                     <  1:343814/150‑1 (MQ=255)
                                                                                                   cTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGa                                                   >  1:484332/1‑150 (MQ=255)
                                                                                                        gcgTCGTTACCCATGTAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAAc                                              <  2:1104148/150‑1 (MQ=255)
                                                                                                        gcgTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAAc                                              <  1:460544/150‑1 (MQ=255)
                                                                                                            cgtTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGc                                          <  1:111779/150‑1 (MQ=255)
                                                                                                                aCCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGAt                                      <  2:821536/150‑1 (MQ=255)
                                                                                                                 cccATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATc                                     <  1:170231/150‑1 (MQ=255)
                                                                                                                 cccATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATc                                     >  2:111810/1‑150 (MQ=255)
                                                                                                                  ccATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATTTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGTCCGATCt                                    <  2:1198447/150‑1 (MQ=255)
                                                                                                                  ccATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCt                                    >  1:1163604/1‑150 (MQ=255)
                                                                                                                       aaGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGcc                               <  1:472604/150‑1 (MQ=255)
                                                                                                                       aaGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGcc                               <  1:468658/150‑1 (MQ=255)
                                                                                                                          cGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAAc                            <  2:273301/150‑1 (MQ=255)
                                                                                                                           ggtaggtaAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACt                           >  2:606619/1‑150 (MQ=255)
                                                                                                                               ggtaAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCgcc                       >  2:273562/1‑150 (MQ=255)
                                                                                                                                 taAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCgccag                     <  1:627892/150‑1 (MQ=255)
                                                                                                                                     tgtgAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATc                                                 >  2:303066/1‑118 (MQ=255)
                                                                                                                                     tgtgAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATc                                                 <  1:303066/118‑1 (MQ=255)
                                                                                                                                     tgtgAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGTCAAc                            <  1:1194264/139‑1 (MQ=255)
                                                                                                                                     tgtgAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGTCAAc                            >  2:1194264/1‑139 (MQ=255)
                                                                                                                                     tgtgAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCGCCAGCCat                 >  1:57604/1‑150 (MQ=255)
                                                                                                                                     tgtgAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCGCCAGCCat                 >  1:307260/1‑150 (MQ=255)
                                                                                                                                     cctagaGCGGTGCCACGAATTTCTGACACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCGCCAGCCat                 >  1:1124125/6‑150 (MQ=255)
                                                                                                                                        agaGCGGTGCCACGAATTTCTGACACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCGCCAGCCAtata              >  2:391685/3‑150 (MQ=255)
                                                                                                                                         aaGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCGCCAGCCATATAg             <  1:606619/150‑1 (MQ=255)
                                                                                                                                              gtgtCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCGCCAGCCATATAGGCGGa        >  2:1241903/1‑150 (MQ=255)
                                                                                                                                                    gAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCGCCAGCCATATAGGCGGATTAcgc  <  2:768642/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGTTTGCTTGTAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAACAGGCCGATCTGAGCCAACTCGCCAGCCATATAGGCGGATTACGC  >  NC_000913/1873546‑1873843

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