Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,725,445 G→A 5.6% A63V (GCG→GTG)  ydhF ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,725,4450GA 5.6% 226.6 / 2.1 71A63V (GCG→GTG) ydhFputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base A (2/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCCATCACAAAACGGGA  >  NC_000913/1725308‑1725590
                                                                                                                                         |                                                                                                                                                 
tctTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCAc                                                                                                                                        >  7:283320/1‑149 (MQ=255)
 ctTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACt                                                                                                                                       >  2:34253/1‑149 (MQ=255)
    aaTGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGa                                                                                                                                     >  5:181394/1‑148 (MQ=255)
    aaTGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGa                                                                                                                                     >  2:183815/1‑148 (MQ=255)
     ttgatgTGATCGCGGTGTGTGATGTAATGACCATGGACTTTTTGTGCACGCGCGGTCGTCGCGATACCGAAGTTACAGCCGATTTCCATCCGTTAAGCCATGGGCGGTGCGAGTTGCAGTGCGTCGCCACACGCCGCTTCggagggg                                                                                                                                     >  6:366810/2‑141 (MQ=255)
     atgatgTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATa                                                                                                                                   <  4:282997/149‑1 (MQ=255)
     atgatgTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATa                                                                                                                                   <  3:322765/149‑1 (MQ=255)
       gatgTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGc                                                                                                                                 >  5:247420/1‑149 (MQ=255)
        atgTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGc                                                                                                                                 >  6:75163/1‑148 (MQ=255)
         tgtgATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGcca                                                                                                                               >  8:195846/1‑149 (MQ=255)
         tgtgATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGcca                                                                                                                               <  6:181394/149‑1 (MQ=255)
             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTTCCTCGCCAAACACCGCTTCGCACTGATAGccacca                                                                                                                            >  3:100796/1‑148 (MQ=255)
                   gTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGccacca                                                                                                                            <  8:347425/142‑1 (MQ=255)
                   gTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGccacca                                                                                                                            >  7:347425/1‑142 (MQ=255)
                      aGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCa                                                                                                                   <  5:75163/148‑1 (MQ=255)
                      aGTGATGTAATGACCAATGACGTGTTCTTCACGCGCGGTCGTCGCGTTACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCa                                                                                                                   <  2:68534/148‑1 (MQ=255)
                         gATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCAt                                                                                                               >  2:276384/1‑149 (MQ=255)
                          aTGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATg                                                                                                              <  7:195846/149‑1 (MQ=255)
                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCa                                                                                                         >  1:275641/1‑149 (MQ=255)
                                 gACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATg                                                                                                              >  1:315549/1‑142 (MQ=255)
                                 gACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATg                                                                                                              <  2:315549/142‑1 (MQ=255)
                                  aCCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACgg                                                                                                      <  6:247420/149‑1 (MQ=255)
                                    cAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACggtg                                                                                                    <  3:9650/149‑1 (MQ=255)
                                     aaTGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTATCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACggtgg                                                                                                   <  8:308595/149‑1 (MQ=255)
                                                cttcACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGcc                                                                                            >  5:103478/1‑145 (MQ=255)
                                                 ttcACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGcc                                                                                            >  2:186725/1‑144 (MQ=255)
                                                 ttcACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACACCGCTTCGCACTGCTAGCCACCATAAATATCAGCATGGTCCACggtg                                                                                                    >  2:195589/1‑136 (MQ=255)
                                                 ttcACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACACCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACggtg                                                                                                    <  1:195589/136‑1 (MQ=255)
                                                    aCGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATc                                                                                      >  5:317187/1‑147 (MQ=255)
                                                    aCGCGCGATCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATaaa                                                                                                                        <  6:256382/113‑1 (MQ=255)
                                                    aCGCGCGATCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATaaa                                                                                                                        >  5:256382/1‑113 (MQ=255)
                                                      gcgcgGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATc                                                                                      >  6:163635/1‑145 (MQ=255)
                                                      gcgcgGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGa                                                                                  <  6:317187/149‑1 (MQ=255)
                                                         cgGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGc                                                                               <  5:163635/149‑1 (MQ=255)
                                                           gtcgtcGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGctc                                                                             <  2:275641/149‑1 (MQ=255)
                                                              gtcgCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCtt                                                                           >  8:79208/1‑148 (MQ=255)
                                                                     tACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATaaaa                                                                   >  8:306180/1‑149 (MQ=255)
                                                                      aCCGCATTTACTGACGATTTCCATCCTTTCAAGCATGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACACCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGGGGTCACGCCGAGATCCAGATGCTCTTCAATaaaa                                                                   <  4:100796/148‑1 (MQ=255)
                                                                       ccGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACt                                                                 <  1:34253/149‑1 (MQ=255)
                                                                        cGCATTTACTGACGATTTCCCTCCGTTCACGCAGGTGAGGTCCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCAAGGTCCACGGTGGTCACGCCGAGATCCAGATTCTCTTCAATAAAACt                                                                 >  4:214998/1‑148 (MQ=255)
                                                                        cGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACt                                                                 >  7:161493/1‑148 (MQ=255)
                                                                            tttACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCa                                                            <  1:276384/149‑1 (MQ=255)
                                                                                 tGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGCTGCTCTTCAATAAAACTGACCAGCTgg                                                       >  2:263134/1‑149 (MQ=255)
                                                                                     gATTTCCATCCGTTCACGCAGGTGAGGGGCCAGTTTCCGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTggcggg                                                   <  3:42048/149‑1 (MQ=255)
                                                                                       tttCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGAc                                                              >  7:160554/1‑136 (MQ=255)
                                                                                       tttCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGAc                                                              <  8:160554/136‑1 (MQ=255)
                                                                                           cATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATg                                                                                <  7:25914/114‑1 (MQ=255)
                                                                                           cATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATg                                                                                >  8:25914/1‑114 (MQ=255)
                                                                                           cATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGAc                                              <  7:306180/148‑1 (MQ=255)
                                                                                           cATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACa                                             <  5:262943/149‑1 (MQ=255)
                                                                                            aTCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACa                                             >  8:196632/1‑148 (MQ=255)
                                                                                                gTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGc                                                                                                                 <  1:357210/76‑1 (MQ=255)
                                                                                                gTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGc                                                                                                                 >  2:357210/1‑76 (MQ=255)
                                                                                                gTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGAcc                                                             >  6:305801/1‑128 (MQ=255)
                                                                                                gTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGAcc                                                             <  5:305801/128‑1 (MQ=255)
                                                                                                gTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTc                                        <  6:103478/149‑1 (MQ=255)
                                                                                                    aCGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCa                                      <  7:250868/147‑1 (MQ=255)
                                                                                                        aggtgaggtgCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCAt                                >  5:221404/1‑149 (MQ=255)
                                                                                                            gaggtgCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGAc                                                              >  1:63304/1‑115 (MQ=255)
                                                                                                            gaggtgCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGAc                                                              <  2:63304/115‑1 (MQ=255)
                                                                                                             aggtgCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATc                           <  7:196632/149‑1 (MQ=255)
                                                                                                                  ccAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAg                      >  4:135766/1‑149 (MQ=255)
                                                                                                                    aGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAg                      >  6:274751/1‑147 (MQ=255)
                                                                                                                        tCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAACACTGCCCAGCTGGCGGGCGGACATAGTCCCGTCCCTCACTCGCCAGTAGcc                 >  5:253962/1‑148 (MQ=255)
                                                                                                                           gTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGt                                    >  5:263983/1‑126 (MQ=255)
                                                                                                                           gTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGt                                    <  6:263983/126‑1 (MQ=255)
                                                                                                                              ccTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCCATcac           >  1:117465/1‑148 (MQ=255)
                                                                                                                                tCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCCATCACaa         >  5:66289/1‑148 (MQ=255)
                                                                                                                                    cAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCTCCAGTAGCCCATCACAAAACgg    >  6:302981/1‑149 (MQ=255)
                                                                                                                                      aaCGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCCATCACAAAACggg   >  1:109479/1‑148 (MQ=255)
                                                                                                                                      aaCGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCCATCACAAAACGGGa  <  6:221404/149‑1 (MQ=255)
                                                                                                                                         gccgcTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGAc                                                              >  7:358191/1‑86 (MQ=255)
                                                                                                                                         gccgcTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGAc                                                              <  8:358191/86‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
TCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTCACGCAGGTGAGGTGCCAGTTTCAGTGCCTCGCCAAACGCCGCTTCGCACTGATAGCCACCATAAATATCAGCATGGTCCACGGTGGTCACGCCGAGATCCAGATGCTCTTCAATAAAACTGACCAGCTGGCGGGCGGACATATTCCAGTCCATCAATCGCCAGTAGCCCATCACAAAACGGGA  >  NC_000913/1725308‑1725590

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: