Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 907,330 C→T 5.8% R318Q (CGA→CAA)  ybjT ← putative NAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913907,3300CT 5.8% 223.3 / 2.4 69R318Q (CGA→CAA) ybjTputative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base T (2/2);  total (36/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.40e-01

CAACGCCCGTGTCTGCCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAG  >  NC_000913/907185‑907475
                                                                                                                                                 |                                                                                                                                                 
cAACGCCCGTGTCTGCCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGc                                                                                                                                                 >  7:1878/1‑148 (MQ=255)
cAACGCCCGTGTCTGCCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGc                                                                                                                                                 >  1:397376/1‑148 (MQ=255)
  aCGCCCGTGTCTGTCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGc                                                                                                                                                 >  8:404412/1‑146 (MQ=255)
  aCGCCCGTGTCTGCCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGCATCACTACCTGCCATAAAGAAACAAGGCTGCCGGACGTTTTAACGGTAAACCCCGCCTGTTTCGCAAAATAACCGTACTCCGGTCGcaa                                                                                                                                               >  6:326096/1‑146 (MQ=255)
     cccGTGTCTGCCACAAAATATTTCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCg                                                                                                                                            >  2:340790/1‑148 (MQ=255)
     cccGTGTCTGCCACAAAATATTGCCAAAGAAATAACGCTATTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGCCGGACGTTTTAACGGTAAACCACGCCTGTTTGGCAAAATAACCGTACTCCGTTCGCCAGCg                                                                                                                                            >  2:207714/1‑148 (MQ=255)
          gtCTGCCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCa                                                                                                                                        >  1:257479/1‑147 (MQ=255)
           tCTGTCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCaa                                                                                                                                       <  7:404412/147‑1 (MQ=255)
           tCTGCCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCaa                                                                                                                                       <  8:16521/147‑1 (MQ=255)
           tCTGCCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCaa                                                                                                                                       <  2:397376/147‑1 (MQ=255)
              gCCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTACCGGTAAACCCCGCCTGGTTGGCAAAATACCCGTACTCCGGTCGCCAGCGGGCAAAGGc                                                                                                                                   >  6:116665/1‑148 (MQ=255)
               ccACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAgg                                                                                                                                    <  3:353397/146‑1 (MQ=255)
                cacaAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCt                                                                                                                                 >  7:273428/1‑148 (MQ=255)
                 acaAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTg                                                                                                                                <  7:296186/148‑1 (MQ=255)
                  caAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAgg                                                                                                                                    <  2:17420/143‑1 (MQ=255)
                  caAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAgg                                                                                                                                    >  1:17420/1‑143 (MQ=255)
                   aaaaTATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAg                                                                                                                              >  1:96620/1‑148 (MQ=255)
                    aaaTATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGcc                                                                                                                                  >  8:110453/1‑143 (MQ=255)
                    aaaTATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGcc                                                                                                                                  <  7:110453/143‑1 (MQ=255)
                       tatTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGcgtc                                                                                                                          <  2:96620/148‑1 (MQ=255)
                         ttGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGcgtcg                                                                                                                         <  4:36824/147‑1 (MQ=255)
                             cAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGcc                                                                                                                    >  5:387812/1‑148 (MQ=255)
                               aaGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCaa                                                                                                                                       <  5:311488/127‑1 (MQ=255)
                               aaGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCaa                                                                                                                                       >  6:311488/1‑127 (MQ=255)
                               aaGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGcccc                                                                                                                  <  6:263634/148‑1 (MQ=255)
                               aaGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGcc                                                                                                                    <  7:78057/146‑1 (MQ=255)
                               aaGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGcc                                                                                                                    >  8:78057/1‑146 (MQ=255)
                               aaGAAATAACGCTCTCTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTTGCCAGCGGGCAAAGGCCTGAGCGTCGTAGcccc                                                                                                                  >  2:150508/1‑148 (MQ=255)
                               aaGAAATAACGCTCTCTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTTGCCAGCGGGCAAAGGCCTGAGCGTCGTAGccc                                                                                                                   >  2:156159/1‑147 (MQ=255)
                                      aaCGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCt                                                                                                           >  4:186560/1‑148 (MQ=255)
                                       aCGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCt                                                                                                           <  8:240287/147‑1 (MQ=255)
                                        cGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTCGCGGACGTTTTAACGGTAACCCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTg                                                                                                          >  6:373743/1‑147 (MQ=255)
                                             tttACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGt                                                                                                      <  8:315712/146‑1 (MQ=255)
                                              ttACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGa                                                                                                   >  4:348060/1‑148 (MQ=255)
                                               tACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGGCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTg                                                                                                    <  2:175852/146‑1 (MQ=255)
                                                aCCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGAc                                                                                                  >  6:174426/1‑147 (MQ=255)
                                                   gccgATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGt                                                                                              >  1:23210/1‑148 (MQ=255)
                                                         ttGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTtc                                                                                         <  2:168014/147‑1 (MQ=255)
                                                             ttCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTcctcct                                                                                    >  4:77102/1‑148 (MQ=255)
                                                                   aCCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTc                                                                               >  6:38851/1‑147 (MQ=255)
                                                                   aCCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTc                                                                               >  6:195510/1‑147 (MQ=255)
                                                                    ccTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTc                                                                               <  8:1878/146‑1 (MQ=255)
                                                                    ccTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCaa                                                                             <  4:74531/148‑1 (MQ=255)
                                                                         cATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGttc                                                                        <  6:387812/148‑3 (MQ=255)
                                                                                cagcaAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTAc                                                                 >  7:80620/1‑148 (MQ=255)
                                                                                  gcaAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGttttt                                                                                          <  2:212356/121‑1 (MQ=255)
                                                                                  gcaAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGttttt                                                                                          >  1:212356/1‑121 (MQ=255)
                                                                                      ggCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTTGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTc                                                           <  1:150508/148‑1 (MQ=255)
                                                                                      ggCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTTGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTc                                                           <  1:156159/148‑1 (MQ=255)
                                                                                      ggCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTc                                                           <  3:77102/148‑1 (MQ=255)
                                                                                              gACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAg                                                   >  5:59035/1‑148 (MQ=255)
                                                                                                 gTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGa                                                <  8:80620/148‑1 (MQ=255)
                                                                                                 gTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGa                                                <  1:340790/148‑1 (MQ=255)
                                                                                                   tttAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGAt                                               >  3:156773/1‑147 (MQ=255)
                                                                                                     tAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCa                                             <  2:257479/147‑1 (MQ=255)
                                                                                                      aaCGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGc                                           >  7:224452/1‑148 (MQ=255)
                                                                                                        cGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCg                                         >  7:88709/1‑148 (MQ=255)
                                                                                                           tAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGTTCGCCAGCGGGCAAAGGCCTGAGCGTCGTATACCCACACGCTGTAGTTTACCAGTTATTACTACTCTTTCAACGTTCAACTTACAGATTAATCGAAAGCGATCATCAGTtg                                      >  4:323709/1‑148 (MQ=255)
                                                                                                               ccccGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTggg                                  <  5:38851/148‑1 (MQ=255)
                                                                                                                  cGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGAt                                >  8:402688/1‑147 (MQ=255)
                                                                                                                           ggCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACg                       >  8:144994/1‑147 (MQ=255)
                                                                                                                                aaTAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAgcgc                 >  1:143367/1‑148 (MQ=255)
                                                                                                                                  tAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCgg               <  2:23210/148‑1 (MQ=255)
                                                                                                                                    aCCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCGGt              <  8:224452/147‑1 (MQ=255)
                                                                                                                                    aCCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCGGt              <  2:143367/147‑1 (MQ=255)
                                                                                                                                     ccGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCGGTat            <  8:273428/148‑1 (MQ=255)
                                                                                                                                            ccGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGc     >  2:74692/1‑148 (MQ=255)
                                                                                                                                             cGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCg      <  5:116665/146‑1 (MQ=255)
                                                                                                                                               ggcGCCAGCGGGCAAAGGCCGGAGCGTCGTAGCCCCAGTCGCTGGAGGTGACCAGTTTTGCCTCCTCTTTCCACGTGCTACGTACCGCGTCAGCGAAAGCGAGCAGCCGTTGGGGGATGAGAGCACGTAGCGCGGTATCATCCGCcag  <  5:9679/146‑1 (MQ=255)
                                                                                                                                               gTCGCCAGCGGGCAAAGGCCTGAGCGTAGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCcag  <  5:326096/148‑1 (MQ=255)
                                                                                                                                                 cGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCcag  >  3:226590/1‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CAACGCCCGTGTCTGCCACAAAATATTGCCAAAGAAATAACGCTCTTTACCGCCGATTTGGTTCACTACCTGCCATAAAGCAGCAAGGCTGGCGGACGTTTTAACGGTAAACCCCGCCTGTTTGGCAAAATAACCGTACTCCGGTCGCCAGCGGGCAAAGGCCTGAGCGTCGTAGCCCCAGTCGCTGGAGTTGACCAGTTTTTCCTCCTCTTTCAACGTGCTACGTACCGCGTCATCGAAAGCGATCAGCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAG  >  NC_000913/907185‑907475

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: ‑