Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,778 C→T 7.1% intergenic (+324/+13) codA → / ← cynR cytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,7780CT 7.1% 244.5 / 2.6 90intergenic (+324/+13)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (54/28);  new base T (3/3);  total (58/32)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG  >  NC_000913/357635‑357925
                                                                                                                                               |                                                                                                                                                   
taaTGTTCAATACTGCCTTCGGAGCGGGTTTCTTGCTTGTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTgtgt                                                                                                                                               <  2:1156331/148‑1 (MQ=14)
    ttgCAAAATTGCCTGCGCATCTGCTTGCTTCCATCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGc                                                                                                                                           <  2:503688/149‑1 (MQ=14)
     tgCAAAAATGCCTTCGCAGCGGCGGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCg                                                                                                                                          <  2:1464048/150‑1 (MQ=255)
      ccAAAAATCCCTTCGCAGCGTCTTTCTCCCTGCTTTTTTTTATGCCTGATGCGGCGTTAACGCCTTATCCGGCCTACGTACAGCACTTAACTCGTAGGCCTGATACGCTTAGCGAATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCgg                                                                                                                                         <  2:757708/149‑1 (MQ=14)
       cAAAAATGCCTTCGCAGCGGCTGTCTGGCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGc                                                                                                                                        <  1:1147488/150‑1 (MQ=255)
       cAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGc                                                                                                                                        <  2:1895974/150‑1 (MQ=255)
       cAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTATATTTAATGCCGCAAAATGCCTGAAAGCCTTAACCGGCCTAAGTAGAGAACGGATATAGTAGGCATGATAAGCGTAGCGCAAGATGCAATTAAAAACTAAGACGTGTGTCAGCggg                                                                                                                                        >  2:1344039/1‑149 (MQ=255)
       cAAAAATGCCTTCGCAGCGGCGGTCTTCCATCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGc                                                                                                                                        <  2:1876051/150‑1 (MQ=255)
         aaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACTCCTTATCCGGCATACGTAGATCACTGAAATCGTAGGCCTCATACGCGTAGAGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTa                                                                                                                                      >  2:524987/1‑150 (MQ=18)
               ccTTCGCAGCGGCTGTCTGCCAGCTGTAGATTATGCAGGATGCGGCGTGAACGCCTTATACGGCCTACGTAGAGCACTGAACTCGTAGTCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGa                                                                                                                                >  2:854013/1‑150 (MQ=255)
                cTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAt                                                                                                                               >  2:1623209/1‑150 (MQ=255)
                cTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGAGTAGCGCATCAGGCAATTCCCGCCGCCGCCCTGTGTCAGCGGCTACCGTGAt                                                                                                                               >  1:1460133/1‑150 (MQ=255)
                 ttCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTAGCGTGAtt                                                                                                                              >  1:669560/1‑150 (MQ=255)
                 ttCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACATGTGTCAGCGGCTACCGTGAtt                                                                                                                              >  2:1521250/1‑150 (MQ=255)
                 ttCGCAGCGGCTGTCTGCCAGCTGCAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAtt                                                                                                                              <  1:1664980/150‑1 (MQ=255)
                   cGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCa                                                                                                                            <  2:827806/150‑1 (MQ=255)
                   cGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCa                                                                                                                            >  2:943184/1‑150 (MQ=255)
                         ggCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAGTCATTTCCg                                                                                                                      >  2:1700916/1‑150 (MQ=255)
                            tGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCca                                                                                                                   >  1:1128668/1‑150 (MQ=255)
                                      cTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAGTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCCGTTCCGCCACCAACCgcgc                                                                                                         >  1:1321723/1‑150 (MQ=255)
                                       cttAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCa                                                                                                        <  2:1882473/148‑1 (MQ=25)
                                        gTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTCCCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                       >  1:828491/1‑150 (MQ=38)
                                        gTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                       >  2:1931561/1‑150 (MQ=38)
                                          aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGCAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAttt                                                                                                     >  2:407729/1‑150 (MQ=38)
                                              tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATcc                                                                                                 >  1:1265345/1‑150 (MQ=40)
                                                tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTCGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCaa                                                                                               >  1:687239/1‑150 (MQ=255)
                                                tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACACCCGCGCATTTATCCAAc                                                                                              >  2:1517263/1‑150 (MQ=255)
                                                  ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGCCCTGTGTCAGCGGCGACCGTGATTCATTTCCGCCAACAACCGCGCAGTTATCCAACg                                                                                             >  1:325164/1‑150 (MQ=255)
                                                    ggataaggcgTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGcc                                                                                           >  1:672495/7‑150 (MQ=14)
                                                         cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAtgtg                                                                                      <  2:1258111/150‑1 (MQ=37)
                                                          ggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGc                                                                                     >  2:1522992/1‑150 (MQ=255)
                                                           gcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCCTGTGCa                                                                                    >  1:1796411/1‑150 (MQ=255)
                                                            cgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaa                                                                                   >  2:1966640/1‑150 (MQ=255)
                                                            cgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaa                                                                                   <  2:532872/150‑1 (MQ=255)
                                                                aaCGCCTTATCCGGACTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAt                                                                               >  2:1867038/1‑150 (MQ=255)
                                                                 aCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCCACGCCATGTGCAAAAATg                                                                              >  1:1910569/1‑150 (MQ=255)
                                                                  tgCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGc                                                                             >  1:1448984/2‑150 (MQ=255)
                                                                  cGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAAACGCGCATTTATCCAACGCCATGTGCAAAAATGc                                                                             >  1:1542846/1‑150 (MQ=255)
                                                                  cGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGATACCGTGATTCAGTTCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGc                                                                             >  2:47816/1‑150 (MQ=255)
                                                                    ccTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCt                                                                           >  2:467930/1‑150 (MQ=255)
                                                                       tATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTcg                                                                        >  2:1148666/1‑150 (MQ=255)
                                                                       tATCCGGCCTACGTAGAGAACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTcg                                                                        >  1:1994404/1‑150 (MQ=255)
                                                                          ctggCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGAcc                                                                                                                                                    <  1:1017004/69‑1 (MQ=38)
                                                                          ctggCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGAcc                                                                                                                                                    >  2:1017004/3‑71 (MQ=38)
                                                                          ccGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCgcgg                                                                     >  2:1239392/1‑150 (MQ=255)
                                                                             gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAtt                                                                                                                              >  1:1079141/1‑90 (MQ=38)
                                                                             gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAtt                                                                                                                              <  2:1079141/90‑1 (MQ=38)
                                                                              ccTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGc                                                                                                                                        >  2:1347502/1‑79 (MQ=39)
                                                                              ccTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGc                                                                                                                                        <  1:1347502/79‑1 (MQ=39)
                                                                              ccTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCggcggc                                                                 <  2:1128668/150‑1 (MQ=255)
                                                                               cTACGTAGAGCTATGATCTCGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCTACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCt                                                                <  1:1027930/150‑1 (MQ=14)
                                                                                tACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTg                                                               >  1:1398310/1‑150 (MQ=255)
                                                                                  cGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTc                                                             <  1:807166/150‑1 (MQ=255)
                                                                                    tAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTGATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTg                                                           >  1:1877447/1‑150 (MQ=12)
                                                                                    tAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTTTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTg                                                           <  1:854013/150‑1 (MQ=255)
                                                                                    tAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGgctggctg                                                           <  1:1615353/150‑1 (MQ=255)
                                                                                        gCACTGAACTCGTAGGCCTGATAAGCGTCGCGCATCAGGCAACTCCAGCCGCAGAGCTGTGGCTGCGTCTACCGTGATGCATTTCCGCCAACAACCGCGCATTCATCCCTCGCCATGTCCAAAAAGGCCTTCGCGGCGGCTGGCTGCCAg                                                       <  1:2001063/150‑1 (MQ=11)
                                                                                        gCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCACCACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAg                                                       >  2:1608209/1‑150 (MQ=255)
                                                                                           cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                    >  2:893122/1‑150 (MQ=255)
                                                                                           cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                    <  1:1522992/150‑1 (MQ=255)
                                                                                             gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGcc                                                                                                                    >  1:166359/1‑84 (MQ=18)
                                                                                             gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGcc                                                                                                                    <  2:166359/84‑1 (MQ=255)
                                                                                                cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                               >  2:1895825/1‑150 (MQ=255)
                                                                                                 tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                              <  2:672495/150‑1 (MQ=255)
                                                                                                     aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                          >  2:410793/1‑150 (MQ=255)
                                                                                                      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                         >  1:583089/1‑150 (MQ=255)
                                                                                                      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                         >  1:486239/1‑150 (MQ=255)
                                                                                                      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATACATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                         >  2:1168623/1‑150 (MQ=255)
                                                                                                        ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTGCCGACGca                                       >  2:1363955/1‑150 (MQ=255)
                                                                                                          tGATAAGCGTCGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACCCCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                     >  2:1801708/1‑150 (MQ=255)
                                                                                                           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCAGCCAGCTATTTTTCCGCCGcaaca                                    >  1:458187/1‑150 (MQ=255)
                                                                                                             tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCACCCACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATGTTTCCGCCGCAACaaa                                  >  1:1993988/1‑150 (MQ=255)
                                                                                                                 cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCCTATCCGCCAACTACCGCGCATTTATCCAACTCCTAGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                              <  1:1344039/150‑1 (MQ=18)
                                                                                                                  gTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCa                                                        <  2:88606/123‑1 (MQ=25)
                                                                                                                  gTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCa                                                        >  1:88606/1‑123 (MQ=12)
                                                                                                                  gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGcc                                                                                                                    <  2:123234/63‑1 (MQ=255)
                                                                                                                  gTAGCGCATCAGGCAATTCCAGCCGCAGAACTGTGTCAGCGGCTACCGTGATTCATTTCCGcc                                                                                                                    >  1:123234/1‑63 (MQ=255)
                                                                                                                    tgcgaATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTAGCGTGATACATTTCCGCGAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                           <  1:1921041/149‑1 (MQ=255)
                                                                                                                             ggCAATTCCAGCCGCAGACCTGTGTCAGTGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa                  >  2:1480332/1‑150 (MQ=255)
                                                                                                                              gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag                 <  2:1265345/150‑1 (MQ=255)
                                                                                                                                  ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATATGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCCATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt             >  1:906388/1‑150 (MQ=255)
                                                                                                                                    ttagcagcAGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACATCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg           <  1:1846021/148‑1 (MQ=255)
                                                                                                                                    ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTACTCTCCAGTAGTgg           >  2:1001854/1‑150 (MQ=255)
                                                                                                                                    ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGAGTCATTTCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCCGCTATTTTTCCGCCGCAACCAAACCGCCGTTCTCGCCAGTAGTgg           >  1:962093/1‑150 (MQ=21)
                                                                                                                                     cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc          >  1:811882/1‑150 (MQ=255)
                                                                                                                                     cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc          >  1:463181/1‑150 (MQ=255)
                                                                                                                                         cGCTGATCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg      <  1:1650335/150‑1 (MQ=255)
                                                                                                                                         cGCAGCCCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACAGCCGCGCATGTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCGGTCTGCGGGCTATTTTTCCGCCGCCCCACACCCGCCGTTCTCTCGCGTCGTGGCgggg      >  1:103689/1‑150 (MQ=18)
                                                                                                                                         cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAACACCGCCGTTCTCTCCAGTAGTGGCgggg      >  1:1481275/1‑150 (MQ=255)
                                                                                                                                         cGCAGACCTGTGGCAGCGGCTACCGTGATTCAATTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCTCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg      <  1:1342090/150‑1 (MQ=255)
                                                                                                                                           cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGAATTTATCCAACGCCATGTGCAAAAATGCATTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGGGGTCTCCCGAAGTGGCGGGGCa    >  2:379392/1‑150 (MQ=21)
                                                                                                                                             gACCTGCCTCACCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag  <  1:1715675/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG  >  NC_000913/357635‑357925

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: