Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,568,919 T→C 9.6% intergenic (‑96/+35) yddW ← / ← gadC liprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,568,9190TC 9.6% 344.0 / 22.2 126intergenic (‑96/+35)yddW/gadCliprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base T (55/58);  new base C (10/2);  total (66/60)
Fisher's exact test for biased strand distribution p-value = 3.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 129 total aligned reads displayed.

AGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTA  >  NC_000913/1568771‑1569068
                                                                                                                                                    |                                                                                                                                                     
  tATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATc                                                                                                                                                    <  2:205857/150‑1 (MQ=255)
   aTCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCg                                                                                                                                                   >  1:1691772/1‑150 (MQ=255)
    tCGCTGGTCTTCTTATCGTTAATTTCTTGTTTCGGGAGCAGATTTCCATATGTGCTGGTTTCCGGTAATCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGAAAGCACGAAAAAGGGAGCGATTATTTATCGc                                                                                                                                                  >  1:304773/1‑150 (MQ=255)
        tGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCTGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTccc                                                                                                                                              <  2:802846/150‑1 (MQ=255)
        tGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTccc                                                                                                                                              >  2:1007696/1‑150 (MQ=255)
        tGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTccc                                                                                                                                              <  1:744855/150‑1 (MQ=255)
           tcttctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTCCCttt                                                                                                                                           <  1:974808/150‑2 (MQ=255)
           tcttctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCACCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTg                                                                                                                                           >  2:500207/1‑150 (MQ=255)
             ttctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATTTGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTc                                                                                                                                         <  2:1463697/150‑1 (MQ=255)
              tctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAAATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCt                                                                                                                                        <  2:612752/150‑1 (MQ=255)
               ctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCtt                                                                                                                                       >  2:1536565/1‑150 (MQ=255)
                tAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTta                                                                                                                                      >  1:488595/1‑150 (MQ=255)
                 aaTCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTtat                                                                                                                                     >  2:1120885/1‑150 (MQ=255)
                    cGTTATTTTGTTGTTTCGGGAGCAGATATCCAAATGGGGTGGGTTCCGGGAAAAAGATGTGCGGGCATTGTCGTACTGATCCCCGTTAAGTCCATACGACACCAAGCCCGAAAAAGGGAGCGATGAAATATCGCTACCTTGTCTTATAAc                                                                                                                                  >  2:1371211/1‑150 (MQ=255)
                    cGTTAATTTCTTGTTTCGGTAGCAGATATCCATATGTGCTGGTTTCCGTTAATCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGTTGAATTATCGCTCCCTTGTCTTATAAc                                                                                                                                  <  2:967889/150‑1 (MQ=255)
                       tAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCAt                                                                                                                               >  1:1591880/1‑150 (MQ=255)
                           ttCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAg                                                                                                                           >  1:1115832/1‑150 (MQ=255)
                             cTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTATGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGAc                                                                                                                         >  1:273909/1‑150 (MQ=255)
                                   tCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGtt                                                                                                                   <  2:961768/150‑1 (MQ=255)
                                     gggAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCATACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAg                                                                                                                 >  1:227124/1‑150 (MQ=255)
                                      ggAGCAGATTTCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCAGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCTCTCCCTTGTCTTATAACCATTCAGACATGGTTAgt                                                                                                                <  1:246394/150‑1 (MQ=255)
                                       gAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTAATTCTCGTACTTATCCCCGTTAAGTCAATACGCCAGCAAGCACGAAAAAGGGAGCGATGAATTATCGATCCCTTGTCTTATAACCATTCAGACATGGTTAgtg                                                                                                               >  1:1083383/1‑150 (MQ=255)
                                        aGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAgtgt                                                                                                              >  1:1499581/1‑150 (MQ=255)
                                         gCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAATGGGAGCGATGAATTATCGCTCCCTTGTCTTATAAGCATTCAGACATGGTTAGTGtt                                                                                                             <  2:1375478/150‑1 (MQ=255)
                                          cAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGttt                                                                                                            >  1:451678/1‑150 (MQ=255)
                                            gATATCCATATGTGCTGGTTTCCGGTAAACATAGGTGGGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGAACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCt                                                                                                          <  2:177535/150‑1 (MQ=255)
                                              tatCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTGCTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTg                                                                                                        <  2:966320/150‑1 (MQ=255)
                                              tatCCATATGTGCTGGTATCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTg                                                                                                        <  2:1461137/150‑1 (MQ=255)
                                               atCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGt                                                                                                       >  2:1246443/1‑150 (MQ=255)
                                                tCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtc                                                                                                      <  1:1632508/150‑1 (MQ=255)
                                                 ccATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtca                                                                                                     >  1:1015923/1‑150 (MQ=255)
                                                    tatGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattc                                                                                                  <  2:273909/150‑1 (MQ=255)
                                                     atGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattca                                                                                                 <  2:1659833/150‑1 (MQ=255)
                                                     atGAGCTGGTTTCCTGTAAACAGATGTTCGCTGATTCTCGTACTTATCCCCTTTAAGTCAATACGACAGCAATCATGAAATAGGGAGCGATGAATTATCGCTCCCTTGTCTTATATCCATTCAGACATGTTTAGTGTTTCTTGtcattca                                                                                                 <  1:130991/150‑1 (MQ=255)
                                                        tgCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatca                                                                                              >  1:505747/1‑150 (MQ=255)
                                                          cTGGTTTCCTGTAAACAGATGTGAGCTAATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGATCCATTGTCATATAACCATTCAGACATTGTTAGTCTTTCTTGTCATTCATcaca                                                                                            >  2:174981/1‑150 (MQ=255)
                                                           tGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCACCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACaa                                                                                           >  1:1052864/1‑150 (MQ=255)
                                                            ggTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAat                                                                                          >  2:78162/1‑150 (MQ=255)
                                                             gTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAata                                                                                         >  1:1622358/1‑150 (MQ=255)
                                                             gTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAata                                                                                         <  1:1532007/150‑1 (MQ=255)
                                                                 ccGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgt                                                                                     >  2:1318336/1‑150 (MQ=255)
                                                                 ccGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgt                                                                                     >  2:691483/1‑150 (MQ=255)
                                                                  cGGTAAACAGATGTGCGCTCGTTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtg                                                                                    >  2:1294831/1‑150 (MQ=255)
                                                                       aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggtg                                                                               <  2:339381/150‑1 (MQ=255)
                                                                       aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggtg                                                                               <  1:1365239/150‑1 (MQ=255)
                                                                        aCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGa                                                                              <  2:1376559/150‑1 (MQ=255)
                                                                           gATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACg                                                                           <  2:1123922/150‑1 (MQ=255)
                                                                             tgtgCTCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTg                                                                         <  2:1659680/150‑1 (MQ=255)
                                                                             tgtgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTg                                                                         >  1:485401/1‑150 (MQ=255)
                                                                             tgtgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCAATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTg                                                                         <  2:1681199/150‑1 (MQ=255)
                                                                               tgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGTACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCa                                                                       <  1:379610/150‑1 (MQ=255)
                                                                                   cTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAcgcg                                                                   >  2:1132739/1‑150 (MQ=255)
                                                                                     cATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCggg                                                                 >  1:1025497/1‑150 (MQ=255)
                                                                                       ttCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTg                                                               >  2:340618/1‑150 (MQ=255)
                                                                                              ccTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAACGGGAGCGATGAATTATCGCTCCCTTGTCTTATACCCATTCAGACAGGGTTACTGTTTCTTGTCATTCATAACAATATAGTGTGGTGAACGTGCACGCGGGGGCAGgaag                                                        >  2:79354/2‑150 (MQ=255)
                                                                                                 tATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGaagaag                                                     <  1:1793132/150‑1 (MQ=255)
                                                                                                  aTCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCCTTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGt                                                    <  1:1120885/150‑1 (MQ=255)
                                                                                                  aTCACCGTTAAGTCATTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGt                                                    <  2:88599/150‑1 (MQ=255)
                                                                                                   tCCCCGTTAAGTGAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTg                                                   <  2:829359/150‑1 (MQ=255)
                                                                                                     cccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAAACATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGaagaagagcc                                                 >  1:627619/1‑146 (MQ=255)
                                                                                                       cGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCt                                               <  2:505747/150‑1 (MQ=255)
                                                                                                          tAAGTCAATACGACAACAAGCACGAAAAAGGGAGCGATGATTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCtttt                                            <  1:413742/150‑1 (MQ=255)
                                                                                                              tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTg                                        <  1:485836/150‑1 (MQ=255)
                                                                                                                aaTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCg                                      >  1:1208070/1‑150 (MQ=255)
                                                                                                                aaTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCg                                      <  2:663142/150‑1 (MQ=255)
                                                                                                                   aCGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTAGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTc                                   <  1:1145490/150‑1 (MQ=255)
                                                                                                                    cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCt                                  <  1:1633213/150‑1 (MQ=255)
                                                                                                                    cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTAATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCATGAAGAAGTGACCTTTTGGTGCGTTCt                                  <  2:451678/150‑1 (MQ=255)
                                                                                                                    cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGCCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCt                                  >  1:1212908/1‑150 (MQ=255)
                                                                                                                     gACATCAATCACGAAAAAGGTATCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGGGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTg                                 <  1:906626/150‑1 (MQ=255)
                                                                                                                       cAGCAAGCACGAAAAAGGGATCGATTAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGAc                               <  1:258057/150‑1 (MQ=255)
                                                                                                                         gcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTg                             >  1:35501/1‑150 (MQ=255)
                                                                                                                         gcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATAACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTg                             >  2:1444385/1‑150 (MQ=255)
                                                                                                                            agcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGGTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTg                          <  2:480897/150‑1 (MQ=255)
                                                                                                                              caCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                        <  1:1253466/150‑1 (MQ=255)
                                                                                                                               aCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAtt                       <  2:580315/150‑1 (MQ=255)
                                                                                                                                cGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATTACCATTCAGACATGGTTTGTGTTTCTTGTCATACATCACAATATAGTGTGGAGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg                      >  1:783449/1‑150 (MQ=255)
                                                                                                                                cGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg                      <  1:1576953/150‑1 (MQ=255)
                                                                                                                                cGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg                      <  2:485401/150‑1 (MQ=255)
                                                                                                                                   acaaGGGAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGcg                   >  2:786098/3‑149 (MQ=255)
                                                                                                                                    caaGGGAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCATACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAAGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTc                  >  1:725802/2‑150 (MQ=255)
                                                                                                                                    caaGGGAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTc                  >  1:968303/2‑150 (MQ=255)
                                                                                                                                      aGGGAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                >  1:922212/1‑150 (MQ=255)
                                                                                                                                       gggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag               >  1:305914/1‑150 (MQ=255)
                                                                                                                                       gggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag               >  1:610740/1‑150 (MQ=255)
                                                                                                                                       gggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag               <  1:1441377/150‑1 (MQ=255)
                                                                                                                                       gggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag               >  1:1030967/1‑150 (MQ=255)
                                                                                                                                        ggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga              <  2:1052864/150‑1 (MQ=255)
                                                                                                                                        ggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga              >  2:1337982/1‑150 (MQ=255)
                                                                                                                                        ggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga              >  1:771695/1‑150 (MQ=255)
                                                                                                                                        ggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTCGGTGCGTTCTGACTGTTGATTGGCTCCaga              >  2:543823/1‑150 (MQ=255)
                                                                                                                                             gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGAc         >  1:1728516/1‑150 (MQ=255)
                                                                                                                                                taaTTATCGAGCCCTTGTCTTATAACCAGTCAGACTTGGGTAGTGTTTCTTGGCATTCAGCACAATATAGTGTGGTGAACGTGCTCGCTGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc      <  1:174981/149‑1 (MQ=255)
                                                                                                                                                gAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc      >  1:1750693/1‑150 (MQ=255)
                                                                                                                                                 aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg     >  1:1152296/1‑150 (MQ=255)
                                                                                                                                                 aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg     >  2:847831/1‑150 (MQ=255)
                                                                                                                                                  aTTATCGCTCCCTTGTCTTATAACCATGCAGCGATGGTTTGTGTTTCTTGTCAGTCATCACAATATAGTGTGTTGAACGTGCACGGGGGGGCAGGATGAAGTGACCATTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg    >  2:275252/1‑150 (MQ=255)
                                                                                                                                                   ttATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGt   <  1:1536565/150‑1 (MQ=255)
                                                                                                                                                    tATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTa  <  1:21738/150‑1 (MQ=255)
                                                                                                                                                    tATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTa  <  2:1025497/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTA  >  NC_000913/1568771‑1569068

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