Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,250,460 T→C 5.7% intergenic (+84/‑102) yqjK → / → yqjF conserved protein YqjK/DoxX family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,250,4600TC 5.7% 214.1 / 2.7 70intergenic (+84/‑102)yqjK/yqjFconserved protein YqjK/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base C (2/2);  total (36/34)
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.11e-01

TGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATG  >  NC_000913/3250313‑3250609
                                                                                                                                                   |                                                                                                                                                     
tGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTGTCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTc                                                                                                                                                     <  2:669024/150‑1 (MQ=255)
     tATGGAGCGCCTGGCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGTCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTcaccac                                                                                                                                                <  1:1517049/150‑1 (MQ=255)
        ggAGCGCCTGGCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGtt                                                                                                                                             >  1:105520/1‑150 (MQ=255)
         gAGCGCCTGGCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGttt                                                                                                                                            >  1:63250/1‑150 (MQ=255)
         gAGCGCCTGGCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGttt                                                                                                                                            <  1:1626845/150‑1 (MQ=255)
         gAGCGCCTGGCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGttt                                                                                                                                            >  2:1248216/1‑150 (MQ=255)
                tGGCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATtc                                                                                                                                     <  1:839054/150‑1 (MQ=255)
                 ggCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATtct                                                                                                                                    >  2:865514/1‑150 (MQ=255)
                  gCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATtctc                                                                                                                                   <  1:705262/150‑1 (MQ=255)
                   cGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATtctct                                                                                                                                  <  2:1163932/150‑1 (MQ=255)
                      cTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCACTCACCACGTTTATGATTCTCTCCa                                                                                                                               >  2:1095426/1‑150 (MQ=255)
                      cTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCACTCACCACGTTTATGATTCTCTCCa                                                                                                                               >  2:1101067/1‑150 (MQ=255)
                      cTGGTTAAAACGAACCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCa                                                                                                                               >  2:636377/1‑150 (MQ=255)
                       tGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCAt                                                                                                                              >  1:1078089/1‑150 (MQ=255)
                        ggTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATc                                                                                                                             >  2:1161842/1‑150 (MQ=255)
                                     cTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAacgacg                                                                                                                <  2:1595085/150‑1 (MQ=255)
                                     cTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAacgacg                                                                                                                <  2:140100/150‑1 (MQ=255)
                                     cTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAacgacg                                                                                                                >  1:1445184/1‑150 (MQ=255)
                                           cagcaACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATa                                                                                                          >  2:1397734/1‑150 (MQ=255)
                                            agcaACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGAAACGACGCTAATAc                                                                                                         >  2:1441907/1‑150 (MQ=255)
                                             gcaACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCACTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATAcc                                                                                                        >  2:681609/1‑150 (MQ=255)
                                               aaCAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACcgc                                                                                                      >  2:1296992/1‑150 (MQ=255)
                                                aCAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACcgcg                                                                                                     >  1:862650/1‑150 (MQ=255)
                                                  aGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCACTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGcc                                                                                                   <  1:1101067/150‑1 (MQ=255)
                                                  aGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCACTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGcc                                                                                                   <  1:1095426/150‑1 (MQ=255)
                                                     ttCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCAtt                                                                                                <  1:169263/150‑1 (MQ=255)
                                                      tCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTg                                                                                               <  2:933197/150‑1 (MQ=255)
                                                       tgcgGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGc                                                                                              >  2:251025/2‑150 (MQ=255)
                                                        gcgGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGca                                                                                             <  1:1083138/150‑1 (MQ=255)
                                                         cgGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGcac                                                                                            <  2:882988/150‑1 (MQ=255)
                                                         cgGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGcac                                                                                            <  1:869641/150‑1 (MQ=255)
                                                          ggTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGcaca                                                                                           <  1:485967/150‑1 (MQ=255)
                                                            ttAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACaaa                                                                                         >  1:420606/1‑150 (MQ=255)
                                                              aaCTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACaaaaa                                                                                       <  1:1428917/150‑1 (MQ=255)
                                                               aCTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACaaaaaa                                                                                      >  1:61922/1‑150 (MQ=255)
                                                               aCTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACaaaaaa                                                                                      >  1:1362461/1‑150 (MQ=255)
                                                                cTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACaaaaaaa                                                                                     >  1:381658/1‑150 (MQ=255)
                                                                     ctctGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATc                                                                                <  2:1319450/150‑1 (MQ=255)
                                                                         ggCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATcagca                                                                            >  1:1365281/1‑150 (MQ=255)
                                                                             ggAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAacgacg                                                                                                                <  1:1176658/110‑1 (MQ=255)
                                                                             ggAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAacgacg                                                                                                                >  2:1176658/1‑110 (MQ=255)
                                                                              gAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTg                                                                       <  1:1029897/150‑1 (MQ=255)
                                                                              gAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTg                                                                       <  1:681609/150‑1 (MQ=255)
                                                                                gCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAg                                                                     >  2:1347340/1‑150 (MQ=255)
                                                                                  cATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTg                                                                   >  2:244659/1‑150 (MQ=255)
                                                                                      ccGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGgcgc                                                               <  1:1476114/150‑1 (MQ=255)
                                                                                        ggCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAg                                                             <  1:1183561/150‑1 (MQ=255)
                                                                                          ccTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTg                                                           >  2:809536/1‑150 (MQ=255)
                                                                                           cTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTgg                                                          <  2:61922/150‑1 (MQ=255)
                                                                                               tCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGagag                                                      <  1:982163/150‑1 (MQ=255)
                                                                                                  cTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTat                                                   >  1:1487627/1‑150 (MQ=255)
                                                                                                  cTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTat                                                   >  2:230695/1‑150 (MQ=255)
                                                                                                  cTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTat                                                   >  1:901114/1‑150 (MQ=255)
                                                                                                    cAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTatga                                                 >  2:1039371/1‑150 (MQ=255)
                                                                                                     aaaTTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTatgat                                                >  1:919320/1‑150 (MQ=255)
                                                                                                             tGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAAt                                        <  2:105520/150‑1 (MQ=255)
                                                                                                                aGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTag                                     <  2:925625/150‑1 (MQ=255)
                                                                                                                aGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTag                                     <  1:1296992/150‑1 (MQ=255)
                                                                                                                aGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTag                                     <  2:63250/150‑1 (MQ=255)
                                                                                                                  aTTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTagaa                                   <  2:1521667/150‑1 (MQ=255)
                                                                                                                         aCAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTg                            >  2:895796/1‑150 (MQ=255)
                                                                                                                          cAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTgg                           >  2:877938/1‑150 (MQ=255)
                                                                                                                          cAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTgg                           >  2:300410/1‑150 (MQ=255)
                                                                                                                            gTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGgtg                         >  1:1319148/1‑150 (MQ=255)
                                                                                                                               ttCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTAc                      <  1:636377/150‑1 (MQ=255)
                                                                                                                                   ttGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGt                  <  2:1445184/150‑1 (MQ=255)
                                                                                                                                         aCAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAgcgcg            <  1:1441907/150‑1 (MQ=255)
                                                                                                                                          cAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAgcgcgc           >  2:68339/1‑150 (MQ=255)
                                                                                                                                          cAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAgcgcgc           >  2:1501404/1‑150 (MQ=255)
                                                                                                                                            atcatcATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCAt         <  1:865514/150‑1 (MQ=255)
                                                                                                                                              catcatTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCAttt       >  2:67587/1‑150 (MQ=255)
                                                                                                                                                tcatTCACCACGTTTATGATTCTCTCCATCGACAGCAACGCCGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTa     >  1:508966/1‑150 (MQ=255)
                                                                                                                                                 catTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTaa    <  1:1161842/150‑1 (MQ=255)
                                                                                                                                                 catTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTaa    >  2:1469465/1‑150 (MQ=255)
                                                                                                                                                   ttCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATg  <  1:230695/150‑1 (MQ=255)
                                                                                                                                                   ttCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATg  <  2:22868/150‑1 (MQ=255)
                                                                                                                                                   ttCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATg  <  1:623522/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGGCGTATGGAGCGCCTGGCGTCTGGTTAAAACGACCCTCAAGCAGCAACAGCTTCGCGGTTAACTTCCCTCTGGCCGGAGCCATTCCGGCCTTATCCCTCAAATTTTTTGAAGATTTTTGACAGTTTTCCTTGCTAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATG  >  NC_000913/3250313‑3250609

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