Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 347,692 C→A 9.5% G251V (GGC→GTC)  prpR ← propionate catabolism operon regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913347,6920CA 9.5% 219.6 / 5.2 73G251V (GGC→GTC) prpRpropionate catabolism operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base A (2/6);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCATAGCCAAACAGTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGGCATCG  >  NC_000913/347546‑347838
                                                                                                                                                  |                                                                                                                                                  
ccTCATAGCCAAACAGTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGcc                                                                                                                                                   >  7:111914/1‑148 (MQ=255)
ccTCATAGCCAAACAGTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGcc                                                                                                                                                   >  2:332455/1‑148 (MQ=255)
         cAAACAGTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGt                                                                                                                                          <  3:233764/148‑1 (MQ=255)
         cAAACAGTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGt                                                                                                                                          <  6:301883/148‑1 (MQ=255)
          aaaCAGTTCTGATTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTTTCGCGCGTCGTGGCGGGCAAAATATTACCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGtt                                                                                                                                         >  1:86395/1‑148 (MQ=255)
           aaCAGTTCTGCTTCCCGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGttt                                                                                                                                        >  3:281011/1‑148 (MQ=255)
           aaCAGTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGttt                                                                                                                                        <  3:263681/148‑1 (MQ=255)
           aaCAGTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCATGGGCCGCAAGCTCTTTGCCCGTCCCCGgtt                                                                                                                                        >  6:295058/1‑148 (MQ=255)
               gTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGc                                                                                                                                                    <  3:170901/132‑1 (MQ=255)
               gTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGc                                                                                                                                                    >  4:170901/1‑132 (MQ=255)
                 tCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCt                                                                                                                                  >  1:433324/1‑148 (MQ=255)
                  cTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCt                                                                                                                                  <  1:350227/147‑1 (MQ=255)
                  cTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCACt                                                                                                                                  >  2:45312/1‑147 (MQ=255)
                     cTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCtcaa                                                                                                                               <  7:282898/147‑1 (MQ=255)
                     cTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCtca                                                                                                                                <  7:208570/146‑1 (MQ=255)
                       tCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCtcaatc                                                                                                                             >  4:137637/1‑147 (MQ=255)
                       tCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCGTTGCCCGTCCCCGTTTCCCCCtcaatc                                                                                                                             >  3:4688/1‑147 (MQ=255)
                        ccAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGACAGCTCTTTGCCCGTCCCCGTTTCCCCCtcaatca                                                                                                                            <  6:184515/147‑1 (MQ=255)
                            cagcGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTGCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAaca                                                                                                                         >  5:150364/1‑146 (MQ=255)
                              gcGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACAccgc                                                                                                                     >  5:78291/1‑148 (MQ=255)
                                      gCAATCGCCCCGCAGTTGACGGCAACACACGGAGGCGACTTTTTGCCCTGT‑‑CGCATCGGGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCTATTTTTTTTAATGTTACGGTGTCCACCTATATCTACACCGCCGCACTGGAg                                                                                                           <  8:158205/148‑1 (MQ=255)
                                         aTCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGc                                                                                                          >  6:41036/1‑148 (MQ=255)
                                          tCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCACCTCAATCAACACCGCCGCACTGGAGCg                                                                                                         <  6:150364/148‑1 (MQ=255)
                                            gCCCCGCAGTTGACGTCAACACACGGCTGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCATTTTTTTTTAATGTTACGGTTTCCACCTATATCTACACCGCCGCACTGGAGCggg                                                                                                       <  2:115545/148‑1 (MQ=255)
                                             ccccGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGc                                                                                                      <  1:332455/148‑1 (MQ=255)
                                               ccGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCAt                                                                                                    <  6:146560/148‑1 (MQ=255)
                                                 gCAGTTGACGGCAACACACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATAc                                                                                                  <  5:245003/148‑1 (MQ=255)
                                                 gCAGTTGACGGCAACAACCGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGGCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATAc                                                                                                  <  3:331224/148‑1 (MQ=255)
                                                       gACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCATTTTTTTGAATGTTACCGTTTCCACCTATATCAACACCGCCGCACTGGAGCGGGCATACAGCaa                                                                                             <  8:179042/147‑1 (MQ=255)
                                                        aCGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCaaa                                                                                            <  4:4688/147‑1 (MQ=255)
                                                          ggCAACAAACGGAGGCGACTTTTGGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCATTTTTTTTTTTTAATGTTACGGTGACCACCGATATCTACACCGCCGCACTGGAGCGGGCATACAGCAAAATa                                                                                         <  2:254305/148‑1 (MQ=255)
                                                            caacaaACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGt                                                                                       >  4:153797/1‑148 (MQ=255)
                                                                 aaCGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCAGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCaa                                                                                             <  1:157224/137‑1 (MQ=255)
                                                                 aaCGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCAGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCaa                                                                                             >  2:157224/1‑137 (MQ=255)
                                                                   cGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGt                                                                                <  3:137637/148‑1 (MQ=255)
                                                                      aTGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACt                                                                             <  5:286415/148‑1 (MQ=255)
                                                                         cGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTg                                                                           <  5:41036/147‑1 (MQ=255)
                                                                         cGACTTTTTGCCCAGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGACCGACCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTATCCCGTACTTga                                                                          >  4:33814/1‑147 (MQ=255)
                                                                          gACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCCCTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGtt                                                                         >  6:339716/1‑148 (MQ=255)
                                                                           aCTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTc                                                                        <  7:280581/148‑1 (MQ=255)
                                                                           aCTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTc                                                                        >  8:280581/1‑148 (MQ=255)
                                                                            cTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTc                                                                        >  1:193358/1‑147 (MQ=255)
                                                                             tttttGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCACGATGACTCGCCTGGGCCGCCAGCTCTTTGCCCGTACCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTCCTGGTTCCa                                                                      >  1:139691/1‑148 (MQ=255)
                                                                              ttttGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACATCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCAt                                                                     >  4:260176/1‑148 (MQ=255)
                                                                                ttGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGTTCTTTGACCGTACCGGTTTCCCCCTCTATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATc                                                                    <  8:111914/147‑1 (MQ=255)
                                                                                   cccTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCtgtg                                                                >  5:256990/1‑148 (MQ=255)
                                                                                        tCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTgtggtg                                                             <  5:172526/146‑1 (MQ=255)
                                                                                        tCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGAt                                                           >  4:115840/1‑148 (MQ=255)
                                                                                        tCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGAt                                                           >  3:268354/1‑148 (MQ=255)
                                                                                        tCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTATGACCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGAt                                                           >  4:76043/1‑148 (MQ=255)
                                                                                         cgcgCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAAt                                                                                          >  1:116712/1‑116 (MQ=255)
                                                                                         cgcgCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAAt                                                                                          <  2:116712/116‑1 (MQ=255)
                                                                                           cgcATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGa                                                        <  2:193358/148‑1 (MQ=255)
                                                                                            gcATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACATCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGAc                                                       <  3:260176/148‑1 (MQ=255)
                                                                                                      gggCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTATTTTACCGTTCCCGTTTCCACCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATc                                             <  8:220897/148‑1 (MQ=255)
                                                                                                        gCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGc                                           >  8:56563/1‑148 (MQ=255)
                                                                                                        gCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGc                                           >  3:20046/1‑148 (MQ=255)
                                                                                                          aaaaTATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGccc                                         <  2:433324/148‑1 (MQ=255)
                                                                                                          aaaaTATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGccc                                         <  4:20046/148‑1 (MQ=255)
                                                                                                            aaTATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCtg                                                                  <  4:267474/121‑1 (MQ=255)
                                                                                                            aaTATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCtg                                                                  >  3:267474/1‑121 (MQ=255)
                                                                                                            aaTATTACCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCACGTTTCCCCCTCAATCAACACCGCCGCACTGGAGAGGGCATAAAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAg                                       >  1:155777/1‑148 (MQ=255)
                                                                                                             atatTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCAt                                                                     >  3:326770/1‑117 (MQ=255)
                                                                                                             atatTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCAt                                                                     <  4:326770/117‑1 (MQ=255)
                                                                                                             atatTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGc                                      >  6:450462/1‑148 (MQ=255)
                                                                                                                 tCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACg                                   >  2:209033/1‑147 (MQ=255)
                                                                                                                    cGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTaa                                <  8:380790/147‑1 (MQ=255)
                                                                                                                    cGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAAc                               >  6:43374/1‑148 (MQ=255)
                                                                                                                    cGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTAAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTaa                                >  7:380790/1‑147 (MQ=255)
                                                                                                                     gATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACg                              <  1:209033/148‑1 (MQ=255)
                                                                                                                            cGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACg                        >  4:90916/1‑147 (MQ=255)
                                                                                                                             gCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCa                      <  6:78291/148‑1 (MQ=255)
                                                                                                                                      gccAGTTCTTTGCCAGTCACCGTTTCCCCCTCAATCTACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGc             <  2:155777/148‑1 (MQ=255)
                                                                                                                                       ccAGCTGTTTTCCTGTCCCCGTTTCCCCCTAAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCg            <  2:139691/148‑1 (MQ=255)
                                                                                                                                          gCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGt         >  1:157650/1‑148 (MQ=255)
                                                                                                                                          gCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGt         >  4:126985/1‑148 (MQ=255)
                                                                                                                                               ttGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGGCAt    >  5:86789/1‑148 (MQ=255)
                                                                                                                                                  cccGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGGCATCg  <  3:33814/147‑1 (MQ=255)
                                                                                                                                                  cccGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGGCATCg  <  5:450462/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCTCATAGCCAAACAGTTCTGCTTCCAGCAGCGATTCGGCAATCGCCCCGCAGTTGACGGCAACAAACGGATGCGACTTTTTGCCCTGTCGCGCATCGTGGCGGGCAAAATATTCCCGATGAATCGCCTGGGCCGCCAGCTCTTTGCCCGTCCCCGTTTCCCCCTCAATCAACACCGCCGCACTGGAGCGGGCATACAGCAAAATAGTCTGCCGTACTTGTTCCATCTGTGGTGATTGACCGAGCATATCGCCCAGCACGTAACGAGTACGCAGGGCGTTGCGGGTGGCATCG  >  NC_000913/347546‑347838

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