Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,881,358 T→A 6.2% D151V (GAC→GTC)  dmlR ← transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,881,3580TA 6.2% 491.1 / 13.3 140D151V (GAC→GTC) dmlRtranscriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base T (56/75);  new base A (5/4);  total (61/79)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

TGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAG  >  NC_000913/1881211‑1881504
                                                                                                                                                   |                                                                                                                                                  
tgatCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTc                                                                                                                                                   >  3:568779/1‑149 (MQ=255)
 gatCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACa                                                                                                                                                  >  3:453251/1‑149 (MQ=255)
  atCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCAt                                                                                                                                                 >  4:5385/1‑149 (MQ=255)
   tCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATc                                                                                                                                                <  1:973128/149‑1 (MQ=255)
   tCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATc                                                                                                                                                >  1:472822/1‑149 (MQ=255)
   tCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGGTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATc                                                                                                                                                <  4:825752/149‑1 (MQ=255)
     gcgcTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGc                                                                                                                                              >  1:571773/1‑149 (MQ=255)
      cgcTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcc                                                                                                                                             <  2:401078/149‑1 (MQ=255)
       gcTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCa                                                                                                                                            >  1:1246257/1‑149 (MQ=255)
       gcTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCa                                                                                                                                            >  1:662339/1‑149 (MQ=255)
        ctctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccc                                                                                                                                           <  2:246543/149‑3 (MQ=255)
         tctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAAt                                                                                                                                          <  2:415915/149‑1 (MQ=255)
          ctTTAATGATCAGGCACGGCAATGCAGATAGATCGGTTAAATGTTTCCGTGCGCCATGTTGCGCAATAAATTCCGTCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAATTTAGGGGCGATGTCATCGCCAATg                                                                                                                                         <  3:914958/149‑1 (MQ=255)
          ctTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATg                                                                                                                                         <  2:774082/149‑1 (MQ=255)
          ctTTAATGACCAGGCACGGCAATGCAGAAAGAACGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTCAATTAGGGGCGATGACATCGccccta                                                                                                                                         <  4:453251/149‑5 (MQ=255)
           tttAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGc                                                                                                                                        >  2:1042466/1‑149 (MQ=255)
            ttAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGACATCGcccctaat                                                                                                                                       <  2:571773/149‑7 (MQ=255)
              aaTGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGaa                                                                                                                                     <  3:1033671/149‑1 (MQ=255)
                 gACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAatat                                                                                                                                  >  1:790871/1‑149 (MQ=255)
                   ccAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCg                                                                                                                                <  2:42491/149‑1 (MQ=255)
                    cAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCga                                                                                                                               >  1:236842/1‑149 (MQ=255)
                     aGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCgag                                                                                                                              >  2:720390/1‑149 (MQ=255)
                     aGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCgag                                                                                                                              <  2:269344/149‑1 (MQ=255)
                     aGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATGAGGGGCGATGTCATCGCCAATGCGAATATCgag                                                                                                                              >  2:720093/1‑149 (MQ=255)
                      ggCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCgaga                                                                                                                             >  1:43662/1‑149 (MQ=255)
                       gCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGAt                                                                                                                            >  3:565890/1‑149 (MQ=255)
                         aCGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCg                                                                                                                          >  2:1031919/1‑149 (MQ=255)
                         aCGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCg                                                                                                                          >  4:309309/1‑149 (MQ=255)
                            gCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGaca                                                                                                                       >  4:654546/1‑149 (MQ=255)
                             cAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGacac                                                                                                                      <  1:1096536/149‑1 (MQ=255)
                             cAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGacac                                                                                                                      <  2:1474320/149‑1 (MQ=255)
                                tGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCtt                                                                                                                   <  2:517558/149‑1 (MQ=255)
                                 gCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTc                                                                                                                  <  4:298851/149‑1 (MQ=255)
                                    gAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACAcc                                                                                                                     >  1:391222/1‑143 (MQ=255)
                                    gAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCAtt                                                                                                               <  2:176090/149‑1 (MQ=255)
                                    gAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCAtt                                                                                                               <  3:792591/149‑1 (MQ=255)
                                       aGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGAc                                                                                                            <  3:531542/149‑1 (MQ=255)
                                        gATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGAcc                                                                                                           <  1:480332/149‑1 (MQ=255)
                                           cGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCaaa                                                                                                        <  2:822834/149‑1 (MQ=255)
                                             gTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCaaa                                                                                                        >  4:827999/1‑147 (MQ=255)
                                             gTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATc                                                                                                      <  4:234725/149‑1 (MQ=255)
                                             gTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATc                                                                                                      <  4:1308407/149‑1 (MQ=255)
                                             gTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATc                                                                                                      <  3:332190/149‑1 (MQ=255)
                                             gTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTGCATTGACCaaa                                                                                                        <  3:827999/147‑1 (MQ=255)
                                              ttAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGGCAATGCGAATTTCGAGATCGACAACTTCATTGACCAAATCGa                                                                                                    >  2:318875/1‑149 (MQ=255)
                                                     tttCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTc                                                                                              >  2:671913/1‑149 (MQ=255)
                                                       tCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCgg                                                                                            <  3:1367179/149‑1 (MQ=255)
                                                            gcgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTc                                                                                       <  4:1373198/149‑1 (MQ=255)
                                                             cgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGAGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTcc                                                                                      <  1:1083722/149‑1 (MQ=255)
                                                              gcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCca                                                                                     >  3:1208974/1‑149 (MQ=255)
                                                                cATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCga                                                                                                                               >  1:1091680/1‑105 (MQ=255)
                                                                cATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCga                                                                                                                               <  2:1091680/105‑1 (MQ=255)
                                                                cATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCAtt                                                                                                               >  4:615319/1‑121 (MQ=255)
                                                                cATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCAtt                                                                                                               <  3:615319/121‑1 (MQ=255)
                                                                cATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCcaca                                                                                   <  3:1357840/149‑1 (MQ=255)
                                                                  tGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATc                                                                                 <  1:938498/149‑1 (MQ=255)
                                                                       cgcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAAc                                                                            <  3:217957/149‑1 (MQ=255)
                                                                        gcAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACg                                                                           >  4:1093287/1‑149 (MQ=255)
                                                                         cAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGt                                                                          >  3:33745/1‑149 (MQ=255)
                                                                          aaTAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTa                                                                         >  2:186352/1‑149 (MQ=255)
                                                                          aaTAAATTCCGACGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTa                                                                         >  4:1219689/1‑149 (MQ=255)
                                                                             aaaTTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAAtt                                                                      >  2:511000/1‑149 (MQ=255)
                                                                                ttCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCa                                                                   <  2:649873/149‑1 (MQ=255)
                                                                                   cGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGtt                                                                <  2:162149/149‑1 (MQ=255)
                                                                                       gACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGtt                                                                >  4:1175857/1‑145 (MQ=255)
                                                                                       gACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGtt                                                                <  3:1175857/145‑1 (MQ=255)
                                                                                        aCGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGa                                                           <  2:43662/149‑1 (MQ=255)
                                                                                             cAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGc                                                      <  1:720093/149‑1 (MQ=255)
                                                                                                 gaATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTc                                                  <  2:73552/149‑1 (MQ=255)
                                                                                                  aaTACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCg                                                 >  1:573479/1‑149 (MQ=255)
                                                                                                        gTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAgcg                                           >  3:993869/1‑149 (MQ=255)
                                                                                                         tAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAgcgc                                          <  1:115453/149‑1 (MQ=255)
                                                                                                          aaTTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCt                                         <  4:273867/149‑1 (MQ=255)
                                                                                                           aTTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCta                                        <  1:1163014/149‑1 (MQ=255)
                                                                                                             tGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCTATGTCATCGCCAATGCGAATATCGAGCTCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCtaat                                      <  1:862722/149‑1 (MQ=255)
                                                                                                                ttgCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTc                                                                                       >  3:1322554/3‑97 (MQ=255)
                                                                                                                tCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTc                                                                                       <  4:1322554/97‑1 (MQ=255)
                                                                                                                tCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAAc                                   <  2:1234213/149‑1 (MQ=255)
                                                                                                                tCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAAc                                   >  1:1240346/1‑149 (MQ=255)
                                                                                                                      gTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGa                             <  2:291545/149‑1 (MQ=255)
                                                                                                                      gTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGa                             <  3:879149/149‑1 (MQ=255)
                                                                                                                           cGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGccac                        <  3:5385/149‑1 (MQ=255)
                                                                                                                           cGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGccac                        <  3:309309/149‑1 (MQ=255)
                                                                                                                                gATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCt                   <  4:1157578/149‑1 (MQ=255)
                                                                                                                                 aTTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTaa                                                                        <  4:1128412/95‑1 (MQ=255)
                                                                                                                                 aTTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTaa                                                                        >  3:1128412/1‑95 (MQ=255)
                                                                                                                                  ttAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgt                 >  2:354063/1‑149 (MQ=255)
                                                                                                                                  ttAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgt                 >  2:122805/1‑149 (MQ=255)
                                                                                                                                  ttAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTg                  <  4:1108841/148‑1 (MQ=255)
                                                                                                                                  ttAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTg                  >  3:1108841/1‑148 (MQ=255)
                                                                                                                                   tAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgtc                <  2:460240/149‑1 (MQ=255)
                                                                                                                                    aaaTTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgtcg               <  1:921620/149‑1 (MQ=255)
                                                                                                                                       ttAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTcc            >  1:1411869/1‑149 (MQ=255)
                                                                                                                                       ttAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTcc            >  4:109286/1‑149 (MQ=255)
                                                                                                                                         attggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCaa          >  4:209060/4‑149 (MQ=255)
                                                                                                                                         attggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCaa          >  4:694818/4‑149 (MQ=255)
                                                                                                                                          ttggCGATGACATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCaaa         >  1:231740/3‑149 (MQ=255)
                                                                                                                                          ggggCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCaccacc                    <  2:181800/138‑1 (MQ=255)
                                                                                                                                                aTGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCaaa   <  1:318875/149‑1 (MQ=255)
                                                                                                                                                aTGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCaaa   >  1:71131/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGATCGCGCTCTTTAATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAG  >  NC_000913/1881211‑1881504

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: