Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,282,866 T→C 6.0% A297A (GCT→GCC agaS → tagatose‑6‑phosphate ketose/aldose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,282,8660TC 6.0% 205.8 / 3.1 67A297A (GCT→GCCagaStagatose‑6‑phosphate ketose/aldose isomerase
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base C (2/2);  total (33/34)
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.70e-01

TGACGGCGGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCC  >  NC_000913/3282721‑3283013
                                                                                                                                                 |                                                                                                                                                   
tGACGGCGGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGTATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGac                                                                                                                                                  >  6:249888/1‑148 (MQ=255)
tGACGGCGGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGaa                                                                                                                                                  >  1:362431/1‑149 (MQ=255)
  aCGGCGGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAAc                                                                                                                                                 <  6:88416/148‑1 (MQ=255)
    ggcggGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCACTGGGCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTc                                                                                                                                              <  2:136791/149‑1 (MQ=255)
       gggTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGcc                                                                                                                                           <  2:362431/149‑1 (MQ=255)
        ggTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATTGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCg                                                                                                                                          >  8:125867/1‑149 (MQ=255)
         gTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGt                                                                                                                                         <  6:221150/149‑1 (MQ=255)
          tAAACTGGCGGCCTTTTATGATTCTCCATCCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTg                                                                                                                                        <  3:68005/149‑1 (MQ=255)
            aaCTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGa                                                                                                                                       <  6:208372/148‑1 (MQ=255)
            aaCTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGAc                                                                                                                                      <  5:88528/149‑1 (MQ=255)
                  cggcCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAg                                                                                                                                >  2:103695/1‑149 (MQ=255)
                  cggcCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAg                                                                                                                                <  1:184394/149‑1 (MQ=255)
                   ggcCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAgg                                                                                                                               >  2:245634/1‑149 (MQ=255)
                      cTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGGCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCaa                                                                                                                            <  3:96685/149‑1 (MQ=255)
                       ttttATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAAt                                                                                                                           >  2:181512/1‑149 (MQ=255)
                          tATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGc                                                                                                                         <  1:164305/148‑1 (MQ=255)
                           aTGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCgt                                                                                                                       >  1:331509/1‑149 (MQ=255)
                           aTGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCgt                                                                                                                       >  6:166018/1‑149 (MQ=255)
                            tGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCgtg                                                                                                                      <  8:68661/149‑1 (MQ=255)
                                 ctcCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAat                                                                                                                  <  1:181512/148‑1 (MQ=255)
                                    cAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAatc                                                                                                                 <  7:302411/146‑1 (MQ=255)
                                     aaCCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGt                                                                                                                                         <  6:54253/121‑1 (MQ=255)
                                        cGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCCGAACTTCGCCGTGACAACCAGGCAATGCGTGTAatcgccatcg                                                                                                          >  6:322025/1‑149 (MQ=255)
                                                 tCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGa                                                                                                   >  5:262570/1‑147 (MQ=255)
                                                    tGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAagca                                                                                              <  7:10586/149‑1 (MQ=255)
                                                      gACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAagcagc                                                                                            <  4:172696/149‑1 (MQ=255)
                                                       aCCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGACcgcagc                                                                                            >  6:194129/1‑148 (MQ=255)
                                                         cAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGAc                                                                                         <  6:262570/149‑1 (MQ=255)
                                                          aaaaTCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACa                                                                                        >  6:129156/1‑149 (MQ=255)
                                                             aTCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCATGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtcg                                                                                     <  1:103695/149‑1 (MQ=255)
                                                                    gTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCgt                                                                                                                       >  5:72347/1‑108 (MQ=255)
                                                                    gTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCgt                                                                                                                       <  6:72347/108‑1 (MQ=255)
                                                                    gTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGc                                                                              >  3:304252/1‑149 (MQ=255)
                                                                          gACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCGGCCGGTcc                                                                        >  2:87215/1‑149 (MQ=255)
                                                                                aCGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACatat                                                                  >  5:166298/1‑149 (MQ=255)
                                                                                  gCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCCGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatca                                                                <  5:322025/149‑1 (MQ=255)
                                                                                     ggtggtggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATcc                                                             >  7:263340/1‑149 (MQ=255)
                                                                                           ggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACat                                                                    <  4:136251/136‑1 (MQ=255)
                                                                                           ggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACat                                                                    >  3:136251/1‑136 (MQ=255)
                                                                                            gtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCAcc                                                      >  1:221675/1‑149 (MQ=255)
                                                                                               tttGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGcctcagt                                                    >  2:254880/1‑143 (MQ=255)
                                                                                               tttGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtc                                                   <  2:331509/149‑1 (MQ=255)
                                                                                                ttGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGTGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtca                                                  >  7:160608/1‑149 (MQ=255)
                                                                                                   tctcCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgt                                               <  5:129156/149‑1 (MQ=255)
                                                                                                   tctcCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgt                                               <  4:304252/149‑1 (MQ=255)
                                                                                                      ccAGCCACCCTTACACCCGTCAGTATTATCTTGATCTGCTGGCTGAACTTCGCCGTCACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCgtcacgtcac                                            <  5:194129/149‑1 (MQ=255)
                                                                                                       cAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGTGACATCGTCGCTGCCGGTCCAAATATCATCCTGCCACCGTCACGTCACt                                           <  7:125867/149‑1 (MQ=255)
                                                                                                       cAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCgtcacgtcac                                            <  5:166018/148‑1 (MQ=255)
                                                                                                          ccACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTc                                                                         <  4:31726/116‑1 (MQ=255)
                                                                                                          ccACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTa                                        >  4:176223/1‑149 (MQ=255)
                                                                                                           cACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTAt                                       >  7:153666/1‑149 (MQ=255)
                                                                                                            aCCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCCGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATc                                      >  2:152596/1‑149 (MQ=255)
                                                                                                             cccTTAAACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATc                                      >  4:74417/1‑148 (MQ=255)
                                                                                                              ccTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGCCAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGa                                    >  1:167328/1‑149 (MQ=255)
                                                                                                                ttACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACg                                  >  5:178644/1‑149 (MQ=255)
                                                                                                                ttACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACg                                  >  7:205879/1‑149 (MQ=255)
                                                                                                                   cacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGtt                                <  8:263340/148‑1 (MQ=255)
                                                                                                                   cacCAGTCAGTATCATCTTGAGCTGCTTGGTGAACTTAGCCGTGACAACCAGGCCATGCGTGTAATCGCCATCGCAGCGGAAAGCAGCGACATCGTCGCTGACGGTCCACATATCCTCCTGCCACCGTCACGTCACTTTTTCGACGtt                                <  6:255491/148‑1 (MQ=255)
                                                                                                                      ccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGc                            >  4:71840/1‑149 (MQ=255)
                                                                                                                            gTATGATCTTGATCTGCTGGCCGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAt                      <  1:152596/149‑1 (MQ=255)
                                                                                                                             tATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCgg                                                                           >  1:151188/1‑95 (MQ=255)
                                                                                                                             tATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCgg                                                                           <  2:151188/95‑1 (MQ=255)
                                                                                                                             tATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAtt                     >  6:81250/1‑149 (MQ=255)
                                                                                                                                 atcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCTTCACGTCACTTTATCGACGTGGAGCAGGCATTTTGc                 >  6:86136/1‑149 (MQ=255)
                                                                                                                                 atcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCAttt                    <  1:245634/146‑1 (MQ=255)
                                                                                                                                    ttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGACACCGTCACGTCACTTTATCGACGTTGATCAGGCATTTTGCTTc              >  6:251077/1‑149 (MQ=255)
                                                                                                                                      gatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTc              <  6:178644/147‑1 (MQ=255)
                                                                                                                                            ctggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGAt         <  5:81250/146‑1 (MQ=255)
                                                                                                                                                 tgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGcc  >  3:103209/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGACGGCGGGTAAACTGGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCC  >  NC_000913/3282721‑3283013

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