Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,568,915 G→A 9.8% intergenic (‑92/+39) 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,9150GA 9.8% 350.0 / 19.4 123intergenic (‑92/+39)yddW/gadCliprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base G (52/57);  new base A (10/2);  total (63/60)
Fisher's exact test for biased strand distribution p-value = 3.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

Only 100 of 126 total aligned reads displayed.

CAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCC  >  NC_000913/1568766‑1569064
                                                                                                                                                     |                                                                                                                                                     
cAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATg                                                                                                                                                       >  2:1532007/1‑150 (MQ=255)
cAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATg                                                                                                                                                       >  2:834288/1‑150 (MQ=255)
 aaCTACTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCACGTTAAGTAAAAACGACAGAAAGCACGAAAAAGGGAGCGATGa                                                                                                                                                      >  1:88599/1‑150 (MQ=255)
  aCTAGTATCGCTGGTCTTGTAATCGTTAATTTCATGTTTCGGGAGCAGATATCCATATGTGCTGGTATCCGGAAAACAGATGTGCTCTCATTATCGTACTTATCCCCGTTAAGTCAATACTACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                     <  1:260580/150‑1 (MQ=255)
       tATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATc                                                                                                                                                <  2:205857/150‑1 (MQ=255)
        aTCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCg                                                                                                                                               >  1:1691772/1‑150 (MQ=255)
         tCGCTGGTCTTCTTATCGTTAATTTCTTGTTTCGGGAGCAGATTTCCATATGTGCTGGTTTCCGGTAATCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGAAAGCACGAAAAAGGGAGCGATTATTTATCGc                                                                                                                                              >  1:304773/1‑150 (MQ=255)
         tCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGc                                                                                                                                              <  1:1371739/150‑1 (MQ=255)
            cTGGGCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTcc                                                                                                                                           <  1:722355/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)
                     tAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTta                                                                                                                                  >  1:1123922/1‑150 (MQ=255)
                     tAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTta                                                                                                                                  >  1:488595/1‑150 (MQ=255)
                     tAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTta                                                                                                                                  <  1:179920/150‑1 (MQ=255)
                      aaTCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTtat                                                                                                                                 >  2:1120885/1‑150 (MQ=255)
                       aTCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTtata                                                                                                                                <  2:1746495/150‑1 (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:973616/1‑150 (MQ=255)
                                ttCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAg                                                                                                                       >  1:1115832/1‑150 (MQ=255)
                                        tCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGtt                                                                                                               <  2:961768/150‑1 (MQ=255)
                                         cGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTa                                                                                                              >  2:380836/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)
                                                   tatCCATATGTGCTGGTATCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTg                                                                                                    <  2:1461137/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)
                                                          atGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattca                                                                                             <  2:1659833/150‑1 (MQ=255)
                                                             tgCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatca                                                                                          >  1:505747/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)
                                                                      ccGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgt                                                                                 >  2:1318336/1‑150 (MQ=255)
                                                                       cGGTAAACAGATGTGCGCTCGTTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtg                                                                                >  2:1294831/1‑150 (MQ=255)
                                                                            aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggtg                                                                           <  1:1365239/150‑1 (MQ=255)
                                                                            aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggtg                                                                           <  2:339381/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)
                                                                                        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)
                                                                                                      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)
                                                                                                         ccccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGa                                              >  1:845849/1‑150 (MQ=255)
                                                                                                            cGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCt                                           <  2:505747/150‑1 (MQ=255)
                                                                                                               tAAGTCAATACGACAACAAGCACGAAAAAGGGAGCGATGATTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCtttt                                        <  1:413742/150‑1 (MQ=255)
                                                                                                                    cAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGc                                   >  2:269097/1‑150 (MQ=255)
                                                                                                                     aaTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCg                                  >  1:1208070/1‑150 (MQ=255)
                                                                                                                     aaTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCg                                  <  2:663142/150‑1 (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)
                                                                                                                            cAGCAAGCACGAAAAAGGGATCGATTAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGAc                           <  1:258057/150‑1 (MQ=255)
                                                                                                                             agcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACt                          <  2:1104255/150‑1 (MQ=255)
                                                                                                                             agcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACt                          <  2:1418308/150‑1 (MQ=255)
                                                                                                                              gcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTg                         >  1:35501/1‑150 (MQ=255)
                                                                                                                                   caCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                    <  1:1253466/150‑1 (MQ=255)
                                                                                                                                   caCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                    >  2:773479/1‑150 (MQ=255)
                                                                                                                                    aCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAtt                   <  2:580315/150‑1 (MQ=255)
                                                                                                                                     cGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg                  <  2:485401/150‑1 (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:432597/2‑150 (MQ=255)
                                                                                                                                         caaGGGAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTc              >  1:968303/2‑150 (MQ=255)
                                                                                                                                            gggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag           <  1:1441377/150‑1 (MQ=255)
                                                                                                                                            gggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag           >  1:610740/1‑150 (MQ=255)
                                                                                                                                            gggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag           >  1:305914/1‑150 (MQ=255)
                                                                                                                                            gggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag           >  1:1030967/1‑150 (MQ=255)
                                                                                                                                             ggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga          >  1:771695/1‑150 (MQ=255)
                                                                                                                                             ggAGCGATAATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga          >  2:1337982/1‑150 (MQ=255)
                                                                                                                                                  gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGAc     >  1:1728516/1‑150 (MQ=255)
                                                                                                                                                     gAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc  >  1:1750693/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCC  >  NC_000913/1568766‑1569064

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: