Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,568,917 A→T 6.6% 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,9170AT 6.6% 318.2 / 9.4 116intergenic (‑94/+37)yddW/gadCliprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base A (58/49);  new base T (2/6);  total (60/56)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

Only 100 of 118 total aligned reads displayed.

ACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGG  >  NC_000913/1568768‑1569066
                                                                                                                                                     |                                                                                                                                                     
   aGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTa                                                                                                                                                    >  2:277071/1‑150 (MQ=255)
    gTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTAt                                                                                                                                                   >  2:858630/1‑150 (MQ=255)
     tATCTCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGAGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATc                                                                                                                                                  >  2:714950/1‑150 (MQ=255)
         gCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTc                                                                                                                                              <  2:1126580/150‑1 (MQ=255)
          cTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTcc                                                                                                                                             >  2:1072160/1‑150 (MQ=255)
          cTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTcc                                                                                                                                             >  2:1711009/1‑150 (MQ=255)
           ttctcttctAATAGTTAATTTCTTGTTTCGGGAGCAGATCTCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTAACCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTccc                                                                                                                                            <  2:292226/147‑1 (MQ=255)
           tGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTccc                                                                                                                                            <  2:1200226/150‑1 (MQ=255)
           tGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTccc                                                                                                                                            <  2:1189065/150‑1 (MQ=255)
            ggTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCt                                                                                                                                           >  1:1556175/1‑150 (MQ=255)
            ggTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTATCCGGTAAACAGATGTGCGCTCATTCTCGTACTTAACCCCGTTAAGTCAATACGACAGCAAACACGAAAAAGCGAGCGATGAATTATCGCTCCCt                                                                                                                                           >  2:1719112/1‑150 (MQ=255)
            ggTCTTCTAATCGTTAAATTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCt                                                                                                                                           <  1:279788/150‑1 (MQ=255)
              tcttctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTg                                                                                                                                         <  1:911425/150‑1 (MQ=255)
               cttctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGt                                                                                                                                        >  1:1587053/1‑150 (MQ=255)
                ttctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTc                                                                                                                                       <  1:719796/150‑1 (MQ=255)
                ttctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTCCCttttt                                                                                                                                       <  2:857838/150‑4 (MQ=255)
                ttctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTCATCGCTCCCttttt                                                                                                                                       <  2:1465292/150‑4 (MQ=255)
                  ctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCtt                                                                                                                                     >  2:576523/1‑150 (MQ=255)
                  ctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCATTGTCtt                                                                                                                                     >  2:1631162/1‑150 (MQ=255)
                     aTCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTGTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTtata                                                                                                                                  <  1:840939/150‑1 (MQ=255)
                       cGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCTTTAAGTCAATACGACAGCTAGCACGAAAAAGGGAGCGATGTATTATCGCTCCCTTTTCTTATAAc                                                                                                                                >  1:332208/1‑150 (MQ=255)
                       cGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAAc                                                                                                                                <  1:396782/150‑1 (MQ=255)
                        gTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAAcc                                                                                                                               >  1:698886/1‑150 (MQ=255)
                         ttAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTATCTTACTTATACACGTTAATTAAATAAGACAGCAAGCACGAAAAAGGGAGCGATGAAGTATCGATACCTTGTCTTATAACCa                                                                                                                              >  1:81037/1‑150 (MQ=255)
                          tAATTTCTTGTTTCGGGAGCTGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTAGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCAtt                                                                                                                            <  2:653434/150‑1 (MQ=255)
                               tCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGa                                                                                                                        >  1:1155200/1‑150 (MQ=255)
                                     ttCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGAACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGt                                                                                                                  >  2:1376622/1‑150 (MQ=255)
                                         ggAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAgt                                                                                                              >  2:1189518/1‑150 (MQ=255)
                                            gCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTAATTCTCGTACTTATGCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCATTGTCTTATAACCATTCAGACATGGTTAGTGtt                                                                                                           >  1:285833/1‑150 (MQ=255)
                                             cAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGttt                                                                                                          >  2:612486/1‑150 (MQ=255)
                                                    ccATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtca                                                                                                   >  1:22353/1‑150 (MQ=255)
                                                     cATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATTCGACAGCAAGCCCGAAACAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcat                                                                                                  >  1:81878/1‑150 (MQ=255)
                                                          gtgCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatc                                                                                             <  2:1331232/150‑1 (MQ=255)
                                                               ggTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAat                                                                                        <  2:1363543/150‑1 (MQ=255)
                                                                 tttCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAatat                                                                                      <  1:967603/150‑1 (MQ=255)
                                                                   tCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAg                                                                                    <  1:733584/150‑1 (MQ=255)
                                                                   tCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAg                                                                                    <  2:1250281/150‑1 (MQ=255)
                                                                    ccGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgt                                                                                   >  2:2074/1‑150 (MQ=255)
                                                                    ccGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgt                                                                                   >  1:974833/1‑150 (MQ=255)
                                                                     cGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtg                                                                                  >  1:774329/1‑150 (MQ=255)
                                                                     cGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtg                                                                                  <  1:1427686/150‑1 (MQ=255)
                                                                      ggTAAACAGATGTGCGCTCATTCTCGTACTTTTCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtgt                                                                                 >  1:33582/1‑150 (MQ=255)
                                                                      ggTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtgt                                                                                 >  1:39232/1‑150 (MQ=255)
                                                                       gTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCAAAATATAGtgtg                                                                                >  1:1268390/1‑150 (MQ=255)
                                                                          aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggtg                                                                             >  1:18508/1‑150 (MQ=255)
                                                                          aaCAGATGTGCGCTCATTCTCGTACTTATCACCGTTAAGTCAATACGACAGCAAGCTCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTAAAACCAATCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggtg                                                                             <  2:339012/150‑1 (MQ=255)
                                                                            cAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAAAATAGTGTGGTGaa                                                                           >  1:1342832/1‑150 (MQ=255)
                                                                              gATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCATACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACg                                                                         >  2:1544289/1‑150 (MQ=255)
                                                                              gATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACg                                                                         <  1:859264/150‑1 (MQ=255)
                                                                              gATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATAAGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACg                                                                         <  1:137928/150‑1 (MQ=255)
                                                                                 gtgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGc                                                                      <  1:1575048/150‑1 (MQ=255)
                                                                                 gtgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGc                                                                      <  1:754798/150‑1 (MQ=255)
                                                                                  tgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATTGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCa                                                                     >  1:813664/1‑150 (MQ=255)
                                                                                            ctcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGGCTTATTACCATTGAGGCAGGGTTTGTGTTTCTTGTCGTTCATCACTATATAGTGTGGTGAACGTGCACGCGGGTGcc                                                           >  2:1156788/1‑149 (MQ=255)
                                                                                            atcGTACTTATCCCCGTTAAGTGAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCa                                                           <  2:814138/149‑1 (MQ=255)
                                                                                             tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAg                                                          <  1:1711009/150‑1 (MQ=255)
                                                                                             tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAg                                                          <  1:881183/150‑1 (MQ=255)
                                                                                              cGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAgg                                                         <  1:117162/150‑1 (MQ=255)
                                                                                              cGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATACAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAgg                                                         <  1:1072160/150‑1 (MQ=255)
                                                                                                   ttATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACCCGCGTGCacgttcac                                                    <  1:714950/150‑8 (MQ=255)
                                                                                                     aTCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGt                                                  >  2:196710/1‑150 (MQ=255)
                                                                                                     aTCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGt                                                  >  2:951825/1‑150 (MQ=255)
                                                                                                     aTCCCCGTTAAGTCAATACGAAAGCAAGCACGAAAAATGGAGCGATTAATTATCGCTCCCTTGTCTTTTAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGt                                                  <  2:1259133/150‑1 (MQ=255)
                                                                                                      tCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTATCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGTGTGCAGGAAGAAGTg                                                 >  1:1239078/1‑150 (MQ=255)
                                                                                                        cccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGAc                                               <  2:409244/150‑1 (MQ=255)
                                                                                                              aaGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTg                                         <  1:1601006/150‑1 (MQ=255)
                                                                                                                 tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTg                                      <  2:1556175/150‑1 (MQ=255)
                                                                                                                 tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTg                                      <  2:716744/150‑1 (MQ=255)
                                                                                                                 tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTg                                      >  2:572702/1‑150 (MQ=255)
                                                                                                                 tCAATAAGACAGCAAGCACGAACAAGGGAGCGATGACTTATCGCTCCCTTGTCTTATAACCATTCAGACTTGGATTGTGTTTCTTGTCATTCATCACAATATTGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTg                                      <  1:834480/150‑1 (MQ=255)
                                                                                                                  cAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGc                                     >  2:1011893/1‑150 (MQ=255)
                                                                                                                  cAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGc                                     >  1:876083/1‑150 (MQ=255)
                                                                                                                       cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCt                                <  2:39232/150‑1 (MQ=255)
                                                                                                                        gACAGCAAGCACGAAAAATGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTg                               <  2:658233/150‑1 (MQ=255)
                                                                                                                        gACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTg                               <  2:802909/150‑1 (MQ=255)
                                                                                                                         aCAGCAAGCACGAAAAAGGGAGCGATGAAATATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGTTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGa                              <  2:1185537/150‑1 (MQ=255)
                                                                                                                            gcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTg                           >  1:1372967/1‑150 (MQ=255)
                                                                                                                             caagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGt                          <  1:1011893/150‑1 (MQ=255)
                                                                                                                               agcaCGAAAAAGGGAGCGATGTATTATAGCTCAGTTGTCTTATAACCTTTCTGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTTGTGAACGTGCACGCGGGTGCAGGAAGAATTGACCTTTTGGTGCGTTCTGACTGTTg                        <  1:1376622/150‑1 (MQ=255)
                                                                                                                               agcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTg                        <  1:1031827/150‑1 (MQ=255)
                                                                                                                               agcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTg                        >  1:644321/1‑150 (MQ=255)
                                                                                                                                gcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGa                       >  2:1080565/1‑150 (MQ=255)
                                                                                                                                  aCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAtt                     <  2:1423722/150‑1 (MQ=255)
                                                                                                                                   cGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATACAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg                    <  2:321044/150‑1 (MQ=255)
                                                                                                                                    gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATACCCCTTCAGACATGGGTAGTGTTTCTTGTCCTTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTgg                   <  1:277071/150‑1 (MQ=255)
                                                                                                                                    gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAANCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTgg                   >  2:427678/1‑150 (MQ=255)
                                                                                                                                    gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTgg                   >  2:940554/1‑150 (MQ=255)
                                                                                                                                    gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTgg                   >  2:421258/1‑150 (MQ=255)
                                                                                                                                     aaaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAAAATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGc                  <  2:1342832/150‑1 (MQ=255)
                                                                                                                                          gggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGCCCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag             >  2:299226/1‑150 (MQ=255)
                                                                                                                                            gAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCagag           <  2:146087/150‑1 (MQ=255)
                                                                                                                                               cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTTGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGAGTTCTGACTGTTTATTGGCTCAAGAGTGa        >  2:205514/1‑150 (MQ=255)
                                                                                                                                               cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa        >  1:676219/1‑150 (MQ=255)
                                                                                                                                               cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa        >  2:707300/1‑150 (MQ=255)
                                                                                                                                                gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGAc       >  1:369994/1‑150 (MQ=255)
                                                                                                                                                gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCAAGAGTGAc       >  2:921434/1‑150 (MQ=255)
                                                                                                                                                  tGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGc     >  2:372500/1‑150 (MQ=255)
                                                                                                                                                   gAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc    <  2:698886/150‑1 (MQ=255)
                                                                                                                                                   gAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGGTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc    <  1:209594/150‑1 (MQ=255)
                                                                                                                                                     aTTATCGCTCCCTTGTCTTATAACCCATCAGACATGGTTAGTGTTTCTTGTCCTTCATCACAAAATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg  <  2:285833/150‑1 (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 < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: