Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,223,594 A→T 5.1% K321I (AAA→ATA)  ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,223,5940AT 5.1% 246.2 / 2.9 79K321I (AAA→ATA) ebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base T (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

AGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATA  >  NC_000913/3223449‑3223737
                                                                                                                                                 |                                                                                                                                               
cgccccAGCTTTGCTATTCCTGCCGAACAGCCGCACCAATGGTCAGCAGAATCCCCATCTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAg                                                                                                                                               <  2:323573/144‑1 (MQ=255)
aGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAg                                                                                                                                               >  7:84598/1‑148 (MQ=255)
        ctttgctttTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGAc                                                                                                                                       <  8:279029/148‑1 (MQ=255)
          ttgctttTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACgg                                                                                                                                     >  7:344886/1‑148 (MQ=255)
          ttgctttTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGAAGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACACCGCGTTGGCTTCCGTGAGATCAAAGTGCGCGACgg                                                                                                                                     >  1:322097/1‑148 (MQ=255)
           tgctttTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACgg                                                                                                                                     <  8:344886/147‑1 (MQ=255)
               tttACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGt                                                                                                                                >  6:343820/1‑148 (MQ=255)
                     gTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGa                                                                                                                          >  8:378352/1‑148 (MQ=255)
                      tCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGa                                                                                                                          <  8:84598/147‑1 (MQ=255)
                             gCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTc                                                                                                                              >  4:242898/1‑136 (MQ=255)
                             gCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTc                                                                                                                              <  3:242898/136‑1 (MQ=255)
                                  agcaATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTa                                                                                                             >  8:12640/1‑148 (MQ=255)
                                  agcaATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTa                                                                                                             <  5:110344/148‑1 (MQ=255)
                                        ggTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGAt                                                                                                       >  7:352744/1‑148 (MQ=255)
                                         gTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGAt                                                                                                       <  1:171912/147‑1 (MQ=255)
                                           cagcagAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGAtgct                                                                                                    >  5:225393/1‑148 (MQ=255)
                                            agcagAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGAtgctg                                                                                                   <  6:146400/148‑1 (MQ=255)
                                             gcagAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGAtgctgc                                                                                                  <  7:24979/148‑1 (MQ=255)
                                              cagAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCa                                                                                                 >  1:172269/1‑148 (MQ=255)
                                                 aaTCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACTGTCTGTTCTGGATCAATAACCGTTATTTAATGCTGCACgg                                                                                              >  2:26251/1‑148 (MQ=255)
                                                      ccTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTc                                                                                          >  2:165934/1‑147 (MQ=255)
                                                       cTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCaa                                                                                        >  5:168013/1‑148 (MQ=255)
                                                         tATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAAcc                                                                                      >  5:358944/1‑148 (MQ=255)
                                                         tATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAAcc                                                                                      >  5:372266/1‑148 (MQ=255)
                                                          aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTAtgtg                                                                                                         <  5:212684/128‑1 (MQ=255)
                                                          aTCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTAtgtg                                                                                                         >  6:212684/1‑128 (MQ=255)
                                                             tttACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCa                                                                                  >  1:196973/1‑148 (MQ=255)
                                                                 ccATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTAtg                                                                                                           <  3:275866/119‑1 (MQ=255)
                                                                 ccATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTAtg                                                                                                           >  4:275866/1‑119 (MQ=255)
                                                                   aTCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTTATATCATAGTGCGAGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCCACCGTCACGAcca                                                                            >  1:312538/1‑146 (MQ=255)
                                                                   aTCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCATAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACaa                                                                            >  1:317598/1‑148 (MQ=255)
                                                                       ggTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGat                                                                        >  3:13596/1‑148 (MQ=255)
                                                                          cATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACATCGatcat                                                                     >  8:131915/1‑148 (MQ=255)
                                                                            tGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCg                                                                   >  8:145673/1‑148 (MQ=255)
                                                                             gACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCg                                                                   >  2:160497/1‑147 (MQ=255)
                                                                             gACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCg                                                                   >  2:236427/1‑147 (MQ=255)
                                                                             gACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCg                                                                   >  8:38717/1‑147 (MQ=255)
                                                                              aCGCTGAAAGACGCCCACGGCAACGTTCTGGAAGTGGTTCCACAACGCGTTGGCGTCCGTTATTTCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCa                                                                 <  5:259361/148‑1 (MQ=255)
                                                                               cGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCaa                                                                <  3:324608/148‑1 (MQ=255)
                                                                               cGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCaa                                                                <  2:172269/148‑1 (MQ=255)
                                                                                   gAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGc                                                            >  2:186057/1‑148 (MQ=255)
                                                                                     aaGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCcg                                                          >  2:35567/1‑148 (MQ=255)
                                                                                       gACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGCCAACGATCATCGCAAAGGCcgcg                                                        >  2:24010/1‑148 (MQ=255)
                                                                                          gCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCTCAAAGGCCGCGCCg                                                     >  8:196438/1‑148 (MQ=255)
                                                                                           ccAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGCTATCATAGTGCGCGACGGTCTGGTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGt                                                    <  2:317598/148‑1 (MQ=255)
                                                                                           ccAACGGCAACGGTCTGGAAGCGGTGCCACCACGCGTTGGCTTCCGTGATATCATAGTGCGCGACGGGCTGGGCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGAGCATCGCAAAGGCCGCGCCGt                                                    <  2:312538/148‑1 (MQ=255)
                                                                                            cAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTc                                                                                   >  2:253821/1‑116 (MQ=255)
                                                                                            cAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTc                                                                                   <  1:253821/116‑1 (MQ=255)
                                                                                             aaCGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTg                                                  >  8:362480/1‑148 (MQ=255)
                                                                                               cGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGa                                                <  7:38717/148‑1 (MQ=255)
                                                                                                    aCGTTCTGGAAGTGTTGCCACAACGCGTTGGCTTCCGTGATATCTAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGTCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGa                                           <  3:138314/148‑1 (MQ=255)
                                                                                                    aCGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGa                                           <  5:343820/148‑1 (MQ=255)
                                                                                                    aCGGTCTGGAAGGGGTGCCACAACGCGTTGGCGTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGa                                           <  2:310736/148‑1 (MQ=255)
                                                                                                       ttCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATcg                                        <  4:13596/148‑1 (MQ=255)
                                                                                                           ggAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTAGGAATGGATCGCGTc                                    >  4:235222/1‑148 (MQ=255)
                                                                                                                tggtgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGc                                                                                             >  2:329263/1‑86 (MQ=255)
                                                                                                                tggtgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGc                                                                                             <  1:329263/86‑1 (MQ=255)
                                                                                                                tggtgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGaa                               <  7:196438/148‑1 (MQ=255)
                                                                                                                tggtgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGaa                               <  7:249421/148‑1 (MQ=255)
                                                                                                                   tgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACTACAACGATCATCGCAAAGGCCGCGACGTTGGAATGGATCGCGTCgaga                                >  1:289376/1‑144 (MQ=255)
                                                                                                                   tgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGa                            <  1:35567/148‑1 (MQ=255)
                                                                                                                      cacaACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCCACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtct                         >  6:144574/1‑148 (MQ=255)
                                                                                                                      cacaACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtct                         <  6:372266/148‑1 (MQ=255)
                                                                                                                       acaACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGGCTGCTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGAtctc                        <  7:378352/148‑1 (MQ=255)
                                                                                                                        caACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTcc                       <  6:225393/148‑1 (MQ=255)
                                                                                                                        caACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTcc                       <  2:196973/148‑1 (MQ=255)
                                                                                                                          aCGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAg                     >  8:84538/1‑148 (MQ=255)
                                                                                                                          aCGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAg                     >  4:308149/1‑148 (MQ=255)
                                                                                                                          aCGCGTTGGCTTCCGTGATATCAAAGTGAGCGACGGTCTGTTCTGGATCAATAACCGTTTTGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAg                     <  6:358944/148‑1 (MQ=255)
                                                                                                                             cgTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTtg                  <  6:168013/148‑1 (MQ=255)
                                                                                                                                     tCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGaa                               >  1:385740/1‑127 (MQ=255)
                                                                                                                                     tCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGaa                               <  2:385740/127‑1 (MQ=255)
                                                                                                                                           atatCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGcaca     <  3:308149/147‑1 (MQ=255)
                                                                                                                                           atatCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGcaca     <  8:352744/147‑1 (MQ=255)
                                                                                                                                           atatCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACaa    <  7:362480/148‑1 (MQ=255)
                                                                                                                                           atatCAAAGTGCGCGACGGTCTGTTCTGGATAAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGACGTTGGAATGGATCGAGTCGAGAAAGATCTCCAGTTGATTAAGAAGcaca     >  6:216644/1‑147 (MQ=255)
                                                                                                                                            tatCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAat   <  1:236427/148‑1 (MQ=255)
                                                                                                                                            tatCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGAGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACaa    >  2:254183/1‑147 (MQ=255)
                                                                                                                                              tCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAata  <  1:186057/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATA  >  NC_000913/3223449‑3223737

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