Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,877,788 A→T 5.5% Q25L (CAG→CTG)  yeaP → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,877,7880AT 5.5% 229.2 / 2.7 73Q25L (CAG→CTG) yeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base T (2/2);  total (35/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 = 9.61e-01

TTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTAC  >  NC_000913/1877644‑1877936
                                                                                                                                                |                                                                                                                                                    
ttgATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCAATTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCtt                                                                                                                                                  >  2:187204/1‑149 (MQ=255)
 tgatgaTGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTTTAAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGGCACTGGAAAGTCTGGTCCGACAGCTTc                                                                                                                                                 >  2:350073/1‑149 (MQ=255)
  gatgatGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTc                                                                                                                                                 <  3:166030/148‑1 (MQ=255)
   atgatgTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCt                                                                                                                                                <  6:382027/148‑1 (MQ=255)
      atgTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGgag                                                                                                                                            <  4:325294/149‑1 (MQ=255)
      atgTTTTAGCATGCCTGCTATGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGga                                                                                                                                             <  6:85038/148‑1 (MQ=255)
        gTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGAt                                                                                                                                          >  1:51331/1‑149 (MQ=255)
         ttttATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATg                                                                                                                                         >  8:3689/1‑149 (MQ=255)
            tATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGGGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCt                                                                                                                                       <  2:345742/148‑1 (MQ=255)
             aTCATGCCTGCAAAGATTAAATAATCAGCGTTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTg                                                                                                                                      <  1:17370/148‑1 (MQ=255)
             aTCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTg                                                                                                                                      >  4:96881/1‑148 (MQ=255)
                 tGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAAt                                                                                                                                 <  1:276662/149‑1 (MQ=255)
                       cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGATTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCAAAAGAACAGTCACTGGAAAGTCTGGTCCGACTGCTTCTGGAGATGCTGGAAATGGTCAc                                                                                                                           >  2:248962/1‑149 (MQ=255)
                       cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACTGCTTCTGGAGATGCTGGAAATGGGCAc                                                                                                                           >  2:242798/1‑149 (MQ=255)
                       cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTc                                                                                                                             <  7:146106/147‑1 (MQ=255)
                       cAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCa                                                                                                                            <  3:43422/148‑1 (MQ=255)
                        aaaGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACt                                                                                                                          >  5:337053/1‑149 (MQ=255)
                         aaGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCAGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACt                                                                                                                          >  6:241470/1‑148 (MQ=255)
                                       aGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCCGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACATTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAacc                                                                                                           <  8:316954/149‑1 (MQ=255)
                                       aGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAac                                                                                                            <  8:240205/148‑1 (MQ=255)
                                          aTTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAacctac                                                                                                        >  8:417269/1‑149 (MQ=255)
                                           tttACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAacctccc                                                                                                       <  3:96881/149‑1 (MQ=255)
                                             tACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAacctacct                                                                                                      <  3:352370/148‑1 (MQ=255)
                                             tACCAGCCGTAACCTGGAGTTGTTCAGTGACAGATCAGATAATCGCCCTCGTCAAGCAAACCCTTGCCAAAGCACAGTCAAAGGAAAGTCAGGTCCGACAGCTTCTTGAGATGCTGGAAATGGTCACTGATATGCAATCAaccaacct                                                                                                      <  2:81225/148‑1 (MQ=255)
                                               cccGCAGTATCCTGGAGTTGTTCCGTTTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGGCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGa                                                                                                    <  1:231411/148‑1 (MQ=255)
                                                   ccgTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCaaa                                                                                               <  5:147438/149‑1 (MQ=255)
                                                      tATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTg                                                                                                                                      >  1:95204/1‑107 (MQ=255)
                                                      tATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTg                                                                                                                                      <  2:95204/107‑1 (MQ=255)
                                                        tCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTgg                                                                                           <  5:241470/148‑1 (MQ=255)
                                                        tCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGa                                                                                          <  4:389738/149‑1 (MQ=255)
                                                             gAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGt                                                                                       <  2:51331/147‑1 (MQ=255)
                                                                   ttCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgc                                                                               >  6:284132/1‑149 (MQ=255)
                                                                          tcagatcagaTTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTg                                                                          >  7:402284/1‑147 (MQ=255)
                                                                          tcagatcagaTTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTg                                                                          >  2:309185/1‑147 (MQ=255)
                                                                             gatcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGCCCACAGTGGATGTCGACGCGCGCCTgcagca                                                                     >  8:343048/1‑149 (MQ=255)
                                                                              atcagatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgcagca                                                                     >  4:131038/1‑148 (MQ=255)
                                                                                 agatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTg                                                                                                                                               <  2:310105/71‑1 (MQ=255)
                                                                                 agatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTg                                                                                                                                               >  1:310105/1‑71 (MQ=255)
                                                                                 agatTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCAtatt                                                                 <  5:73775/149‑1 (MQ=255)
                                                                                       tCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGtt                                                            <  1:350073/148‑1 (MQ=255)
                                                                                         gCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGc                                                         >  8:123664/1‑149 (MQ=255)
                                                                                         gCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGc                                                         <  1:309185/149‑1 (MQ=255)
                                                                                           ccGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGccc                                                       >  1:298637/1‑149 (MQ=255)
                                                                                           ccGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGccc                                                       <  1:187204/149‑1 (MQ=255)
                                                                                                tctcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTaa                                                   >  5:392518/1‑148 (MQ=255)
                                                                                                  tcGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAg                                                >  8:282954/1‑149 (MQ=255)
                                                                                                   cGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGc                                               >  7:112935/1‑149 (MQ=255)
                                                                                                    gCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGc                                               >  7:112665/1‑148 (MQ=255)
                                                                                                      aaTCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAg                                            <  3:131038/149‑1 (MQ=255)
                                                                                                          ccTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGaaaa                                        >  1:380840/1‑149 (MQ=255)
                                                                                                              gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATg                                      <  6:392518/147‑1 (MQ=255)
                                                                                                              gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGt                                     <  7:3689/148‑1 (MQ=255)
                                                                                                              gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTa                                    >  2:387815/1‑149 (MQ=255)
                                                                                                              gCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCACGTAACAGCCAGAAAATGt                                     <  6:337053/148‑1 (MQ=255)
                                                                                                                cAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACa                                  >  8:201436/1‑149 (MQ=255)
                                                                                                                       cAGTCACTGGAAAGTCTGGTCCGACTGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGga                           <  1:242798/149‑1 (MQ=255)
                                                                                                                       cAGTCACTGGAAAGTCTGGTCCGACTGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGga                           <  1:248962/149‑1 (MQ=255)
                                                                                                                          tCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAATGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGaa                        >  4:149915/1‑149 (MQ=255)
                                                                                                                          tCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAG‑AAATGTACATCCCGGAGAAt                       >  1:421041/1‑149 (MQ=255)
                                                                                                                            aCTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtt                      >  8:64293/1‑149 (MQ=255)
                                                                                                                            aCTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAtt                      <  8:106612/149‑1 (MQ=255)
                                                                                                                              tGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAAttt                     >  6:248670/1‑148 (MQ=255)
                                                                                                                                gAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACTTCCCGGAGAATTTTAc                  >  5:152564/1‑149 (MQ=255)
                                                                                                                                gAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTa                   <  5:284132/148‑1 (MQ=255)
                                                                                                                                gAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTAc                  <  8:112665/149‑1 (MQ=255)
                                                                                                                                  aaGTCTGGTCCGACAGCTTCTGGATATGCTGTATATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGCATGTAGAAGCGATACTGCAGCATATTATTTTTGCCCGTAACAGCAATAAAATCTAAATCACGGAGAATTTTAccc                >  1:152047/1‑148 (MQ=255)
                                                                                                                                         gTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGtctc            <  2:298637/146‑1 (MQ=255)
                                                                                                                                         gTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGtctc            <  1:387815/146‑1 (MQ=255)
                                                                                                                                         gTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTg         <  8:112935/149‑1 (MQ=255)
                                                                                                                                            cGACAGCTTCTCGAGATGCAGGAAATGGTCACAGATATCGAATCAACCTACCAGACCAAAGTGGAAGTCCAAGCGCGCCTGCACCATATTATGTTTGCCCGTAACAGCCAGAAAATGTAGATCCCGGAGAATTTTACCGTCTGGTGCGa      >  4:165768/1‑149 (MQ=255)
                                                                                                                                             gACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGAt     >  4:936/1‑149 (MQ=255)
                                                                                                                                             gACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGAt     >  8:210510/1‑149 (MQ=255)
                                                                                                                                               cAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTa   <  7:193384/149‑1 (MQ=255)
                                                                                                                                                aGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTAc  <  7:201436/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCCCGGAGAATTTTACCGTCTCGTGGGATTAC  >  NC_000913/1877644‑1877936

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