Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,156,913 G→A G46D (GGT→GAT)  ycfH → putative DNase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,156,9130GA96.3% 1338.1 / 26.1 402G46D (GGT→GAT) ycfHputative DNase
Reads supporting (aligned to +/- strand):  ref base G (5/10);  new base A (185/202);  total (190/212)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 409 total aligned reads displayed.

AAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGG  >  NC_000913/1156765‑1157063
                                                                                                                                                    |                                                                                                                                                      
aacagagaCCTCATGTTTTTAGTCGAATCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGACCGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGAtt                                                                                                                                                      <  1:519664/148‑1 (MQ=255)
aaGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGAtt                                                                                                                                                      >  8:1290748/1‑151 (MQ=255)
 agagagACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGACACAACATTAACGGAttt                                                                                                                                                     >  4:262247/1‑150 (MQ=255)
    gagaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTAttt                                                                                                                                                  <  6:997543/151‑1 (MQ=255)
     agaCATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTa                                                                                                                                                 <  5:133563/151‑1 (MQ=255)
        catcatGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACat                                                                                                                                              >  4:1022533/1‑151 (MQ=255)
        catcatGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACat                                                                                                                                              >  8:276459/1‑151 (MQ=255)
           catGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGGACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATg                                                                                                                                           >  8:374405/1‑151 (MQ=255)
              gTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCg                                                                                                                                         <  1:145051/150‑1 (MQ=255)
               tttttAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCggg                                                                                                                                       >  8:1204009/1‑151 (MQ=255)
                ttttAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCggg                                                                                                                                       >  7:1623663/1‑150 (MQ=255)
                   tAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATcg                                                                                                                                   >  5:69425/1‑150 (MQ=255)
                   tAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATc                                                                                                                                    <  4:673594/150‑1 (MQ=255)
                   tAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATc                                                                                                                                    <  1:681327/150‑1 (MQ=255)
                    aGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCt                                                                                                                                   <  6:450818/150‑1 (MQ=255)
                       cGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTa                                                                                                                               <  5:898347/151‑1 (MQ=255)
                        gACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAg                                                                                                                              <  1:889871/151‑1 (MQ=255)
                        gACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAg                                                                                                                              <  1:396015/151‑1 (MQ=255)
                         aCTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAgg                                                                                                                             >  7:56145/1‑151 (MQ=255)
                          cTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGc                                                                                                                            >  5:976383/1‑151 (MQ=255)
                          cTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGc                                                                                                                            >  7:1280141/1‑151 (MQ=255)
                           tCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCg                                                                                                                           >  1:820551/1‑151 (MQ=255)
                            cacaCTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCg                                                                                                                           <  4:1117761/150‑1 (MQ=255)
                               acTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACg                                                                                                                       <  4:347120/151‑1 (MQ=255)
                               acTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACg                                                                                                                       <  8:1569157/151‑1 (MQ=255)
                                 tGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTg                                                                                                                     >  5:179490/1‑151 (MQ=255)
                                  gCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGt                                                                                                                      <  2:412808/149‑1 (MQ=255)
                                  gCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTg                                                                                                                     <  4:1208479/150‑1 (MQ=255)
                                        tcGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgt                                                                                                              >  6:81817/1‑151 (MQ=255)
                                          gATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAACGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtcg                                                                                                            <  1:1267372/151‑1 (MQ=255)
                                          gATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtcg                                                                                                            >  3:159479/1‑151 (MQ=255)
                                            tGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAAcgtcgt                                                                                                           >  8:844427/1‑150 (MQ=255)
                                              gTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAtt                                                                                                        >  8:597825/1‑151 (MQ=255)
                                               gcTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAttt                                                                                                       <  1:198023/150‑1 (MQ=255)
                                                cTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAtttt                                                                                                      >  3:374851/1‑151 (MQ=255)
                                                cTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAtttt                                                                                                      >  8:432517/1‑151 (MQ=255)
                                                  ggATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAtttt                                                                                                      <  4:616323/149‑1 (MQ=255)
                                                  ggATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTAttt                                                                                                       <  1:1107153/148‑1 (MQ=255)
                                                   gATTATCAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACCCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCTGATTATTTACATATGCGCGATCTGGTACGAGAACGTGACAACGTCGTATTTTCtt                                                                                                   >  2:650871/1‑151 (MQ=255)
                                                      tATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTtgt                                                                                                 <  8:1110089/150‑1 (MQ=255)
                                                      tATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTtgt                                                                                                 <  3:1515170/150‑1 (MQ=255)
                                                       aTGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTtgtg                                                                                                <  5:1175947/150‑1 (MQ=255)
                                                       aTGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTtgtg                                                                                                <  4:82801/150‑1 (MQ=255)
                                                        tGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGc                                                                                              >  5:103635/1‑151 (MQ=255)
                                                         gAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCtt                                                                                                   <  4:1603757/145‑1 (MQ=255)
                                                         gAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGc                                                                                              <  1:504760/150‑1 (MQ=255)
                                                          aaTCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACCTGACAACGTCGTATTTTCTTGTGGCGt                                                                                            >  5:425118/1‑151 (MQ=255)
                                                                tGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATcc                                                                                      >  8:679515/1‑151 (MQ=255)
                                                                 gCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATcc                                                                                      <  6:1144475/150‑1 (MQ=255)
                                                                 gCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCa                                                                                     <  4:752185/151‑1 (MQ=255)
                                                                     aaGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCa                                                                                 <  2:116426/151‑1 (MQ=255)
                                                                      aGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCaa                                                                                <  6:1329905/151‑1 (MQ=255)
                                                                          cGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAAcc                                                                              >  6:1586283/1‑149 (MQ=255)
                                                                             ggATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGaa                                                                          <  3:874582/150‑1 (MQ=255)
                                                                             ggATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAAt                                                                         >  2:1427569/1‑151 (MQ=255)
                                                                               aTGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATg                                                                        >  8:1026554/1‑150 (MQ=255)
                                                                                tGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGAt                                                                      >  5:1481439/1‑151 (MQ=255)
                                                                                 gACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCACCCAGAATGATc                                                                     >  7:777478/1‑151 (MQ=255)
                                                                                  aCGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATcc                                                                    <  6:1292143/151‑1 (MQ=255)
                                                                                   cGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTCTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATccc                                                                   <  4:982795/151‑1 (MQ=255)
                                                                                    gTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCt                                                                  >  7:413665/1‑151 (MQ=255)
                                                                                         ggCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGa                                                              <  2:1051856/150‑1 (MQ=255)
                                                                                             aaaGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTag                                                         <  3:613125/151‑1 (MQ=255)
                                                                                               aGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTaga                                                        >  7:849813/1‑150 (MQ=255)
                                                                                                  cgccgcACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTaga                                                        >  4:1292819/1‑147 (MQ=255)
                                                                                                       cACGCGATGTGAAATTTTGTCTGNCAGTCGNNACAACATTACNNNANTATTTACATATGNGGGATNTNGTAGGNGAACGTGACAACGNNGTATTTTCTNGTGGNGTNCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACg                                               >  2:830006/1‑151 (MQ=255)
                                                                                                       cACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACg                                               <  5:74059/151‑1 (MQ=255)
                                                                                                        aCGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGc                                              >  1:1036745/1‑151 (MQ=255)
                                                                                                           cgATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGt                                           <  8:1673428/151‑1 (MQ=255)
                                                                                                             aTGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCt                                         <  5:1599712/151‑1 (MQ=255)
                                                                                                               gtgAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTgg                                       >  4:200461/1‑151 (MQ=255)
                                                                                                               gtgAAATTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCAACCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTg                                        <  3:1029598/150‑1 (MQ=255)
                                                                                                                   aaTTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTggcggc                                   <  1:1166657/151‑1 (MQ=255)
                                                                                                                    aTTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTggcggc                                   >  7:932594/1‑150 (MQ=255)
                                                                                                                    aTTTTGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCa                                  <  3:1319938/151‑1 (MQ=255)
                                                                                                                       ttGTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCAACGATGTAGAAGATTTACGCCGTCTGGCGGCaga                                >  3:422359/1‑150 (MQ=255)
                                                                                                                         gTCTGGCAGTCGCTACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCagaaga                             <  5:577198/151‑1 (MQ=255)
                                                                                                                         gTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCagaaga                             <  3:941014/151‑1 (MQ=255)
                                                                                                                         gTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCagaag                              <  1:1522162/150‑1 (MQ=255)
                                                                                                                         gTCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCaga                                >  4:479912/1‑148 (MQ=255)
                                                                                                                          tCTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCagaag                              >  4:852497/1‑149 (MQ=255)
                                                                                                                          gcTGGCCGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAagag                            <  7:1455180/150‑1 (MQ=255)
                                                                                                                           cTGGCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAgg                           >  2:1591973/1‑151 (MQ=255)
                                                                                                                             ggCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGgt                         <  6:373791/151‑1 (MQ=255)
                                                                                                                              gCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGgtg                        <  3:1123580/151‑1 (MQ=255)
                                                                                                                              gCAGTCGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATATACGCCGTCTGGCGGCAGAAGAGGgtg                        <  2:1325861/151‑1 (MQ=255)
                                                                                                                                   cGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGtgttgt                    >  3:1458123/1‑150 (MQ=255)
                                                                                                                                   cGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGtgttgt                    <  3:1417902/150‑1 (MQ=255)
                                                                                                                                   cGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGtgttgt                    >  6:1661973/1‑150 (MQ=255)
                                                                                                                                   cGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTa                   <  6:657831/151‑1 (MQ=255)
                                                                                                                                   cGCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTa                   <  4:1446367/151‑1 (MQ=255)
                                                                                                                                    gCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGtgttgt                    >  8:1381243/1‑149 (MQ=255)
                                                                                                                                    gCCACAACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAg                  <  7:1532390/151‑1 (MQ=255)
                                                                                                                                      cacaACATTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAgcg                >  2:294973/1‑151 (MQ=255)
                                                                                                                                        caacaTTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCt              <  1:1664560/151‑1 (MQ=255)
                                                                                                                                            aTTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGt          >  8:1517710/1‑151 (MQ=255)
                                                                                                                                            aTTACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGt          <  4:1330685/151‑1 (MQ=255)
                                                                                                                                              tACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTg         <  2:1499117/150‑1 (MQ=255)
                                                                                                                                              tACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTg         >  3:134062/1‑150 (MQ=255)
                                                                                                                                              tACCGGATTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGACGTCTGGCGGCAGAAGACGCTGTTGTAGCGCTGGGTg         >  7:1047964/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AAGAGAGACATCATGTTTTTAGTCGACTCACACTGCCATCTCGATGGTCTGGATTATGAATCTTTGCATAAGGACGTGGATGACGTTCTGGCGAAAGCCGCCGCACGCGATGTGAAATTTTGTCTGGCAGTCGCCACAACATTACCGGGTTATTTACATATGCGGGATCTGGTAGGCGAACGTGACAACGTCGTATTTTCTTGTGGCGTGCATCCACTCAACCAGAATGATCCCTACGATGTAGAAGATTTACGCCGTCTGGCGGCAGAAGAGGGTGTTGTAGCGCTGGGTGAAACCGG  >  NC_000913/1156765‑1157063

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