Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,677,301 T→C 6.1% K201K (AAA→AAG pntA ← pyridine nucleotide transhydrogenase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,677,3010TC 6.1% 211.5 / 3.8 66K201K (AAA→AAGpntApyridine nucleotide transhydrogenase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (35/27);  new base C (2/2);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

ATCGACCTCTTTTGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAG  >  NC_000913/1677154‑1677448
                                                                                                                                                   |                                                                                                                                                   
atCGACCTCTTTTGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCtt                                                                                                                                                    <  4:514204/149‑1 (MQ=255)
atCGACCTCTTTTGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCtt                                                                                                                                                    >  1:605816/1‑149 (MQ=255)
 tCGACCTCTTTTGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCttt                                                                                                                                                   <  2:89916/149‑1 (MQ=255)
   gACCTCTTTTGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCa                                                                                                                                                 >  3:684662/1‑149 (MQ=255)
    aCCTCTTTTGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCAc                                                                                                                                                >  3:459309/1‑149 (MQ=255)
    aCCTCTTTTGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCAc                                                                                                                                                >  4:436972/1‑149 (MQ=255)
         ttttGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCt                                                                                                                                                     >  1:382029/1‑139 (MQ=255)
         ttttGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCt                                                                                                                                                     <  2:382029/139‑1 (MQ=255)
              ccTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCcgggcgg                                                                                                                                      <  2:605816/149‑1 (MQ=255)
                   gcggcAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTc                                                                                                                                 <  4:459309/149‑1 (MQ=255)
                   gcggcAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTc                                                                                                                                 <  1:361447/149‑1 (MQ=255)
                      gcAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGaa                                                                                                                              <  2:172751/149‑1 (MQ=255)
                       cAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAAt                                                                                                                             <  4:730539/149‑1 (MQ=255)
                           gagTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATgcac                                                                                                                         <  3:300874/149‑1 (MQ=255)
                              ttCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATgcacgca                                                                                                                      >  4:185575/1‑149 (MQ=255)
                              ttCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATgcacgca                                                                                                                      >  4:581864/1‑149 (MQ=255)
                                  aTTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAAt                                                                                                                  >  4:489397/1‑149 (MQ=255)
                                  aTTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAAt                                                                                                                  >  4:151644/1‑149 (MQ=255)
                                     tCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATcgc                                                                                                               <  3:489397/149‑1 (MQ=255)
                                       cGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTCGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCgcgc                                                                                                             <  2:695950/149‑1 (MQ=255)
                                         cTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCg                                                                                                           >  4:233857/1‑149 (MQ=255)
                                            tgatgaACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGcc                                                                                                            >  4:630505/1‑145 (MQ=255)
                                            tgatgaACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGcc                                                                                                            <  3:630505/145‑1 (MQ=255)
                                             gatgaACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGAc                                                                                                       >  4:659826/1‑149 (MQ=255)
                                                  aCGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGttt                                                                                                  >  1:458421/1‑149 (MQ=255)
                                                   cgcgTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTtg                                                                                                 <  3:362884/149‑1 (MQ=255)
                                                      gTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTtgctg                                                                                              <  3:659826/149‑1 (MQ=255)
                                                              tCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATg                                                                                      >  4:97136/1‑149 (MQ=255)
                                                              tCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATg                                                                                      >  2:47189/1‑149 (MQ=255)
                                                              tCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATg                                                                                      >  1:815337/1‑149 (MQ=255)
                                                                        aTAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGAc                                                                            >  2:147853/1‑149 (MQ=255)
                                                                        aTAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGAc                                                                            <  2:289723/149‑1 (MQ=255)
                                                                            ccATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGc                                                                        >  2:341548/1‑149 (MQ=255)
                                                                             cATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCa                                                                       >  3:498005/1‑149 (MQ=255)
                                                                              aTCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATgcac                                                                                                                         >  3:465490/1‑98 (MQ=255)
                                                                              aTCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCaa                                                                      >  3:541193/1‑149 (MQ=255)
                                                                              aTCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTCTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATgcac                                                                                                                         <  4:465490/98‑1 (MQ=255)
                                                                               tCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAac                                                                     <  2:350585/149‑1 (MQ=255)
                                                                                          aGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCaa                                                          <  2:270895/149‑1 (MQ=255)
                                                                                          aGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCaa                                                          >  3:726260/1‑149 (MQ=255)
                                                                                                 tctTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCAt                                                   >  4:119001/1‑149 (MQ=255)
                                                                                                 tctTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCAt                                                   <  3:97136/149‑1 (MQ=255)
                                                                                                       aaaTCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACttt                                             <  4:498005/149‑1 (MQ=255)
                                                                                                        aaTCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGcc                                                                                >  3:158156/1‑113 (MQ=255)
                                                                                                        aaTCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGcc                                                                                <  4:158156/113‑1 (MQ=255)
                                                                                                          tCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGc                                          <  4:289399/149‑1 (MQ=255)
                                                                                                            cAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCg                                        <  3:581864/149‑1 (MQ=255)
                                                                                                             aGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCgg                                       >  4:220079/1‑149 (MQ=255)
                                                                                                              gCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGt                                      <  2:458421/149‑1 (MQ=255)
                                                                                                                tCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCCTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTgg                                    >  4:147686/1‑149 (MQ=255)
                                                                                                                tCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCCTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTgg                                    >  4:162515/1‑149 (MQ=255)
                                                                                                                  gAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCa                                  >  2:511081/1‑149 (MQ=255)
                                                                                                                    ggAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGAc                                                                            <  1:309843/105‑1 (MQ=255)
                                                                                                                    ggAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGAc                                                                            >  2:309843/1‑105 (MQ=255)
                                                                                                                    ggAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACt                                >  3:102677/1‑149 (MQ=255)
                                                                                                                     gAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAAt                                                                                       >  4:2202/1‑93 (MQ=255)
                                                                                                                     gAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAAt                                                                                       <  3:2202/93‑1 (MQ=255)
                                                                                                                         tCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTccc                           <  2:159385/149‑1 (MQ=255)
                                                                                                                                 cATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAg                   >  4:433744/1‑149 (MQ=255)
                                                                                                                                   tACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTa                 <  1:47189/149‑1 (MQ=255)
                                                                                                                                      tttGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAAtt              >  3:357236/1‑149 (MQ=255)
                                                                                                                                       ttGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAAttt             <  1:341548/149‑1 (MQ=255)
                                                                                                                                       ttGAACTTGTTCCTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAAttt             <  3:162515/149‑1 (MQ=255)
                                                                                                                                       ttGAACTTGTTCCTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAAttt             <  3:147686/149‑1 (MQ=255)
                                                                                                                                           aCTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGccc         >  3:393540/1‑149 (MQ=255)
                                                                                                                                            cTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCg        >  1:367845/1‑149 (MQ=255)
                                                                                                                                             ttGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCgg       >  2:208176/1‑149 (MQ=255)
                                                                                                                                             ttGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCgg       >  4:791962/1‑149 (MQ=255)
                                                                                                                                               gTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGcc          >  4:311571/1‑144 (MQ=255)
                                                                                                                                               gTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGcc          <  3:311571/144‑1 (MQ=255)
                                                                                                                                                  cTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAaag  >  1:639665/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCGACCTCTTTTGCCTGGGCGGCAAAGAGTTCCATTTCCGCTTTGATGAACGCGTCCGACATCACTTTGGCATAGCCATCGCCGCTGCCAGCTTCCTCTTTAAAATCCAGCTCGAGGAATTCCGCGCCCATACTTTGAACTTGTTCTTTCACTTCCGGGCGGGTGTCGAATGCACGCACAATCGCGCCGAGACTGTTTGCTGCGCCAATGGCGGCCAGACCTGCAACACCCGCACCAATCACCATCACTTTTGCCGGTGGCACTTTCCCGGCCGCAGTAATTTGCCCGGTAAAG  >  NC_000913/1677154‑1677448

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