Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,302,836 G→C 8.7% intergenic (+23/‑63) oppA → / → oppB oligopeptide ABC transporter periplasmic binding protein/murein tripeptide ABC transporter/oligopeptide ABC transporter inner membrane subunit OppB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,302,8360GC 8.7% 288.6 / 15.8 92intergenic (+23/‑63)oppA/oppBoligopeptide ABC transporter periplasmic binding protein/murein tripeptide ABC transporter/oligopeptide ABC transporter inner membrane subunit OppB
Reads supporting (aligned to +/- strand):  ref base G (47/37);  new base C (5/3);  total (52/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 = 5.50e-01

CGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCG  >  NC_000913/1302689‑1302981
                                                                                                                                                   |                                                                                                                                                 
cgcccATTGTTCCTGTTTATTACTACGAGAATGCGCGTCTGCTGAAACCGTGGGTTGGTGGCTATACAGGCAAAGATCAGCTGGATAAAACCTATAAAAGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGCGGAAGGAGTGTc                                                                                                                                                 >  2:510647/4‑150 (MQ=255)
cGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTc                                                                                                                                                 >  2:1247107/1‑150 (MQ=255)
          tCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACg                                                                                                                                       >  2:754539/1‑150 (MQ=255)
          tCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGACACTCCTGCCCCACg                                                                                                                                       <  1:1194734/150‑1 (MQ=255)
           ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACGGTGGGTTGGGGGCTATACCGTCAAAGTTCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATAGGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                      <  2:783013/150‑1 (MQ=255)
           ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                      <  2:223684/150‑1 (MQ=255)
           ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                      >  2:1136107/1‑150 (MQ=255)
           ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAAAACCTATACCCGTAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                      <  2:1303544/150‑1 (MQ=255)
                 tattaCTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTg                                                                                                                                <  1:915024/150‑1 (MQ=255)
                  attaCTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGAGCAGGACACTCCTGCCCCAc                                                                                                                                        <  2:1078271/141‑1 (MQ=255)
                  attaCTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGAGCAGGACACTCCTGCCCCAc                                                                                                                                        >  1:1078271/1‑141 (MQ=255)
                            gcaTGCGCGGCGGGCGACACCGTGGGGCGGCGGCCATACCGGCAAAGCTCCGCCGGATAATACCTCTACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGGGGGGCAGGAGGGTCCTGCTCCACGGTGTCTGATTTTTATCGc                                                                                                                     <  1:1543667/148‑1 (MQ=255)
                            gAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGc                                                                                                                     <  2:526097/150‑1 (MQ=255)
                             aaTGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCa                                                                                                                    >  2:380056/1‑150 (MQ=255)
                              aTGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCAt                                                                                                                   >  1:824371/1‑150 (MQ=255)
                                 cgcgTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTAc                                                                                                                >  1:743446/1‑150 (MQ=255)
                                      cTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTcc                                                                                                                                          >  2:1454321/1‑119 (MQ=255)
                                      cTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTcc                                                                                                                                          <  1:1454321/119‑1 (MQ=255)
                                       tGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAgg                                                                                                          >  1:944957/1‑150 (MQ=255)
                                       tGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAgg                                                                                                          <  1:1559666/150‑1 (MQ=255)
                                         gTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGca                                                                                                        >  1:1502951/1‑150 (MQ=255)
                                         gTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGca                                                                                                        >  1:255750/1‑150 (MQ=255)
                                               ccGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGcc                                                                                                  <  1:1421676/150‑1 (MQ=255)
                                                   gggTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGaa                                                                                              >  2:1186443/1‑150 (MQ=255)
                                                    ggTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAg                                                                                             <  1:23564/150‑1 (MQ=255)
                                                    ggTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAg                                                                                             >  1:1075778/1‑150 (MQ=255)
                                                     gTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAgg                                                                                            <  1:754539/150‑1 (MQ=255)
                                                      ttGGTGTCTATAGCGTCAAAGATCCGCTTGATAATACCTATAGCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAATGCACAGGCCAGAAGGt                                                                                           <  2:248651/150‑1 (MQ=255)
                                                         gtggCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCTCAGGCCAGAAGGTAgg                                                                                        <  1:1217810/150‑1 (MQ=255)
                                                           ggCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGc                                                                                      <  1:175083/150‑1 (MQ=255)
                                                           ggCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGc                                                                                      >  1:420776/1‑150 (MQ=255)
                                                               ataCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATg                                                                                  >  2:750652/1‑150 (MQ=255)
                                                               ataCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATg                                                                                  >  2:756381/1‑150 (MQ=255)
                                                                taCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGt                                                                                 >  1:1006914/1‑150 (MQ=255)
                                                                taCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGt                                                                                 >  1:1358686/1‑150 (MQ=255)
                                                                  ccGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTa                                                                               >  2:985297/1‑150 (MQ=255)
                                                                  ccGGCAAAGATCCGCTGGATAATAACTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTa                                                                               >  2:911177/1‑150 (MQ=255)
                                                                   cGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGTCCAGAAGGTAGGGCAATGTTaa                                                                              <  2:172990/150‑1 (MQ=255)
                                                                    ggCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTaaa                                                                             >  1:837600/1‑150 (MQ=255)
                                                                     tcAAAGATCCGCGTGATAATACCTATTGCCGGAATATTTACATTGTGACGCACTATTGGCAATTCGTGGGGCAGGAGTGGCCTGCTCCACGGTGTCTGATTTTTATCGCATTAATGAAGGCACAGGCCAGACGGTAGGGCAATGTTaaaa                                                                            <  2:715544/149‑1 (MQ=255)
                                                                      cAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAgg                                                                                        >  2:902177/1‑137 (MQ=255)
                                                                      cAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAgg                                                                                        <  1:902177/137‑1 (MQ=255)
                                                                      cAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAAt                                                                           >  2:1519616/1‑150 (MQ=255)
                                                                          gATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTAt                                                                       <  2:940429/150‑1 (MQ=255)
                                                                          gATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTAt                                                                       >  1:1486177/1‑150 (MQ=255)
                                                                             ccGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCa                                                                                                 >  1:1272935/1‑121 (MQ=255)
                                                                             ccGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCa                                                                                                 <  2:1272935/121‑1 (MQ=255)
                                                                             ccGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCt                                                                    <  2:346105/150‑1 (MQ=255)
                                                                              cGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTa                                                                   >  2:1609229/1‑150 (MQ=255)
                                                                              cGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTa                                                                   >  2:650494/1‑150 (MQ=255)
                                                                               gCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTAc                                                                  >  2:1230757/1‑150 (MQ=255)
                                                                                    ataataCCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGc                                                             <  1:2904/150‑1 (MQ=255)
                                                                                         aCCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCt                                                        <  1:1020109/150‑1 (MQ=255)
                                                                                          catataACCGTAATATGTACATTGAGAATCACTAATTGAAAAACGAGGGGCAGGAGTGAACTGATCCAAGATGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTg                                                       >  2:1315291/3‑150 (MQ=255)
                                                                                            tataCCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGa                                                     >  1:805795/1‑150 (MQ=255)
                                                                                               aCCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGc                                                  <  1:1136107/150‑1 (MQ=255)
                                                                                                cccGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCg                                                 >  1:858884/1‑150 (MQ=255)
                                                                                                 ccGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCTCAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGa                                                <  1:1247107/150‑1 (MQ=255)
                                                                                                         tGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACg                                        <  1:1269304/150‑1 (MQ=255)
                                                                                                         tGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACg                                        >  2:1139163/1‑150 (MQ=255)
                                                                                                          gTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGc                                       <  1:1609229/150‑1 (MQ=255)
                                                                                                               ttGTGAAGCACTAATGGCAATACGTGGGGAAGGAGTGTCCTGATCAACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGATGTCTGGAAGCGATTCAGACGCtattt                                  >  2:178161/1‑150 (MQ=255)
                                                                                                                tgtgAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTACTCTACGTAGCTGTCTGGAAGCGATTCCGACGCtattta                                 <  2:783511/150‑1 (MQ=255)
                                                                                                                 gtgAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCtatttat                                >  2:1582145/1‑150 (MQ=255)
                                                                                                                     aGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAg                                                   >  2:1198715/1‑127 (MQ=255)
                                                                                                                     aGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAg                                                   <  1:1198715/127‑1 (MQ=255)
                                                                                                                     aGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCtt                            >  2:398198/1‑150 (MQ=255)
                                                                                                                     aGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCtt                            >  2:381859/1‑150 (MQ=255)
                                                                                                                        aCTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCttatt                         >  2:1184840/1‑150 (MQ=255)
                                                                                                                        aCTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCttatt                         <  1:650494/150‑1 (MQ=255)
                                                                                                                         cTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCttatta                        >  2:533248/1‑150 (MQ=255)
                                                                                                                         cTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTCATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCttatta                        >  1:1518462/1‑150 (MQ=255)
                                                                                                                          tAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTAc                       <  2:342923/150‑1 (MQ=255)
                                                                                                                           aaTGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACt                      >  2:1400992/1‑150 (MQ=255)
                                                                                                                              ggCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTAtt                   >  1:216888/1‑150 (MQ=255)
                                                                                                                              ggCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTAtt                   <  2:1358686/150‑1 (MQ=255)
                                                                                                                               gCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTAttt                  <  1:256312/150‑1 (MQ=255)
                                                                                                                               agacaccGTGGAGCAGGACACTCCTGCCCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCtt                            <  2:475768/134‑1 (MQ=255)
                                                                                                                               agacaccGTGGAGCAGGACACTCCTGCCCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCtt                            >  1:475768/7‑140 (MQ=255)
                                                                                                                                gacaccGTGGAGCAGGACACTCCTGCCCAACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTGCGACGCTATTTATTCTTATTACTATTTc                 >  1:1093074/6‑150 (MQ=255)
                                                                                                                                cAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTc                 >  2:520146/1‑150 (MQ=255)
                                                                                                                                 aaTACGTGGGGCAGGAGTGTCCTGCTCCACTGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCg                <  1:1101027/150‑1 (MQ=255)
                                                                                                                                 aaTACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTTGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCg                >  2:35611/1‑150 (MQ=255)
                                                                                                                                  caccGTGGAGCAGGACACTCCTGCCCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGt               >  1:656064/4‑150 (MQ=255)
                                                                                                                                  aTACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGt               >  1:270945/1‑150 (MQ=255)
                                                                                                                                   accGTGGAGCAGGACACTCCTGCCCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAAAGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGtt              >  2:1596156/3‑150 (MQ=255)
                                                                                                                                     cGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCt            >  2:760539/1‑150 (MQ=255)
                                                                                                                                        ggggCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTa         >  1:145835/1‑150 (MQ=255)
                                                                                                                                           gCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTatga      <  1:852539/150‑1 (MQ=255)
                                                                                                                                            cAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTatgat     <  1:756381/150‑1 (MQ=255)
                                                                                                                                            cAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTatgat     <  1:750652/150‑1 (MQ=255)
                                                                                                                                             aGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTatgatg    <  2:420776/150‑1 (MQ=255)
                                                                                                                                                   gtCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATgcg  <  2:1364391/146‑1 (MQ=255)
                                                                                                                                                   gtCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATgcg  >  1:1364391/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCG  >  NC_000913/1302689‑1302981

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: