Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,580,619 T→C 5.6% K282E (AAA→GAA)  hsdS ← specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,580,6190TC 5.6% 227.2 / 3.7 71K282E (AAA→GAA) hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base T (30/37);  new base C (2/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCA  >  NC_000913/4580471‑4580764
                                                                                                                                                    |                                                                                                                                                 
ggTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAAttt                                                                                                                                                   >  2:20911/1‑149 (MQ=255)
        aTCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGt                                                                                                                                           >  3:291325/1‑149 (MQ=255)
            ttGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAg                                                                                                                                       <  2:132093/149‑1 (MQ=255)
             tGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATGTGTGGCGGTTTAGt                                                                                                                                      >  3:435443/1‑149 (MQ=255)
              ggTTAATCGAGCTCGAATAAGTTTATCAGGATATAGTAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGtt                                                                                                                                     <  3:209590/149‑1 (MQ=255)
                  aaTCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACt                                                                                                                                 <  3:150134/149‑1 (MQ=255)
                    tCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTACAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACttt                                                                                                                               <  1:460163/149‑1 (MQ=255)
                      gAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCt                                                                                                                             <  2:554819/149‑1 (MQ=255)
                        gCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAATCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGa                                                                                                                           <  2:299338/149‑1 (MQ=255)
                           cGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACa                                                                                                                        >  1:71423/1‑149 (MQ=255)
                             aaTAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACAtt                                                                                                                      >  1:140975/1‑149 (MQ=255)
                             aaTAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACAtt                                                                                                                      >  1:342245/1‑149 (MQ=255)
                              aTAAGTTTATCAGGATATAGCAAATTTTGATGCTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTCCTGAACATTc                                                                                                                     >  4:485723/1‑149 (MQ=255)
                                aaGTTTATCAGGATATAGCAAATTTTGATGCTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTa                                                                                                                   >  4:437955/1‑149 (MQ=255)
                                  gTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTGGTTCACTTTCTGAACATTCTAGa                                                                                                                 >  1:287520/1‑149 (MQ=255)
                                         aGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGaa                                                                                                          <  3:472355/149‑1 (MQ=255)
                                                 gCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGttt                                                                                                  <  4:97930/149‑1 (MQ=255)
                                                 gCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGttt                                                                                                  <  4:208645/149‑1 (MQ=255)
                                                  cAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGtttt                                                                                                 <  2:474826/149‑1 (MQ=255)
                                                       tttGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCt                                                                                            <  3:473779/149‑1 (MQ=255)
                                                         tGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTAc                                                                                          <  4:278986/149‑1 (MQ=255)
                                                         tGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTAc                                                                                          <  2:270895/149‑1 (MQ=255)
                                                          gATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACa                                                                                         <  3:393179/149‑1 (MQ=255)
                                                             gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATgg                                                                                      >  1:431088/1‑149 (MQ=255)
                                                              ttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCTCTTTCTGAACATTTTAGAAACCGAATATCGTTTTGATCTACATGGc                                                                                     >  3:210280/1‑149 (MQ=255)
                                                              ttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGc                                                                                     <  2:46206/149‑1 (MQ=255)
                                                              ttgcAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGc                                                                                     <  1:224857/149‑1 (MQ=255)
                                                                   aTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTAATCTACATGGCCAGCa                                                                                >  1:232422/1‑149 (MQ=255)
                                                                           aaTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGa                                                                        >  3:311408/1‑149 (MQ=255)
                                                                           aaTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGa                                                                        >  1:451372/1‑149 (MQ=255)
                                                                             tAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAAc                                                                      >  3:368860/1‑149 (MQ=255)
                                                                                cccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTaa                                                                   <  3:485723/149‑1 (MQ=255)
                                                                                   acaAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGc                                                                <  2:140975/149‑1 (MQ=255)
                                                                                      aaCACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgta                                                             <  2:291774/149‑1 (MQ=255)
                                                                                           caacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTAtt                                                        <  3:325495/149‑1 (MQ=255)
                                                                                           caacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTAtt                                                        >  3:195536/1‑149 (MQ=255)
                                                                                           caacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTAtt                                                        >  1:234578/1‑149 (MQ=255)
                                                                                            aacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTg                                                       >  1:85946/1‑149 (MQ=255)
                                                                                             acaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTgg                                                      <  3:274227/149‑1 (MQ=255)
                                                                                               aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTc                                                                                                                     >  3:391127/1‑84 (MQ=255)
                                                                                               aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTc                                                                                                                     <  4:391127/84‑1 (MQ=255)
                                                                                                   tCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGAcc                                                >  3:136347/1‑149 (MQ=255)
                                                                                                   tCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTCGTGGCGGTTTAGTTCGCTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGAcc                                                >  4:326795/1‑149 (MQ=255)
                                                                                                   tCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTCGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGAcc                                                >  2:255668/1‑149 (MQ=255)
                                                                                                        aCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAacac                                           <  1:20911/149‑1 (MQ=255)
                                                                                                          ttCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAAcacca                                         >  2:189426/1‑149 (MQ=255)
                                                                                                          ttCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAAcacca                                         <  4:392213/149‑1 (MQ=255)
                                                                                                           tCCGATATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGGCAGCACTTACAGAACTAATGCGTAGTATTGGATGACCAAcaccac                                        >  4:488984/1‑149 (MQ=255)
                                                                                                             cGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACtt                                      <  4:136347/149‑1 (MQ=255)
                                                                                                                  taGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCAtt                                 >  1:444270/1‑149 (MQ=255)
                                                                                                                  taGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCAtt                                 >  4:209403/1‑149 (MQ=255)
                                                                                                                       aGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCt                            <  2:234578/149‑1 (MQ=255)
                                                                                                                       aGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACATACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCt                            >  1:356680/1‑149 (MQ=255)
                                                                                                                        gTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCtt                           <  2:71423/149‑1 (MQ=255)
                                                                                                                        gTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCtt                           >  4:36737/1‑149 (MQ=255)
                                                                                                                           cttaaaAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTtga                        <  3:202750/147‑1 (MQ=255)
                                                                                                                           aataaaAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTtga                        <  2:342245/149‑1 (MQ=255)
                                                                                                                           aataaaAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTtga                        >  1:324129/1‑149 (MQ=255)
                                                                                                                              aaaaGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGGATTGGATGACCAACACCACTTTCATTTGGCCTttatag                     >  4:202750/1‑147 (MQ=255)
                                                                                                                               aaaGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACAGCACTTTCATTTGGCTTTGATGaa                    >  3:100748/1‑149 (MQ=255)
                                                                                                                                aaGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGaaa                   <  2:421997/149‑1 (MQ=255)
                                                                                                                                aaGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGaaa                   <  4:105088/149‑1 (MQ=255)
                                                                                                                                aaGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGaaa                   >  2:446326/1‑149 (MQ=255)
                                                                                                                                aaGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGaaa                   <  1:446326/149‑1 (MQ=255)
                                                                                                                                aaGATCTCCATCTTGTAATTCGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGaaa                   <  1:255668/149‑1 (MQ=255)
                                                                                                                                aaGATCTCCATCTTGTAATTCGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGaaa                   <  3:326795/149‑1 (MQ=255)
                                                                                                                                   aTCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGAc                <  1:3966/149‑1 (MQ=255)
                                                                                                                                    tctcCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGAcc               <  3:36737/149‑1 (MQ=255)
                                                                                                                                     ctcCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCa              <  2:437511/149‑1 (MQ=255)
                                                                                                                                            ttGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTa       >  2:372589/1‑149 (MQ=255)
                                                                                                                                              gcaaTTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAAt     <  3:191085/147‑1 (MQ=255)
                                                                                                                                                 aTTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCa  <  3:503349/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCA  >  NC_000913/4580471‑4580764

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