Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,302,843 T→C 8.5% intergenic (+30/‑56) 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,8430TC 8.5% 266.1 / 12.3 94intergenic (+30/‑56)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 T (48/38);  new base C (5/3);  total (53/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

TCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCG  >  NC_000913/1302699‑1302988
                                                                                                                                                |                                                                                                                                                 
tCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACg                                                                                                                                              >  2:754539/1‑150 (MQ=255)
tCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGACACTCCTGCCCCACg                                                                                                                                              <  1:1194734/150‑1 (MQ=255)
 ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACGGTGGGTTGGGGGCTATACCGTCAAAGTTCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATAGGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                             <  2:783013/150‑1 (MQ=255)
 ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                             >  2:1136107/1‑150 (MQ=255)
 ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                             <  2:223684/150‑1 (MQ=255)
 ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAAAACCTATACCCGTAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                             <  2:1303544/150‑1 (MQ=255)
       tattaCTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTg                                                                                                                                       <  1:915024/150‑1 (MQ=255)
        attaCTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGAGCAGGACACTCCTGCCCCAc                                                                                                                                               >  1:1078271/1‑141 (MQ=255)
        attaCTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGAGCAGGACACTCCTGCCCCAc                                                                                                                                               <  2:1078271/141‑1 (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                                                                                                                                                 <  1:1454321/119‑1 (MQ=255)
                            cTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTcc                                                                                                                                                 >  2:1454321/1‑119 (MQ=255)
                             tGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAgg                                                                                                                 <  1:1559666/150‑1 (MQ=255)
                             tGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAgg                                                                                                                 >  1:944957/1‑150 (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:420776/1‑150 (MQ=255)
                                                 ggCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGc                                                                                             <  1:175083/150‑1 (MQ=255)
                                                     ataCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATg                                                                                         >  2:750652/1‑150 (MQ=255)
                                                     ataCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATg                                                                                         >  2:756381/1‑150 (MQ=255)
                                                      taCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGt                                                                                        >  1:1358686/1‑150 (MQ=255)
                                                      taCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGt                                                                                        >  1:1006914/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:650494/1‑150 (MQ=255)
                                                                    cGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTa                                                                          >  2:1609229/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                                               >  2:1139163/1‑150 (MQ=255)
                                                                                               tGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACg                                               <  1:1269304/150‑1 (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:381859/1‑150 (MQ=255)
                                                                                                           aGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCtt                                   >  2:398198/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                          <  2:1358686/150‑1 (MQ=255)
                                                                                                                    ggCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTAtt                          >  1:216888/1‑150 (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         >  1:1364391/1‑146 (MQ=255)
                                                                                                                                         gtCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATgcg         <  2:1364391/146‑1 (MQ=255)
                                                                                                                                          tCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTcg    >  1:1190933/1‑150 (MQ=255)
                                                                                                                                            cTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTcgcg  >  2:286449/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCG  >  NC_000913/1302699‑1302988

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: