Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,187,092 G→C 7.8% intergenic (+22/+27) pepT → / ← roxA peptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,187,0920GC 7.8% 301.8 / 13.2 90intergenic (+22/+27)pepT/roxApeptidase T/50S ribosomal protein L16 arginine hydroxylase; 2‑oxoglutarate oxygenase
Reads supporting (aligned to +/- strand):  ref base G (42/41);  new base C (5/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

CGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGTGCATC  >  NC_000913/1186943‑1187241
                                                                                                                                                     |                                                                                                                                                     
cggtcccGCACCTGTTCACTGGCGGTTCCAACTATCCTGGTACGCATGAGTTTGTGACTGTGGCAGGTATGGTCAAAGCGGGGCAGGTGCTCGTCCGTAATGCCGAGATAACGGCGCAACGGAAGTAAGCGAAAATGGATGCGGCAtgtg                                                                                                                                                       <  2:668061/146‑1 (MQ=255)
    cccGAACGTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGGGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                                   <  1:651530/150‑1 (MQ=255)
     ccGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGcc                                                                                                                                                  <  1:1140208/150‑1 (MQ=255)
      cGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCg                                                                                                                                                 >  1:465501/1‑150 (MQ=255)
        aaCCTGTTCACTGGCGGCTACAACTATCATGGTAAGCATGAGTGTGTGACTCTGGAAGGTATGGAAAATGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCa                                                                                                                                               <  2:189968/150‑1 (MQ=255)
          ccTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                                                                             <  2:596777/150‑1 (MQ=255)
           cTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATcc                                                                                                                                            <  1:1590492/150‑1 (MQ=255)
           cTGTTCACTGGCGGTTACAACTATCATGGGAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATcc                                                                                                                                            <  2:261681/150‑1 (MQ=255)
            tGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCg                                                                                                                                           >  1:343813/1‑150 (MQ=255)
              ttCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGc                                                                                                                                         <  1:405077/150‑1 (MQ=255)
                 aCTGGCGGTTACAACTATCATGGTAAGCATGATTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTa                                                                                                                                      <  1:967795/150‑1 (MQ=255)
                 aCTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAACGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTa                                                                                                                                      <  1:1122318/150‑1 (MQ=255)
                          tACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAAc                                                                                                                             <  1:650635/150‑1 (MQ=255)
                              aCTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATCGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGAAACGGAAGTAAGCGAAAAGGGATGCGGAATGTGATGCCGCATCCGGCTTAAATCCAAACTTAc                                                                                                                         >  1:327052/1‑150 (MQ=255)
                                tATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGTGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAccc                                                                                                                       >  1:534772/1‑150 (MQ=255)
                                tATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGTGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAccc                                                                                                                       >  1:638796/1‑150 (MQ=255)
                                 aTCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCt                                                                                                                      <  2:929736/150‑1 (MQ=255)
                                     tGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCga                                                                                                                  >  1:1256050/1‑150 (MQ=255)
                                     tGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCga                                                                                                                  >  1:409115/1‑150 (MQ=255)
                                      ggTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaa                                                                                                                 >  2:192174/1‑150 (MQ=255)
                                        tAAGCATGAGTTTGTGACTCTGGAACGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaaga                                                                                                               >  1:959106/1‑150 (MQ=255)
                                           gCATGAGTTTGTGACTCTGGAAGGTATTGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                                                                            <  2:301436/150‑1 (MQ=255)
                                           gCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                                                                            <  1:888680/150‑1 (MQ=255)
                                           gCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                                                                            >  2:4309/1‑150 (MQ=255)
                                           gCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                                                                            <  1:623392/150‑1 (MQ=255)
                                            cATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCAGGCTTAAATCCAAACTTACCCTTCGAAGAACCa                                                                                                           <  1:1259456/150‑1 (MQ=255)
                                               gAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATa                                                                                                        <  2:899992/150‑1 (MQ=255)
                                               gAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATa                                                                                                        >  1:992042/1‑150 (MQ=255)
                                                aGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAc                                                                                                       >  2:650913/1‑150 (MQ=255)
                                                 gTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAcc                                                                                                      <  2:375230/150‑1 (MQ=255)
                                                 gTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAcc                                                                                                      >  2:29697/1‑150 (MQ=255)
                                                  tttGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAccc                                                                                                     <  1:1297847/150‑1 (MQ=255)
                                                     gtgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTACGAAGAACCAATACCCGCt                                                                                                  >  2:1674413/1‑150 (MQ=255)
                                                      tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTa                                                                                                 >  1:634789/1‑150 (MQ=255)
                                                         ctctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTg                                                                                              <  2:409115/150‑1 (MQ=255)
                                                         ctctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTg                                                                                              <  2:1256050/150‑1 (MQ=255)
                                                         ctctGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTg                                                                                              <  1:1472418/150‑1 (MQ=255)
                                                         ctctGGAAGGAATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGATAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTg                                                                                              <  1:1303638/150‑1 (MQ=255)
                                                             ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                                                                          >  2:8296/1‑150 (MQ=255)
                                                             ggAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                                                                          >  1:855795/1‑150 (MQ=255)
                                                             ggAAGGTAGGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCa                                                                                          <  1:119924/150‑1 (MQ=255)
                                                               aaGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAgc                                                                                        >  1:89915/1‑150 (MQ=255)
                                                                 ggTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAgcgc                                                                                      >  2:35996/1‑150 (MQ=255)
                                                                 ggTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAgcgc                                                                                      >  2:107528/1‑150 (MQ=255)
                                                                        aaaaaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGc                                                                               >  2:365316/1‑150 (MQ=255)
                                                                                gTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGTGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                                                                       <  2:534772/150‑1 (MQ=255)
                                                                                gTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGTGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                                                                       <  2:638796/150‑1 (MQ=255)
                                                                                gTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAg                                                                       <  2:465501/150‑1 (MQ=255)
                                                                                  gCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGa                                                                     <  2:343813/150‑1 (MQ=255)
                                                                                    aGGTGATCGGCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAAt                                                                   <  1:339572/150‑1 (MQ=255)
                                                                                     tgTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATg                                                                  >  2:1489620/2‑150 (MQ=255)
                                                                                       tGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGAc                                                                <  2:623730/150‑1 (MQ=255)
                                                                                        gATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACg                                                               >  2:613050/1‑150 (MQ=255)
                                                                                           cGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGAt                                                            >  2:891055/1‑150 (MQ=255)
                                                                                            gTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATc                                                           >  1:1644652/1‑150 (MQ=255)
                                                                                                  aTTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAg                                                     <  1:1033809/150‑1 (MQ=255)
                                                                                                  aTTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAg                                                     <  1:192174/150‑1 (MQ=255)
                                                                                                   aatccGAGATACCGGCGCACCGGATGTAAGCGAAACGGGCGGCGGCATGTGAAGCCGCATCCGGCTTAAATCCAAACTTACCGTTCGACGAACCAAAACCCTCTATTTACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgc                                                    <  1:1131804/147‑1 (MQ=255)
                                                                                                    tGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgcg                                                   <  1:1306256/150‑1 (MQ=255)
                                                                                                    tGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgcg                                                   >  2:456453/1‑150 (MQ=255)
                                                                                                            tAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCa                                           >  2:898689/1‑150 (MQ=255)
                                                                                                                  cgcAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtc                                     >  2:185444/1‑150 (MQ=255)
                                                                                                                    cAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCCGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtctc                                   >  2:1621336/1‑150 (MQ=255)
                                                                                                                    cAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATtctc                                   >  1:983818/1‑150 (MQ=255)
                                                                                                                     aaCGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTcc                                  <  1:650913/150‑1 (MQ=255)
                                                                                                                      aCGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCg                                 >  1:1162331/1‑150 (MQ=255)
                                                                                                                      aCGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCg                                 >  2:1276913/1‑150 (MQ=255)
                                                                                                                         gAAGTAAGCGAAAAGGGATGTGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAg                              <  1:1582882/150‑1 (MQ=255)
                                                                                                                            gTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCa                           >  1:539499/1‑150 (MQ=255)
                                                                                                                              aaGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGGCCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCCTCGCCAAAATTATCCGCAGTCAgc                         >  2:176878/1‑150 (MQ=255)
                                                                                                                              aaGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGCCCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAATATTCTCCGCCGTCAgc                         >  2:354967/1‑150 (MQ=255)
                                                                                                                              aaGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGCGGAATGACGGATCTTCCCGCGCATCGCCAAAATTCTCCGCAGTCAgc                         >  2:182734/1‑150 (MQ=255)
                                                                                                                                gCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAgcgc                       <  1:992329/150‑1 (MQ=255)
                                                                                                                                gCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAgcgc                       <  1:376213/150‑1 (MQ=255)
                                                                                                                                 ttaagccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCa                      >  2:206070/8‑150 (MQ=255)
                                                                                                                                 cGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCa                      >  1:1224514/1‑150 (MQ=255)
                                                                                                                                  gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCaa                     >  1:701156/1‑150 (MQ=255)
                                                                                                                                    aaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATg                   <  2:98593/150‑1 (MQ=255)
                                                                                                                                     gccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGt                  >  2:604727/4‑150 (MQ=255)
                                                                                                                                      aGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGtt                 <  1:8296/150‑1 (MQ=255)
                                                                                                                                       cggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAGTGACGGATCTTCCAGCGCATCGCCAAACTTCTCCGCAGTCAGCGCAAGGTTg                >  1:490091/2‑150 (MQ=255)
                                                                                                                                        ggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGc               >  1:963167/1‑150 (MQ=255)
                                                                                                                                         gATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCt              >  2:1697794/1‑150 (MQ=255)
                                                                                                                                           tGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTgg            <  2:963167/150‑1 (MQ=255)
                                                                                                                                             cGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCg          <  2:1625767/150‑1 (MQ=255)
                                                                                                                                              ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGa         >  1:1331011/1‑150 (MQ=255)
                                                                                                                                               gCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAg        <  2:992042/150‑1 (MQ=255)
                                                                                                                                               gCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAg        >  2:172584/1‑150 (MQ=255)
                                                                                                                                                cATGTGATGCCGCATCCGGCTTAAATCCAACCTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGCCGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGt       >  2:1064124/1‑150 (MQ=255)
                                                                                                                                                cATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGt       <  1:172584/150‑1 (MQ=255)
                                                                                                                                                cATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGt       >  1:1131496/1‑150 (MQ=255)
                                                                                                                                                     gATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGTGCATc  >  2:1358060/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAGTGCATC  >  NC_000913/1186943‑1187241

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: