Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,356,571 G→A 5.1% R616Q (CGG→CAG)  gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,356,5710GA 5.1% 265.0 / 2.8 78R616Q (CGG→CAG) gltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base A (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

GATGAAAGAGGAGCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACTTCGCCGTGTT  >  NC_000913/3356425‑3356718
                                                                                                                                                  |                                                                                                                                                   
gaTGAAAGAGGAGCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCgg                                                                                                                                                    <  3:28320/148‑1 (MQ=255)
gaTGAAAGAGGAGCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCgg                                                                                                                                                    <  3:32352/148‑1 (MQ=255)
   gAAAGAGGAGCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTCGCGGCACCGTGCTGCAGGTGCTCTCCGACCGGat                                                                                                                                                  >  8:305378/1‑146 (MQ=255)
      agagGAGCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGTAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATc                                                                                                                                              >  2:84263/1‑148 (MQ=255)
      agagGAGCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATc                                                                                                                                              <  4:228307/148‑1 (MQ=255)
       gaggagCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCg                                                                                                                                             <  3:152362/148‑1 (MQ=255)
          gagCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCt                                                                                                                                           >  8:304915/1‑147 (MQ=255)
           agCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTaa                                                                                                                                         <  3:251970/148‑1 (MQ=255)
              aCTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTa                                                                                                                                          >  1:221923/1‑144 (MQ=255)
              aCTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTa                                                                                                                                          >  2:129041/1‑144 (MQ=255)
              aCTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACAACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCAACTTGATGCTGTTGCTCTCCGACAAGAATATCGCTTAAg                                                                                                                                       >  1:259939/1‑147 (MQ=255)
               cTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGAt                                                                                                                                     <  7:106319/148‑1 (MQ=255)
                 aCCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAg                                                                                                                                       <  2:221923/144‑1 (MQ=255)
                 aCCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCg                                                                                                                                   >  3:130664/1‑148 (MQ=255)
                   cgcgCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGcc                                                                                                                                 >  4:18875/1‑148 (MQ=255)
                   cgcgCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGcc                                                                                                                                 <  5:140379/148‑1 (MQ=255)
                    gcgcAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCt                                                                                                                                <  1:115827/148‑1 (MQ=255)
                         gATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCAGAATATCGCTAAAGATCGCCTGCCgg                                                                                                                           >  8:32252/1‑148 (MQ=255)
                           tACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGtt                                                                                                                         <  1:310406/148‑1 (MQ=255)
                                    tatCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCg                                                                                                                <  4:233717/148‑1 (MQ=255)
                                         ccTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGc                                                                                                           >  2:48393/1‑148 (MQ=255)
                                          cTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCg                                                                                                          >  5:2549/1‑148 (MQ=255)
                                            ttGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGt                                                                                                        >  1:38083/1‑148 (MQ=255)
                                             tGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGt                                                                                                        <  7:305378/147‑1 (MQ=255)
                                              gACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGtt                                                                                                       <  8:23124/147‑1 (MQ=255)
                                                  tCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCACGATGGCGGTTGgcgc                                                                                                  >  3:292888/1‑148 (MQ=255)
                                                  tCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCACGATGGAGGTTGgcgc                                                                                                  >  6:92592/1‑148 (MQ=255)
                                                   cACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGa                                                                                                                                                   >  6:257899/1‑98 (MQ=255)
                                                   cACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGa                                                                                                                                                   <  5:257899/98‑1 (MQ=255)
                                                   cACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGgcgc                                                                                                  >  6:174285/1‑147 (MQ=255)
                                                    aCTAAAACCACGCTCGTAGCGACAGGCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGa                                                                                                <  1:84263/148‑1 (MQ=255)
                                                      tAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATc                                                                                              >  3:294258/1‑148 (MQ=255)
                                                       aaaaCCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATc                                                                                              >  2:181871/1‑147 (MQ=255)
                                                       aaaaCCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATc                                                                                              >  4:94942/1‑147 (MQ=255)
                                                        aaaCCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATcc                                                                                             <  3:94942/147‑1 (MQ=255)
                                                        aaaCCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATc                                                                                              <  7:5048/146‑1 (MQ=255)
                                                          aCCACGCTCGAAGCGACAGTCAACCAGCTGTGCGACAAAGCCGAAAACATCGTACGTAGCGGCACCGTGCTGCAGGCGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGGTCCAGCCCCGATGGCGGTTGGCGCGATCCAGa                                                                                          <  1:138472/148‑1 (MQ=255)
                                                          aCCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCa                                                                                                                      <  1:204279/120‑1 (MQ=255)
                                                          aCCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCa                                                                                                                      >  2:204279/1‑120 (MQ=255)
                                                            cACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGAcc                                                                                        >  3:246076/1‑148 (MQ=255)
                                                            cACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGAcc                                                                                        >  1:213427/1‑148 (MQ=255)
                                                              cGCTCTAAGCGACAGTCAAAGAGCTTTGCTACAAAGCCGAAAAAATGGTACGTAGCGTCACCGTGCTGCTGGTGCTCTCCGACCAGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCACCGATGTCGGTTGGCGCGATCCAGACCCg                                                                                      <  2:259939/148‑1 (MQ=255)
                                                               gCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGt                                                                                     <  1:173063/148‑1 (MQ=255)
                                                                  cGAAGCGACAGTCAAATAGCTGTGCGACAAAGACGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTg                                                                                  >  5:239877/1‑148 (MQ=255)
                                                                       cGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGa                                                                             <  1:48393/148‑1 (MQ=255)
                                                                            gTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCa                                                                          >  2:129680/1‑146 (MQ=255)
                                                                               aaaGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGcc                                                                     >  7:132748/1‑148 (MQ=255)
                                                                                atgagCTGGGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGGGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGTTCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGAGCAAAGCCt                                                                    <  4:292888/146‑1 (MQ=255)
                                                                                 agagCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTg                                                                   <  2:38083/148‑1 (MQ=255)
                                                                                    gCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGt                                                                <  5:174285/148‑1 (MQ=255)
                                                                                        tgCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTg                                                              <  2:56293/146‑1 (MQ=255)
                                                                                        tgCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCg                                                            <  4:246076/148‑1 (MQ=255)
                                                                                         gCGACAAAGCCGAAACAATGGTACGTAGCGGCACCGTGCGGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGa                                                           <  6:116500/148‑1 (MQ=255)
                                                                                          cGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGa                                                           <  7:304915/147‑1 (MQ=255)
                                                                                            aCAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGc                                                        >  4:19902/1‑148 (MQ=255)
                                                                                              aaaGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGa                                                           <  5:35046/143‑1 (MQ=255)
                                                                                              aaaGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGa                                                           >  6:35046/1‑143 (MQ=255)
                                                                                                 gCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAAca                                                   >  6:88788/1‑148 (MQ=255)
                                                                                                   cGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCAGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAAcatc                                                 <  7:32252/148‑1 (MQ=255)
                                                                                                    gAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAAcatca                                                >  4:144536/1‑148 (MQ=255)
                                                                                                     aaaaaaTGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAg                                                                       >  8:306685/1‑124 (MQ=255)
                                                                                                     aaaaaaTGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAg                                                                       <  7:306685/124‑1 (MQ=255)
                                                                                                     aaaaaaTGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAAcatcat                                               <  2:213427/148‑1 (MQ=255)
                                                                                                       aaaaTGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCAtcg                                             <  3:101855/148‑1 (MQ=255)
                                                                                                        aaaTGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCAtcgt                                            <  3:144536/148‑1 (MQ=255)
                                                                                                                   aGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGAt                                                                            >  6:253509/1‑105 (MQ=255)
                                                                                                                   aGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGAt                                                                            <  5:253509/105‑1 (MQ=255)
                                                                                                                                 gctgGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCg                   >  3:55772/1‑148 (MQ=255)
                                                                                                                                     gTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGca                 <  4:55772/146‑1 (MQ=255)
                                                                                                                                     gTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGca                 <  1:129680/146‑1 (MQ=255)
                                                                                                                                      tGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGCTCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGcacca              >  1:297663/1‑148 (MQ=255)
                                                                                                                                        ctctCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACt            <  1:129041/148‑1 (MQ=255)
                                                                                                                                          ctcCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGc                  <  1:97279/140‑1 (MQ=255)
                                                                                                                                          ctcCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGc                  >  2:97279/1‑140 (MQ=255)
                                                                                                                                           tcCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACtt           <  4:294258/146‑1 (MQ=255)
                                                                                                                                           tcCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACTTCg         >  1:45891/1‑148 (MQ=255)
                                                                                                                                            ccGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACTTCGc        >  5:191168/1‑148 (MQ=255)
                                                                                                                                             cGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACTTCGcc       <  6:2549/148‑1 (MQ=255)
                                                                                                                                              gACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACTTCGCCg      >  7:332701/1‑148 (MQ=255)
                                                                                                                                                ccGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACTTCGCCgtg    >  2:282966/1‑148 (MQ=255)
                                                                                                                                                  ggAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACTTCGCCGTGtt  >  4:194517/1‑148 (MQ=255)
                                                                                                                                                  ggAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACTTCGCCGTGtt  <  3:18875/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATGAAAGAGGAGCACTACCGCGCAGATACGCTGGATATCACCTTTGACGTCACTAAAACCACGCTCGAAGCGACAGTCAAAGAGCTGTGCGACAAAGCCGAAAAAATGGTACGTAGCGGCACCGTGCTGCTGGTGCTCTCCGACCGGAATATCGCTAAAGATCGCCTGCCGGTTCCAGCCCCGATGGCGGTTGGCGCGATCCAGACCCGTCTGGTCGATCAAAGCCTGCGTTGCGATGCCAACATCATCGTCGAAACCGCCAGCGCCCGCGATCCGCACCACTTCGCCGTGTT  >  NC_000913/3356425‑3356718

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: