Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,206,629 C→G 6.1% pseudogene (499/792 nt) jayE → e14 prophage; putative protein JayE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,206,6290CG 6.1% 321.4 / 7.4 95pseudogene (499/792 nt)jayEe14 prophage; putative protein JayE
Reads supporting (aligned to +/- strand):  ref base C (48/41);  new base G (3/3);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTAAAGAAGTGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTGATGGTTATCC  >  NC_000913/1206480‑1206777
                                                                                                                                                     |                                                                                                                                                    
cTAAAGAAGTGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCAc                                                                                                                                                      >  2:618465/1‑150 (MQ=255)
cTAAAGAAGTGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCAc                                                                                                                                                      >  2:799445/1‑150 (MQ=255)
 tAAAGAAGTGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACt                                                                                                                                                     <  1:162893/150‑1 (MQ=255)
  aaaGAAGTGCCCGGCATTACCCGAGCATGGACATACCGTAACTTGATCGGAACGGGAACTCTCGGTGTGATGATTCCCAGCAGTGACCTGATTAATACCATTCCGGAAGAATCAACGGAAACGGCTGCAAGACAACATATCGGACCACtt                                                                                                                                                    >  2:299736/1‑149 (MQ=255)
     gaagTGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGTGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAATGGCGGCAAGACAACATATCGGGCCACTGGcc                                                                                                                                                 <  2:968501/150‑1 (MQ=255)
     gaagTGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGcc                                                                                                                                                 >  1:779721/1‑150 (MQ=255)
       agTGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGcccc                                                                                                                                               <  1:555146/150‑1 (MQ=255)
        gTGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCg                                                                                                                                              >  1:1015070/1‑150 (MQ=255)
         tGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCgg                                                                                                                                             >  1:435488/1‑150 (MQ=255)
         tGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCgg                                                                                                                                             >  1:513574/1‑150 (MQ=255)
          gCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGt                                                                                                                                            >  1:344245/1‑150 (MQ=255)
          gCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACCggggccagt                                                                                                                                            >  1:475540/1‑143 (MQ=255)
           cccGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCAGTGGcccgatat                                                                                                                                           <  1:284749/150‑6 (MQ=255)
           cccGGCATTACCCGCGCATGCACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGTCAAGACAACATATCGGGCCAGTGGcccgatat                                                                                                                                           <  1:13710/150‑6 (MQ=255)
            ccGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCAGTGGcccgatatg                                                                                                                                          <  1:709408/150‑7 (MQ=255)
            ccGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTgg                                                                                                                                          >  1:275195/1‑150 (MQ=255)
              ggCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACCGGGGCCAGTGGcc                                                                                                                                        <  2:475540/150‑2 (MQ=255)
               gCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAg                                                                                                                                       <  1:967819/150‑1 (MQ=255)
                 aTTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGc                                                                                                                                     >  2:22368/1‑150 (MQ=255)
                 aTTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGc                                                                                                                                     <  2:275195/150‑1 (MQ=255)
                 aTTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGc                                                                                                                                     <  1:926250/150‑1 (MQ=255)
                    aCCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGctct                                                                                                                                  >  1:364337/1‑150 (MQ=255)
                     cccGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGAGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTg                                                                                                                                 >  2:1183089/1‑150 (MQ=255)
                       cgcgCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGAt                                                                                                                               >  2:33057/1‑150 (MQ=255)
                         cgcATGGACATACCGTCAATTGATGGGAACGGCAACTGTCGGTGTGATGATAGCAAGCAGTGACCTGATTAATCCCATTCAGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGAttt                                                                                                                             >  2:1472333/1‑150 (MQ=255)
                            aTGGACATACCGTCACTTGATGGGATCGTGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGAAAGACAACTTATCGGTCCACTGGCCCCGGTGGCAGGCTCTGATTTGTa                                                                                                                          <  2:696895/150‑1 (MQ=255)
                                aCATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTAtgtg                                                                                                                      >  1:977218/1‑150 (MQ=255)
                                 cATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATgtgt                                                                                                                     >  2:388451/1‑150 (MQ=255)
                                  aTACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGTCGGCAAGAGAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGtt                                                                                                                    <  1:463042/150‑1 (MQ=255)
                                   tACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTc                                                                                                                   >  2:346295/1‑150 (MQ=255)
                                    aCCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCa                                                                                                                  <  2:1277635/150‑1 (MQ=255)
                                    aCCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCa                                                                                                                  <  2:364337/150‑1 (MQ=255)
                                     ccGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAg                                                                                                                 <  1:909311/150‑1 (MQ=255)
                                     ccGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAg                                                                                                                 <  1:152200/150‑1 (MQ=255)
                                      catCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAgg                                                                                                                >  2:874111/3‑150 (MQ=255)
                                       gTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGc                                                                                                               <  1:33057/150‑1 (MQ=255)
                                       gTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGc                                                                                                               >  2:451055/1‑150 (MQ=255)
                                       gTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGACAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGc                                                                                                               <  2:325825/150‑1 (MQ=255)
                                         cACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCg                                                                                                             >  2:722903/1‑150 (MQ=255)
                                           cTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGt                                                                                                           >  1:395649/1‑150 (MQ=255)
                                           cTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGt                                                                                                           >  1:639416/1‑150 (MQ=255)
                                            ttGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTg                                                                                                          <  1:799445/150‑1 (MQ=255)
                                             tGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCCCTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTgg                                                                                                         >  2:934660/1‑150 (MQ=255)
                                              gATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGc                                                                                                        >  2:1293358/1‑150 (MQ=255)
                                                     acgggaacTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACg                                                                                                 >  1:1194249/1‑150 (MQ=255)
                                                     acgggaacTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAAGCCGGTGGCACATACg                                                                                                 >  2:295487/1‑150 (MQ=255)
                                                      cgggaacTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACgg                                                                                                >  1:188732/1‑150 (MQ=255)
                                                        ggaacTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTg                                                                                              <  1:832459/150‑1 (MQ=255)
                                                            cTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGAtt                                                                                          >  2:617434/1‑150 (MQ=255)
                                                                cGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCa                                                                                      <  1:874111/150‑1 (MQ=255)
                                                                cGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCa                                                                                      <  1:598020/150‑1 (MQ=255)
                                                                cGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGAGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCa                                                                                      <  2:1207864/150‑1 (MQ=255)
                                                                   tgtgATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCatat                                                                                   <  2:1181725/150‑1 (MQ=255)
                                                                     tgatgaTTGCTAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATcc                                                                                 >  1:780737/1‑150 (MQ=255)
                                                                     tgatgaTTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATcc                                                                                 <  1:1084232/150‑1 (MQ=255)
                                                                             gCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACg                                                                         >  1:1556575/1‑150 (MQ=255)
                                                                               cagcagTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGTTTTTGATATCCGCGTGACGcc                                                                       <  2:397848/150‑1 (MQ=255)
                                                                               cagcagTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGcc                                                                       >  1:1443496/1‑150 (MQ=255)
                                                                                 gcagTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCgg                                                                     >  1:339600/1‑150 (MQ=255)
                                                                                    gTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGaca                                                                  >  2:1151103/1‑150 (MQ=255)
                                                                                           gATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGa                                                           <  2:669516/150‑1 (MQ=255)
                                                                                                 tCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTACATATCCGCGTGACGCCGGACACACCAGAAATACg                                                     >  1:1604022/1‑150 (MQ=255)
                                                                                                   ccATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACggg                                                   >  2:505273/1‑150 (MQ=255)
                                                                                                     aTTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCt                                                 <  2:1197211/150‑1 (MQ=255)
                                                                                                      ttCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCa                                                                                      <  2:1305665/112‑1 (MQ=255)
                                                                                                      ttCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCa                                                                                      >  1:1305665/1‑112 (MQ=255)
                                                                                                       tCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGc                                               <  2:513574/150‑1 (MQ=255)
                                                                                                        ccGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGcc                                              <  1:451055/150‑1 (MQ=255)
                                                                                                        ccGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGcc                                              >  2:1052980/1‑150 (MQ=255)
                                                                                                         cGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCa                                             >  2:389177/1‑150 (MQ=255)
                                                                                                          ggAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCAt                                            <  1:299736/150‑1 (MQ=255)
                                                                                                            aagaaTCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGGTCTGATTTGTTTGTGTTCAGGCCGGTGTCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTa                                          <  2:435488/150‑1 (MQ=255)
                                                                                                                aTCAACGGAAACGGCGGCAAGACAACATATCGGGACACTGGCCCCGGTGGCAGGCTCTGATTTGTATCTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACcgc                                      >  1:561848/1‑150 (MQ=255)
                                                                                                                  cAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCgg                                    <  1:22368/150‑1 (MQ=255)
                                                                                                                  cAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCgg                                    >  2:955036/1‑150 (MQ=255)
                                                                                                                      ggAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGtt                                >  1:1079644/1‑150 (MQ=255)
                                                                                                                       gAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTg                               >  1:156289/1‑150 (MQ=255)
                                                                                                                       gAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTg                               <  1:1183089/150‑1 (MQ=255)
                                                                                                                       gAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTg                               >  1:595363/1‑150 (MQ=255)
                                                                                                                            ggcggcAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGcgttc                          >  2:1323098/1‑150 (MQ=255)
                                                                                                                            ggcggcAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGcgttc                          >  2:1317557/1‑150 (MQ=255)
                                                                                                                                 cAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTAGCGTTCGTTcc                     >  1:1041913/1‑150 (MQ=255)
                                                                                                                                  aaGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCct                    <  2:779721/150‑1 (MQ=255)
                                                                                                                                   aGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCctg                   <  2:156289/150‑1 (MQ=255)
                                                                                                                                        acaTATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGTCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCg              <  1:722903/150‑1 (MQ=255)
                                                                                                                                        acaTATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCg              <  2:639416/150‑1 (MQ=255)
                                                                                                                                         caTATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGTCACTTTCGGTTTATTTTCATATGGGCTTGACGCCGTTCTCACCAGTAATACTGTCTGTCATTACAGCTGTGTTTAGTTGTTTACTGCTGCGt             <  2:331620/150‑1 (MQ=255)
                                                                                                                                         caTATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACTCCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGt             <  1:1052980/150‑1 (MQ=255)
                                                                                                                                          caccggggCCAGTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTAACGCGGAGTTGCGTTCGTTCCTGCTGCGTg            >  2:795636/6‑150 (MQ=255)
                                                                                                                                          atatCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTg            >  2:1015674/1‑150 (MQ=255)
                                                                                                                                           tatCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTGa           >  2:1564851/1‑150 (MQ=255)
                                                                                                                                           accggggCCAGTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTGa           >  2:402277/5‑150 (MQ=255)
                                                                                                                                           accggggCCAGTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTGa           >  2:396349/5‑150 (MQ=255)
                                                                                                                                            atCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTGAt          <  2:780737/150‑1 (MQ=255)
                                                                                                                                             tCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTGATg         <  2:60/150‑1 (MQ=255)
                                                                                                                                             tCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTGATg         <  1:346295/150‑1 (MQ=255)
                                                                                                                                                    aCTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTGATGGTTATcc  >  1:73482/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTAAAGAAGTGCCCGGCATTACCCGCGCATGGACATACCGTCACTTGATGGGAACGGGAACTGTCGGTGTGATGATTGCCAGCAGTGACCTGATTAATCCCATTCCGGAAGAATCAACGGAAACGGCGGCAAGACAACATATCGGGCCACTGGCCCCGGTGGCAGGCTCTGATTTGTATGTGTTCAGGCCGGTGGCACATACGGTGGATTTTCATATCCGCGTGACGCCGGACACACCAGAAATACGGGCTGCCATTACCGCGGAGTTGCGTTCGTTCCTGCTGCGTGATGGTTATCC  >  NC_000913/1206480‑1206777

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