Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 848,773 T→C 5.5% T47A (ACC→GCC)  dps ← Fe‑binding and storage protein; stress‑inducible DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913848,7730TC 5.5% 226.8 / 2.9 73T47A (ACC→GCC) dpsFe‑binding and storage protein; stress‑inducible DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base C (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

AGTGGTCCCCAGAGCTACACCGCCCAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAGCGGTACTCATAATTTCAT  >  NC_000913/848627‑848919
                                                                                                                                                  |                                                                                                                                                  
aGTGGTCCCCAGAGCTACACCGCCCAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTaa                                                                                                                                                  <  2:205248/149‑1 (MQ=255)
     tCCCCAGAGCTACACCGCCCAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCa                                                                                                                                               >  3:176896/1‑147 (MQ=255)
       cccAGAGCTACACCGCCCAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCaaag                                                                                                                                            >  5:342836/1‑148 (MQ=255)
         cAGAGCTACACCGCCCAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCCTCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGGGCGCTTGTTTGGTAATCaaaga                                                                                                                                           >  1:122307/1‑147 (MQ=255)
            agCTACACCGCCCAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGAt                                                                                                                                      <  2:199234/149‑1 (MQ=255)
                 cacCGCCCAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAAt                                                                                                                                  <  8:331278/148‑1 (MQ=255)
                       ccAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTgg                                                                                                                           <  2:170140/149‑1 (MQ=255)
                        cAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATc                                                                                                                                     >  4:5445/1‑138 (MQ=255)
                        cAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATc                                                                                                                                     <  3:5445/138‑1 (MQ=255)
                               actgcacGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTg                                                                                                                   <  4:98394/149‑1 (MQ=255)
                                ctgcacGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTgg                                                                                                                  >  2:403737/1‑149 (MQ=255)
                                   cacGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGa                                                                                                               >  1:230516/1‑149 (MQ=255)
                                   cacGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGa                                                                                                               >  4:170482/1‑149 (MQ=255)
                                      tttCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGACATACAGCATTTCATGTACGGCTATGAAGTTAGCGCCGCGCATGTTGCAGTACGCATGTTTGGTAATCAAAGAAAGATCAATTCACAGGATAACCTGTCGttt                                                                                                             >  6:18/2‑148 (MQ=255)
                                            ccATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTACAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGAGATTCAGCAAc                                                                                                      >  7:391168/1‑149 (MQ=255)
                                             cATTGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGGCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAgaatct                                                                                                     <  2:66210/149‑6 (MQ=255)
                                              aTGGTCTCCAGATGATCGATCAGTGCCGGGCGGAAGCCATCCAGCATTTCATGTACGGCAAAGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTGTGGTCATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAActc                                                                                                    <  6:111068/149‑1 (MQ=255)
                                              aTGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAActc                                                                                                    >  2:305533/1‑149 (MQ=255)
                                               tGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAActc                                                                                                    >  8:35299/1‑148 (MQ=255)
                                                    tCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGCAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGt                                                                                              >  2:156185/1‑149 (MQ=255)
                                                    tCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGCAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGt                                                                                              >  2:161660/1‑149 (MQ=255)
                                                     ccAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGtt                                                                                             >  2:242781/1‑149 (MQ=255)
                                                      cAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTg                                                                                            >  4:276113/1‑149 (MQ=255)
                                                       aGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGc                                                                                           >  1:311194/1‑149 (MQ=255)
                                                        gatgatCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCt                                                                                          <  4:254644/149‑1 (MQ=255)
                                                             tcgatcAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCtttttt                                                                                     >  5:343426/1‑149 (MQ=255)
                                                             tcgatcAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCtttttt                                                                                     >  1:224157/1‑149 (MQ=255)
                                                                atcAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTAAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTtctc                                                                                  >  8:150706/1‑149 (MQ=255)
                                                                    gtgcggtgcgGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGCAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTg                                                                              <  1:156185/149‑1 (MQ=255)
                                                                    gtgcggtgcgGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGCAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTg                                                                              <  1:161660/149‑1 (MQ=255)
                                                                     tgcggtgcggAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACt                                                                                                <  6:255576/130‑1 (MQ=255)
                                                                     tgcggtgcggAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACt                                                                                                >  5:255576/1‑130 (MQ=255)
                                                                       cggtgcggAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCg                                                                           <  7:144599/149‑1 (MQ=255)
                                                                         gtgcggAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGga                                                                         <  6:343426/149‑1 (MQ=255)
                                                                         gtgcggAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGga                                                                         <  1:403737/149‑1 (MQ=255)
                                                                         gtgcggAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGga                                                                         <  7:35299/149‑1 (MQ=255)
                                                                             ggAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACa                                                                     <  2:144551/149‑1 (MQ=255)
                                                                              gAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACAt                                                                    >  6:97196/1‑149 (MQ=255)
                                                                                aGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCg                                                                  <  2:230516/149‑1 (MQ=255)
                                                                                   cATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGAt                                                                                                                                      <  6:327334/78‑1 (MQ=255)
                                                                                   cATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGAt                                                                                                                                      >  5:327334/1‑78 (MQ=255)
                                                                                   cATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTg                                                               >  5:230093/1‑149 (MQ=255)
                                                                                    aTCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGc                                                              >  2:356982/1‑149 (MQ=255)
                                                                                    aTCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGc                                                              <  3:276113/149‑1 (MQ=255)
                                                                                     tCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCg                                                             <  1:356982/149‑1 (MQ=255)
                                                                                       cAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTTGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCgg                                                            <  8:1490/148‑1 (MQ=255)
                                                                                       cAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCgg                                                            <  5:97196/148‑1 (MQ=255)
                                                                                       cAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCgg                                                            <  1:220150/148‑1 (MQ=255)
                                                                                         gCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCGTGGTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGta                                                         >  7:300450/1‑149 (MQ=255)
                                                                                            tttCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGtata                                                       >  7:91768/1‑148 (MQ=255)
                                                                                             ttCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGtata                                                       >  4:35221/1‑147 (MQ=255)
                                                                                              tCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATGGTTGCGGGTATAAAg                                                    >  8:105686/1‑149 (MQ=255)
                                                                                                 tGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAg                                                    >  2:336817/1‑146 (MQ=255)
                                                                                                       gCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGt                                           <  5:171976/149‑1 (MQ=255)
                                                                                                         aaTGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCg                                         >  5:166453/1‑149 (MQ=255)
                                                                                                         aaTGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCg                                         >  7:373222/1‑149 (MQ=255)
                                                                                                           tGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCt                                       >  8:26270/1‑149 (MQ=255)
                                                                                                           tGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGGTGCGGGTATAAAGCAGATTgggagct                                       >  1:38492/1‑144 (MQ=255)
                                                                                                               gTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTTGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGTGTATAAAGGAGATTGGTCGCTTTTg                                   <  8:300450/149‑1 (MQ=255)
                                                                                                               gTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCtttt                                    >  2:308814/1‑148 (MQ=255)
                                                                                                               gTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTg                                   <  7:150706/149‑1 (MQ=255)
                                                                                                                   gcgcCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGCTGGGTCGCTTTTGAtt                                >  6:323962/1‑148 (MQ=255)
                                                                                                                       cgcgCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACt                           <  1:242781/149‑1 (MQ=255)
                                                                                                                              gTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAg                    >  8:94071/1‑149 (MQ=255)
                                                                                                                                 ccAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAGCg                  <  1:122939/148‑1 (MQ=255)
                                                                                                                                    gTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGTAGATTGGTCGCTTTTGATTTAACTAATTTAGCGGTAc              >  7:107286/1‑149 (MQ=255)
                                                                                                                                       cgcTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAGCGGTACTCa           <  8:393945/149‑1 (MQ=255)
                                                                                                                                        gcTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAGCGGTACTCa           >  6:170511/1‑148 (MQ=255)
                                                                                                                                         gttGTTTGGTAATCAAAGAAAGATCAATAAACTGGATTACCTGGCGATTCAGCAACTCAACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAGCGGTACTCAt          >  8:240667/2‑148 (MQ=255)
                                                                                                                                             ttttGTAATCAAAGAAAGATCAATAATCTGGATAACCTGGCGATTCAGAAACTCTACTGTTGCTTTTTTCTCGCTGTCTGAGAAATCGTTGCGGGTATAAAGCAGATTGGTCTCTTTTGATTTAACTAATTTAGCGGTACTCATAAttt     <  8:391168/149‑1 (MQ=255)
                                                                                                                                                 gTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAGCGGTACTCATAATTTCat  >  7:348424/1‑148 (MQ=255)
                                                                                                                                                 gTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAGCGGTACTCATAATTTCat  <  8:348424/148‑1 (MQ=255)
                                                                                                                                                 gTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAGCGGTACTCATAATTTCat  <  5:170511/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGTGGTCCCCAGAGCTACACCGCCCAGCTGCACTGCACGTTCTGCCATGGTATCCAGATGATCGATCAGTGCGGTGCGGAAGCCATCCAGCATTTCATGTACGGCAATGAAGTTAGCGCCGCGCATGTTCCAGTGCGCTTGTTTGGTAATCAAAGAAAGATCAATAAACTGGATAACCTGGCGATTCAGCAACTCTACTGTTGCTTTTTTCTCGCTGTCGGAGACATCGTTGCGGGTATAAAGCAGATTGGTCGCTTTTGATTTAACTAATTTAGCGGTACTCATAATTTCAT  >  NC_000913/848627‑848919

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: