Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,775 A→T 5.1% intergenic (+321/+16) 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,7750AT 5.1% 238.9 / 2.8 93intergenic (+321/+16)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (50/36);  new base T (3/2);  total (55/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

ACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC  >  NC_000913/357632‑357922
                                                                                                                                               |                                                                                                                                                   
aCGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTAGTAGGCCTGCTAAGCGTAGCGCAACGGGCAATTCCAGCCGCAGACCt                                                                                                                                                >  1:1139400/1‑149 (MQ=255)
  gCCATGTGCAAAAATGGCCTCGCGGCGGCTGTCTGCCAGCTGCAGTTTATGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtgt                                                                                                                                              <  2:415357/149‑1 (MQ=18)
   ccATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTCCGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCtgtg                                                                                                                                             <  2:1046600/149‑1 (MQ=255)
          cAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGGTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTTTcggcgg                                                                                                                                      >  1:389402/1‑149 (MQ=255)
          cAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCgg                                                                                                                                      >  1:476645/1‑149 (MQ=255)
           aaaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCCTAGGCGGGATAAGCGTAGCGCAGCAGGCAATTCCAGCCGCAGCCCGGTGTCAGGGGc                                                                                                                                     >  1:1072868/1‑149 (MQ=18)
                  ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTg                                                                                                                              >  1:1144511/1‑149 (MQ=255)
                   cTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGa                                                                                                                             <  1:507334/149‑1 (MQ=255)
                    ttCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTAGCGTGAt                                                                                                                            >  1:462861/1‑149 (MQ=255)
                    ttCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAt                                                                                                                            >  1:1168744/1‑149 (MQ=255)
                    ttCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAt                                                                                                                            <  2:112062/149‑1 (MQ=255)
                    ttCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTAAGCGGCTACCGGGAt                                                                                                                            >  2:579197/1‑149 (MQ=255)
                        cAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGAATGCGGCGTGAACGCCTTATCCGACCTACGGAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGGCCTGTGTCAGCGGCTACCGGGAATCAt                                                                                                                        >  2:351361/1‑149 (MQ=17)
                          gCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAttt                                                                                                                      >  1:151697/1‑149 (MQ=255)
                          gCGGCTGTCTGCCAGCTGTAGTTTATGCCGGAGGCGGCGCGAACGGCTTCTCGGGCCTACGTAGGGCCTTGACCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAttt                                                                                                                      <  1:1094446/149‑1 (MQ=17)
                           cGGCTGTCTGCCCGCTGTAGTTTATGCCGGATGCCGCGTGAACGCCTTATCCGGCCTCCGTAGCGCACTGACCTCGTAGGCCTGCTGAGCGTCGCGCATCCGACAATTCCAGCCGCCCCCGTCTGTCAGCGGCCAGCGGGAGTCAtgtg                                                                                                                     >  2:684106/1‑146 (MQ=9)
                           cGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCGGACCTGCGTCAGGGGCTACCGTGAGTCATATc                                                                                                                     >  2:703746/1‑149 (MQ=14)
                               tGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGcc                                                                                                                 >  2:483763/1‑149 (MQ=255)
                               cgTCTGCCAGCTGTAGTTTATGCCGGATGCGCCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCCATTCCAGCCGCAGACCGGTGTCAGCGGCTACCGTGATCCATTTCCGcc                                                                                                                 >  2:948698/2‑149 (MQ=255)
                                 tCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGTCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCcaa                                                                                                               >  1:461710/1‑149 (MQ=255)
                                 tCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCcaa                                                                                                               <  1:354507/149‑1 (MQ=255)
                                 tCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTAGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGAGTCATTTCCGCcaa                                                                                                               >  1:716236/1‑149 (MQ=255)
                                   tGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCcaaca                                                                                                             >  1:209494/1‑149 (MQ=255)
                                   tGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCcaaca                                                                                                             >  1:1204362/1‑149 (MQ=255)
                                   tGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGCCCTACGTAGCGCCCTCCACTCGTACGCCTGATACGCGTAGCGCATCAGGCAATTCAAGCCGTAGACCTGTGTCAGCGGCTACCGGGATTCAGTTCCGccacca                                                                                                             >  2:758727/1‑149 (MQ=14)
                                          tGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCgcgc                                                                                                      >  2:553405/1‑149 (MQ=40)
                                            tAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCGACAACCgcgcct                                                                                                    >  2:160468/1‑147 (MQ=255)
                                            tAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                    >  2:160037/1‑149 (MQ=38)
                                              gTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAttt                                                                                                  <  2:226496/149‑1 (MQ=38)
                                               tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAttca                                                                                                 >  1:1171927/1‑147 (MQ=14)
                                               tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTa                                                                                                 >  1:678376/1‑149 (MQ=39)
                                               tgtATGCCGGAGGCTGCGTGAATGCCTTATCCGGCCTACGTAGAGCACTGAACTAGTCGGCCAGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCTAACAACCGCGCATTTa                                                                                                 <  1:820382/147‑1 (MQ=11)
                                                    gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCt                                                                                                                                                >  2:264406/1‑97 (MQ=37)
                                                    gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCt                                                                                                                                                <  1:264406/97‑1 (MQ=37)
                                                    gCCGGAGGCGGCCTGAACGCCTTAGCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCaa                                                                                            <  1:1136113/149‑1 (MQ=255)
                                                       gtttGCGGTGTGGACGCCTCATCCGGCCAACGTTACGCTCTGACCACCTCGGCCTGCTAAGCGTAGCCCATCACGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGc                                                                                         <  2:191167/146‑1 (MQ=14)
                                                       gggcgcggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCGACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGc                                                                                         <  1:630011/145‑1 (MQ=255)
                                                        gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGcc                                                                                        <  1:147617/149‑1 (MQ=255)
                                                        gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGcc                                                                                        <  2:1171927/149‑1 (MQ=255)
                                                            cggcgTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAtgt                                                                                    >  1:982350/1‑149 (MQ=255)
                                                                  gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCCACGCCATGTGCaaaa                                                                              >  2:638234/1‑149 (MQ=255)
                                                                  gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaa                                                                              >  2:637886/1‑149 (MQ=255)
                                                                   aaCGGCTTACCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCACCAACAACCGCGTATTTATCCAACGCCATGTGCaaaaa                                                                             <  1:670268/149‑1 (MQ=255)
                                                                   aaCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaaa                                                                             <  2:101943/149‑1 (MQ=255)
                                                                   aaCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaaa                                                                             <  1:335076/149‑1 (MQ=255)
                                                                     cGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATg                                                                           >  1:144399/1‑149 (MQ=255)
                                                                        cTTATCCGGCCTACGTAGAGCACTTAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCt                                                                        >  1:967599/1‑149 (MQ=255)
                                                                            tCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTcgc                                                                    >  1:836360/1‑149 (MQ=255)
                                                                            tCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTcgc                                                                    >  1:645301/1‑149 (MQ=255)
                                                                             ccGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTcgcg                                                                   >  2:1194761/1‑149 (MQ=255)
                                                                              cGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCgcgg                                                                  >  1:620758/1‑149 (MQ=255)
                                                                                  cTACGTAGAGCACTGAACTCGTAGGCCTGCTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCGTGTGCAAAAATGCCTTCGcggcgga                                                              >  1:172204/1‑148 (MQ=255)
                                                                                     cGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGt                                                           <  2:678376/149‑1 (MQ=255)
                                                                                          agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGcc                                                      <  1:418442/149‑1 (MQ=255)
                                                                                             aCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGc                                                   >  1:76510/1‑149 (MQ=255)
                                                                                               tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTCATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                 >  1:1061519/1‑149 (MQ=25)
                                                                                                 aaCTCGTAGGCCTGCTAAGCGTAGCGCATCAGGCAATTCCAGCAGCAGACCTGGGGCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCCTGTGCAAAAAGCCCTTCGCGGCGGCTGTCGGACAACTagt                                               >  2:505158/1‑147 (MQ=14)
                                                                                                 aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCGGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtt                                               >  2:313619/1‑149 (MQ=255)
                                                                                                 aaCTCGTAGGCCTGATAAGCGTAACGCAGATGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtt                                               >  2:924291/1‑149 (MQ=255)
                                                                                                   cTCGTAGGCCTGATAAGCGTAGCGCATCAGGAAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtttt                                             <  1:1165171/149‑1 (MQ=255)
                                                                                                    tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTGttttt                                            >  1:1493/1‑149 (MQ=255)
                                                                                                      gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                          >  2:730822/1‑149 (MQ=255)
                                                                                                      gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                          <  2:1144511/149‑1 (MQ=255)
                                                                                                       caGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                         <  2:826945/148‑1 (MQ=255)
                                                                                                          gcttGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCTTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                      <  2:811389/146‑1 (MQ=255)
                                                                                                          gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAttt                                                                                                                      <  2:201981/69‑1 (MQ=255)
                                                                                                          gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCAttt                                                                                                                      <  1:1247535/69‑1 (MQ=255)
                                                                                                             tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                   >  1:77649/1‑149 (MQ=255)
                                                                                                              gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                  >  2:144296/1‑149 (MQ=255)
                                                                                                              gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCCcaac                                  >  1:478561/1‑149 (MQ=255)
                                                                                                              gATAAGCGTAACGCATCAGGCAATTCCAGCCACAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                  >  2:763008/1‑149 (MQ=255)
                                                                                                                tacgCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                <  1:517034/146‑1 (MQ=255)
                                                                                                                   gCGGAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                             <  2:447549/149‑1 (MQ=255)
                                                                                                                     gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTGCCGGCGCAAAAAAAccg                           >  1:296388/1‑149 (MQ=255)
                                                                                                                            aTCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                     <  2:399273/132‑1 (MQ=255)
                                                                                                                            aTCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                     >  1:399273/1‑132 (MQ=255)
                                                                                                                              cAGGCAATTCCAGCTGCAGACCTGTGTCAGCGGCTATCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCATTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAACACCGCCGTtctct                  <  1:915613/149‑1 (MQ=255)
                                                                                                                              cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                  <  1:1297396/149‑1 (MQ=255)
                                                                                                                               aGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                 <  1:295185/149‑1 (MQ=255)
                                                                                                                                ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCTTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                >  2:853591/1‑149 (MQ=255)
                                                                                                                                 ccAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa               <  2:461710/148‑1 (MQ=255)
                                                                                                                                    tttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCTACACCCGCGCATTTATCCAACGCCATGTGCAAAAAGGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagta            <  2:1061519/148‑1 (MQ=255)
                                                                                                                                    aTTCCAGGCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagta            <  2:656856/149‑1 (MQ=255)
                                                                                                                                    aTTCCAGCCGCAGACCTGTGTCAGCAGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCTCCGCAACAAAACCGCCGTTCTCTCCagta            <  1:790566/149‑1 (MQ=255)
                                                                                                                                           ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCCGTCTGCCAGCTATTTTTCCGCCGCAACAAACCGCCGTTCTCTCCAGTAGTGGCggg    >  2:1125985/1‑149 (MQ=255)
                                                                                                                                            cGCAGACCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg    <  1:351361/149‑1 (MQ=255)
                                                                                                                                             gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg   >  1:1116629/1‑149 (MQ=255)
                                                                                                                                             gCAGACCTGGGTCAGCGGCTACCGTGATTCCTTTCCGCCAACATCCGCGCATTCATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGGCTGCCAGCTATTTTTCCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg   <  2:780417/149‑1 (MQ=21)
                                                                                                                                              cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTCTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGGAGTGGCGGGGc  >  1:533442/1‑149 (MQ=255)
                                                                                                                                              cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCCACAAAACCGCCGTTCACTCCGGTTGTGGCGGGGc  >  2:684060/1‑149 (MQ=255)
                                                                                                                                              cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc  <  1:1069749/149‑1 (MQ=255)
                                                                                                                                              cAGACCTGTGTCAGCGGCTACCGTGATTCACTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACACAACCGCCGTTCTCTCCAGTAGTGGCGGGGc  >  1:673005/1‑149 (MQ=255)
                                                                                                                                              cAGACCTGTGTCAGCGGCTACCGTGATTCAATTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGACGCTCTCTCCAGTAGTGGCGGGGc  >  1:672944/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
ACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC  >  NC_000913/357632‑357922

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