Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,568,917 A→T 11.1% intergenic (‑94/+37) 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,9170AT11.1% 307.9 / 21.4 122intergenic (‑94/+37)yddW/gadCliprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base A (52/56);  new base T (11/3);  total (63/59)
Fisher's exact test for biased strand distribution p-value = 4.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

ACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGG  >  NC_000913/1568768‑1569066
                                                                                                                                                     |                                                                                                                                                     
aCTAGTATCGCTGGTCTTGTAATCGTTAATTTCATGTTTCGGGAGCAGATATCCATATGTGCTGGTATCCGGAAAACAGATGTGCTCTCATTATCGTACTTATCCCCGTTAAGTCAATACTACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                       <  1:260580/150‑1 (MQ=255)
   aGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTa                                                                                                                                                    >  1:1749419/1‑150 (MQ=255)
       tCGCTGGTCTTCTTATCGTTAATTTCTTGTTTCGGGAGCAGATTTCCATATGTGCTGGTTTCCGGTAATCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGAAAGCACGAAAAAGGGAGCGATTATTTATCGc                                                                                                                                                >  1:304773/1‑150 (MQ=255)
       tCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGc                                                                                                                                                <  1:1371739/150‑1 (MQ=255)
           tGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTccc                                                                                                                                            >  2:1007696/1‑150 (MQ=255)
           tGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTccc                                                                                                                                            <  1:744855/150‑1 (MQ=255)
           tGGTCTTCAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCt                                                                                                                                           >  2:485836/1‑150 (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:1123922/1‑150 (MQ=255)
                   tAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTta                                                                                                                                    <  1:179920/150‑1 (MQ=255)
                     aTCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTtata                                                                                                                                  <  2:1746495/150‑1 (MQ=255)
                       cGTTAATTTCTTGTTTCGGTAGCAGATATCCATATGTGCTGGTTTCCGTTAATCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGTTGAATTATCGCTCCCTTGTCTTATAAc                                                                                                                                <  2:967889/150‑1 (MQ=255)
                              ttCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAg                                                                                                                         >  1:1115832/1‑150 (MQ=255)
                                cTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTATGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGAc                                                                                                                       >  1:273909/1‑150 (MQ=255)
                                      tCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGtt                                                                                                                 <  2:961768/150‑1 (MQ=255)
                                       cGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTa                                                                                                                >  2:380836/1‑150 (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)
                                               gATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCt                                                                                                        >  2:1494661/1‑150 (MQ=255)
                                               gATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCt                                                                                                        <  2:609250/150‑1 (MQ=255)
                                                 tatCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTGCTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTg                                                                                                      <  2:966320/150‑1 (MQ=255)
                                                  atCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGt                                                                                                     >  2:1246443/1‑150 (MQ=255)
                                                   tCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtc                                                                                                    >  1:339381/1‑150 (MQ=255)
                                                   tCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtc                                                                                                    <  1:1632508/150‑1 (MQ=255)
                                                    ccATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtca                                                                                                   >  1:1015923/1‑150 (MQ=255)
                                                    ccATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtca                                                                                                   <  1:163027/150‑1 (MQ=255)
                                                       tatGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattc                                                                                                <  2:273909/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)
                                                                          aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggtg                                                                             <  1:1365239/150‑1 (MQ=255)
                                                                           aCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGa                                                                            <  2:1376559/150‑1 (MQ=255)
                                                                              gATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAGGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACg                                                                         <  1:121690/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)
                                                                                 gtgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGc                                                                      <  1:1561079/150‑1 (MQ=255)
                                                                                  tgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGTACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCa                                                                     <  1:379610/150‑1 (MQ=255)
                                                                                      cTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAcgcg                                                                 >  2:1132739/1‑150 (MQ=255)
                                                                                      cTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCCATGAATAATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTTGTGTTTCTTGTCATTCATCACAATATAGTGTGCTGAACGTGCAcgcg                                                                 >  2:120038/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)
                                                                                                    tATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTTAACGTGCACGCGGGTGCAGGaagaag                                                   <  1:500207/150‑1 (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)
                                                                                                                  cAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGc                                     >  2:269097/1‑150 (MQ=255)
                                                                                                                   aaTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCg                                    <  2:663142/150‑1 (MQ=255)
                                                                                                                   aaTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCg                                    >  1:1208070/1‑150 (MQ=255)
                                                                                                                      aCGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCACTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTc                                 >  2:827401/1‑150 (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)
                                                                                                                           agcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACt                            <  2:1418308/150‑1 (MQ=255)
                                                                                                                           agcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACt                            <  2:1104255/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                      >  2:773479/1‑150 (MQ=255)
                                                                                                                                   cGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATTACCATTCAGACATGGTTTGTGTTTCTTGTCATACATCACAATATAGTGTGGAGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg                    >  1:783449/1‑150 (MQ=255)
                                                                                                                                   cGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg                    <  1:1576953/150‑1 (MQ=255)
                                                                                                                                    agacaaGGGAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACCTGGTGAGTGTTTCTTGTAATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGAACTTTTGGTGCGTTCTGACTGTTGATTgg                   >  1:263996/5‑150 (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:610740/1‑150 (MQ=255)
                                                                                                                                          gggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag             >  1:1030967/1‑150 (MQ=255)
                                                                                                                                           ggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga            <  2:1052864/150‑1 (MQ=255)
                                                                                                                                           ggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga            >  1:771695/1‑150 (MQ=255)
                                                                                                                                           ggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga            >  2:1337982/1‑150 (MQ=255)
                                                                                                                                           ggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTCGGTGCGTTCTGACTGTTGATTGGCTCCaga            >  2:543823/1‑150 (MQ=255)
                                                                                                                                            gAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCagag           <  2:1366595/150‑1 (MQ=255)
                                                                                                                                                gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGAc       >  1:1728516/1‑150 (MQ=255)
                                                                                                                                                   taaTTATCGAGCCCTTGTCTTATAACCAGTCAGACTTGGGTAGTGTTTCTTGGCATTCAGCACAATATAGTGTGGTGAACGTGCTCGCTGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc    <  1:174981/149‑1 (MQ=255)
                                                                                                                                                    aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg   >  2:847831/1‑150 (MQ=255)
                                                                                                                                                    aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg   >  1:1152296/1‑150 (MQ=255)
                                                                                                                                                     aTTATCGCTCCCTTGTCTTATAACCATGCAGCGATGGTTTGTGTTTCTTGTCAGTCATCACAATATAGTGTGTTGAACGTGCACGGGGGGGCAGGATGAAGTGACCATTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg  >  2:275252/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGG  >  NC_000913/1568768‑1569066

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: