Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,614,424 T→A 6.7% intergenic (+163/+257) yjjV → / ← yjjW putative DNase/putative pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,614,4240TA 6.7% 289.2 / 4.5 112intergenic (+163/+257)yjjV/yjjWputative DNase/putative pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base T (49/50);  new base A (3/5);  total (53/59)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATA  >  NC_000913/4614277‑4614572
                                                                                                                                                   |                                                                                                                                                    
cccgccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTc                                                                                                                                                     <  1:928415/147‑1 (MQ=255)
 aCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGAcg                                                                                                                                                    <  1:426925/149‑1 (MQ=21)
    ccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGt                                                                                                                                                 <  1:587627/149‑1 (MQ=255)
        aTCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTa                                                                                                                                             <  1:900483/149‑1 (MQ=255)
         tCCGGCAGTTGGCGCACAATGCCTGATGGGACGCTTGAGGCGTCTTATCGTGCCGACAGGGTTGTGCCGAACCGTAGGCCGGGTGAGGCGTTCACGCCGCATCCGGCAGTTGGCGCGCAATGCCTGATGCGACGGTTGTCGAGTCTTAt                                                                                                                                            >  2:1040340/1‑149 (MQ=14)
            ggCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                                                                         <  2:225552/149‑1 (MQ=255)
             gCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTTGGCGCAGAATGCCGGATGCGACGCGTGTCGCGTCTTATCATg                                                                                                                                        >  1:880549/1‑149 (MQ=255)
              cggTTGGCGCCCACTGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCCACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGc                                                                                                                                       <  2:721440/147‑1 (MQ=255)
              cAGTGGGCGCACAACGCCGGAGGCGACGCTCACCGCTTCTTATCATGCCCACAGGTTTGTGCCGAGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGc                                                                                                                                       <  2:589937/149‑1 (MQ=14)
                 ttGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTAACAGTAGGCCGCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGCCTTATCATGCCTa                                                                                                                                    <  1:924711/149‑1 (MQ=14)
                   ggCGCACACTGCGGGTTGCGACGCGTAACGCGTCTTATCCTGCCGATAGGTTTGGGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACa                                                                                                                                  <  2:142161/149‑1 (MQ=14)
                   ggCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTCATGCGACGCTTGTCGCGTCTTATCCTGCCTACa                                                                                                                                  <  1:989308/149‑1 (MQ=255)
                    gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACaa                                                                                                                                 >  1:1303995/1‑149 (MQ=255)
                      gcACAATGCCTGATGCGCCGCGTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTACGGGGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGt                                                                                                                               <  2:538182/149‑1 (MQ=14)
                        acaATGCCTGATGCGGCGCTTAACGCGTCTTATCATGCCTACCGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTACCATGCCTACAAGTCt                                                                                                                             <  2:753321/149‑1 (MQ=255)
                         cgaTGCCTGGTGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGAATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtg                                                                                                                            <  2:1008784/147‑1 (MQ=255)
                          gaggCCGGATGCGACGCTGAGCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGAGGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtgt                                                                                                                           <  2:1239864/146‑1 (MQ=14)
                          aaTGCCTGATGCGCCGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtgt                                                                                                                           <  2:785527/149‑1 (MQ=255)
                            tGCTTGATGTGAGGCTTAGCGCGTCTTATCATGGCTACAGGTTTGTGCCGAGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGc                                                                                                                         <  2:280275/149‑1 (MQ=14)
                            tGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACCGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGc                                                                                                                         <  2:1083388/149‑1 (MQ=255)
                            tGCCGGATGCGACGCTGAGCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGc                                                                                                                         <  2:1076055/149‑1 (MQ=255)
                            ggccgggtGCGACGCTTACCCGGTCTTTTCATGCCCACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGc                                                                                                                         <  2:196153/142‑1 (MQ=14)
                              ccTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCg                                                                                                                       >  2:440305/1‑149 (MQ=18)
                                tGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGaa                                                                                                                     >  1:13722/1‑149 (MQ=255)
                                    gCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGt                                                                                                                 <  1:1095611/149‑1 (MQ=255)
                                      gACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTCGGCCGGATCAGGCGTACACGCCGCATCCGGCGGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATAATGCCTACAAGTCTGTGCCGTACCGTAg                                                                                                               >  1:800210/1‑149 (MQ=14)
                                         gCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGcc                                                                                                            <  2:1225892/149‑1 (MQ=255)
                                         ccTTAACGCGTCTTATCATGCCTACAGGTTTTTGCCGAGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGcc                                                                                                            <  2:1014441/148‑1 (MQ=255)
                                          cTTAACGGGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCg                                                                                                           <  2:696532/149‑1 (MQ=255)
                                          cTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCg                                                                                                           <  2:177807/149‑1 (MQ=255)
                                              aCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATa                                                                                                       >  1:1108767/1‑149 (MQ=255)
                                              aCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTAGCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATa                                                                                                       >  2:398163/1‑149 (MQ=9)
                                              aCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGATGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCAGAACCGGAGGGCGGAtc                                                                                                       >  2:714596/1‑148 (MQ=255)
                                               cgcgTTTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATaa                                                                                                      <  2:1332072/149‑1 (MQ=255)
                                                  gTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAGTTGGCGCACAATGCCTGAAGCGACGCTTGTCGCGTCTTATCGTGCCTCCAAGTCTGTGCCGAACCGTAGGCCGGATAAGGc                                                                                                   >  2:635130/1‑149 (MQ=255)
                                                   tCTTATCATGCCTACAGGTTTGTGACGAGCCGTAGGCCGGATAAGGCGTTCGCGCCGAATCCGGCAGTTGGCGGACGATGCATGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCg                                                                                                  >  1:1011453/1‑149 (MQ=14)
                                                    cTTATCATGCCTTCAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGt                                                                                                 <  2:741855/149‑1 (MQ=255)
                                                     ttATCACGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGtt                                                                                                <  2:736368/149‑1 (MQ=255)
                                                       accATGCCTACAGGTTTGTGGCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCa                                                                                              <  1:133191/147‑1 (MQ=255)
                                                       aTCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCa                                                                                              <  2:941136/149‑1 (MQ=255)
                                                         cATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTTCCGACCCGTAGGCCGGATAAGGCGTTCAcg                                                                                            <  1:697467/149‑1 (MQ=255)
                                                                 cAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATCAGGCGTTCACGCCGCATcg                                                                                    >  1:129312/1‑148 (MQ=255)
                                                                   ggTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                  <  1:707951/149‑1 (MQ=255)
                                                                      ttGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTcgc                                                                                                                                                   >  1:15708/1‑81 (MQ=38)
                                                                       tgtgCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTGTCATGCCTACAAGTCTGTGCCGAACTGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                              <  1:876553/149‑1 (MQ=255)
                                                                         tgCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGGTAAGGCGTTCACGCCGCCTCCGGCAGTCg                                                                            >  1:200688/1‑149 (MQ=255)
                                                                         tgCCGAACCGTAGGCCGGATAAGGCGTTCACGCAGCATCCGGCAGTTGACGCACAATGCCTGATGCGACGCTTGTCACGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCCGTCg                                                                            >  1:865989/1‑149 (MQ=14)
                                                                           ccGCACCGTAGGCCGGATAAGGCGTTCACGCCGCCCCCGGCAGTGGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATGCGTCAggtcgc                                                                          <  2:13722/149‑6 (MQ=11)
                                                                               gccGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCAt                                                                      <  1:7521/148‑1 (MQ=255)
                                                                                 cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTCGGCGCaaaa                                                                    >  2:1158215/1‑149 (MQ=255)
                                                                                      gCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGGATCCGGCAGTCGGCGCCTAATGCCt                                                               >  1:1048907/1‑149 (MQ=255)
                                                                                       ccGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAGTCGGCGCATAATGGCTg                                                              >  2:601789/1‑149 (MQ=255)
                                                                                          gATAAGGCGTTCCCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATACGGCGTTCACGCCGCAGCGGGATCCCGGCGCATAATGCCtcgtg                                                           >  2:946143/1‑145 (MQ=14)
                                                                                             aaGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGgc                                                                          <  1:988434/131‑1 (MQ=18)
                                                                                               ggCGTTCCCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGCCGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACg                                                      <  1:901803/149‑1 (MQ=255)
                                                                                                gCGTTCACGGCGCTTCCGGCAGTTGGCGCACAATGCCTGATGCGCCGCTTGTCGCGTCTTATCATGCCTACAGGTCTGTGCCGAGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGc                                                     <  2:399880/149‑1 (MQ=14)
                                                                                                gCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGc                                                     >  1:1023226/1‑149 (MQ=255)
                                                                                                gCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCCGGCGGCGCATAATGCCTGGAGCGACGc                                                     >  2:95706/1‑149 (MQ=14)
                                                                                                gCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTCTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCCCCCCGCATCCGGCCGTCGCCGCATAATGCCTGACGcgccgc                                                     >  2:188016/1‑149 (MQ=14)
                                                                                                   ttCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTCGACGCCTACTGCCTGATGCGCCGCtcg                                                  >  2:927977/1‑147 (MQ=14)
                                                                                                   ttCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACCCCGCATCCGGCAGTCGGGGCATGAGGCCTGATGCGACGCttt                                                  >  2:204047/1‑148 (MQ=255)
                                                                                                   caagcgtcgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGCCTTATCATGCCTACAAGTCTGTGCCGAAGCGTCTGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTACGCGCACAATGCCTGGTGCGACGCGTg                                                  >  2:1151693/8‑149 (MQ=11)
                                                                                                      acgacgCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTcg                                               <  1:725061/149‑1 (MQ=255)
                                                                                                         ccgcATCCGGCAGTTGGCGCGCAATGCCCGATGCGACGCTTGCCGCGTCTTATCATGCCTACACGTCTGTGCCGAACCGTAGGCCGGATAAGGCGATCACGCCGCATCCGGCAGTCGGCTCATAATGCCTGATGCGACGCTTGTCGCGt                                            <  2:744048/149‑1 (MQ=14)
                                                                                                         ccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGt                                            <  1:581917/149‑1 (MQ=255)
                                                                                                         ccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGt                                            <  1:1243479/149‑1 (MQ=255)
                                                                                                         ccgcATCCGGCAGTGGGCGCACAATGCCTGATGCGACGCTTGGCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGt                                            <  1:286092/149‑1 (MQ=255)
                                                                                                               ccGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCa                                                                                                                                          >  1:266289/1‑49 (MQ=40)
                                                                                                               ccGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAAGTCTGAGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATc                                      <  2:982417/149‑1 (MQ=255)
                                                                                                                cggcggTTGGCACCCAAGGCCTGTTGCGACCCGTGACGCGTCTTAACCGGGCTATAGCGTTGTGTCGAACCGGGGGCCTTATTAGGCGTTCGCTCCTCATCCTGCAGTTGGCTCATAATGCCTGATGCGACGCGTGTCGCGTCTTATCa                                     <  2:656452/149‑1 (MQ=255)
                                                                                                                cGGCAGTTGGCGCACAATGCCTGATGCGACGCTTCTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCa                                     <  1:431406/149‑1 (MQ=255)
                                                                                                                 ggCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACCAGTCTGTGCCCAACCGTAGGCCGGATACGGCGTTCACGACGCATCCGGCGGTCGACGCAGAATGCCTGATCCGGCGCATGTCGCGTCTTAtgag                                    >  2:350929/1‑146 (MQ=14)
                                                                                                                 ggCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACGGTAGGCGGGATAAGGCGTTAACGCCGCATCCGGAAGTCGGCGCATAATGCCTGATGCGACGCTTGTGGCGTCTTATCAt                                    >  2:836843/1‑149 (MQ=14)
                                                                                                                 ggCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGCCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGGTCGGCGCTTACTGCCTGATGCGCCGCTCCTCGCGCCtcatcac                                    >  2:729465/1‑148 (MQ=14)
                                                                                                                  gCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATCCGGCGTTCACGCCGCATCCGGCCGGCGGCGCATAATGCCAGATGCGACGCTTGCCGCGACTGAGCATg                                   >  2:47839/1‑149 (MQ=14)
                                                                                                                  gCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATCAGGCGTTCCCGCCGCATCCGGCCGTCGGCGCACAATGCGTGATGCGACGCTTGTGGCGTCTGATCATg                                   >  2:647964/1‑149 (MQ=14)
                                                                                                                   cAGTTGGCGCACAATGCCTGATGCGACGCTTGGCGCGTCTTACCCTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGc                                  <  2:1299744/149‑1 (MQ=255)
                                                                                                                    aGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGcc                                 <  1:282442/149‑1 (MQ=255)
                                                                                                                    aGTTGGAGCACCATGCCTGATGCGACGCTTGCCGCGTCTTATCATGCCTACAGGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGcc                                 <  2:61661/149‑1 (MQ=255)
                                                                                                                     gTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCt                                >  1:241102/1‑149 (MQ=255)
                                                                                                                        ggCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCGTACCGACAAGACTGTGCCGAACCGTAGGCCGTATAAGGCGTTCCCGCCGCCTCCGGCCGTCGGCGCATAATGCCTGCTGCGACGCGTGGCGCGTCCTATCCTGCCcaca                             >  2:1161878/1‑149 (MQ=14)
                                                                                                                        ggCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCCTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACa                             >  2:1062436/1‑149 (MQ=255)
                                                                                                                         gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGCCTTATCAGGCCTACaa                            >  2:239740/1‑148 (MQ=255)
                                                                                                                         gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCCTGCCTACAg                            >  2:201532/1‑149 (MQ=255)
                                                                                                                         gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAg                            <  2:137199/149‑1 (MQ=255)
                                                                                                                         gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGGCTTATCATGCGTACAg                            >  1:346729/1‑149 (MQ=255)
                                                                                                                           gcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGt                          >  2:539374/1‑149 (MQ=255)
                                                                                                                           gcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCAGAATGCCTGATGCGACGCTTGTCGCGGCTTATCATGCCTACCggc                          >  2:1330761/1‑148 (MQ=255)
                                                                                                                             acaATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTa                                                                                                                >  2:373744/1‑61 (MQ=255)
                                                                                                                             acaATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGCCGCTTGTCGCGCCTTATCATGCCTACAGCttt                        >  1:1105937/1‑149 (MQ=255)
                                                                                                                                  gCCTGACGCGCCGCTTGTCGCCTCTAAACATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGcc                   <  1:422015/149‑1 (MQ=17)
                                                                                                                                     tGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCGTATCATGCCTACAGGTTTGTGCCGac                >  1:1200507/1‑148 (MQ=255)
                                                                                                                                          cGACGCTTGTCGCGTTTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTa           <  1:165226/149‑1 (MQ=255)
                                                                                                                                          cGACCCTTGTGGCGTTTTATCATGCTTATAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTa           <  1:650622/149‑1 (MQ=17)
                                                                                                                                           gACGCTTGGCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGCCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGGAg          >  2:834574/1‑149 (MQ=255)
                                                                                                                                              gCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGcc       >  2:1193029/1‑149 (MQ=255)
                                                                                                                                               cTTGTCGCGTCTTATCGTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCg      >  2:593688/1‑149 (MQ=255)
                                                                                                                                                ttGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCgg     >  1:476699/1‑149 (MQ=255)
                                                                                                                                                ttGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCgg     >  2:489250/1‑149 (MQ=255)
                                                                                                                                                   tCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGAGGCGGGCCTTGTCGCGCCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATa  >  2:1157250/1‑149 (MQ=14)
                                                                                                                                                   |                                                                                                                                                    
CACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATA  >  NC_000913/4614277‑4614572

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