Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,971,075 2 bp→GC coding (942‑943/1068 nt) rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,0750CG99.8% 1956.9 / NA 568C314W (TGC→TGG) rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (255/311);  total (255/313)
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
*NC_0009133,971,0760GC100.0% 1961.7 / NA 567V315L (GTG→CTG) rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (254/313);  total (254/313)

Only 100 of 590 total aligned reads displayed.

TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162
                                                                                       ||                                                                                      
   ccTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 <  1:8310081/93‑1 (MQ=255)
     tGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAgg                                                                               >  2:8548814/1‑93 (MQ=255)
        tCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGaaa                                                                            <  2:20989/93‑1 (MQ=255)
         cACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAc                                                                           <  1:4659725/93‑1 (MQ=255)
         cACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAc                                                                           <  1:724738/93‑1 (MQ=255)
         cACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAc                                                                           >  1:6841738/1‑93 (MQ=255)
           ccTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGACACCt                                                                         >  2:1369595/1‑93 (MQ=255)
           ccTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCt                                                                         <  1:1599358/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  2:2786215/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  1:6369115/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  1:5980656/93‑1 (MQ=255)
               tGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGa                                                                     <  2:6499828/93‑1 (MQ=255)
                gTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaa                                                                    <  1:861445/93‑1 (MQ=255)
                 tCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaaa                                                                   >  1:3719549/1‑93 (MQ=255)
                   gCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGt                                                                 <  1:1015090/93‑1 (MQ=255)
                    cTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTg                                                                >  1:5184072/1‑93 (MQ=255)
                      gACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGa                                                              >  2:5580229/1‑93 (MQ=255)
                       gccGTCCGGGGCATGATCTGCGTCATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  2:4835480/92‑1 (MQ=255)
                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  1:4879143/93‑1 (MQ=255)
                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             >  1:5480222/1‑93 (MQ=255)
                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             >  1:1529460/1‑93 (MQ=255)
                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  2:665670/93‑1 (MQ=255)
                        ccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAt                                                            <  2:3997877/93‑1 (MQ=255)
                        ccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAt                                                            <  1:7453565/93‑1 (MQ=255)
                         cGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATg                                                           >  2:3706107/1‑93 (MQ=255)
                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGc                                                          <  2:287537/93‑1 (MQ=255)
                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGc                                                          <  1:8359847/93‑1 (MQ=255)
                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGc                                                          <  1:6080263/93‑1 (MQ=255)
                          gTCCGGGGCATAATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCATGAAACCTTTGAAAGTGGAGTGc                                                          >  2:6080548/1‑93 (MQ=255)
                          gTCCGGCGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGc                                                          <  2:1682703/93‑1 (MQ=255)
                           tCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCg                                                         <  2:2488954/93‑1 (MQ=255)
                            ccGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGt                                                        >  2:5736082/1‑93 (MQ=255)
                              ggggTATGATCTGCGTTATGCCTTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAGCCTTTGAAAGTGGAATGCGTaa                                                      <  2:4691404/93‑1 (MQ=255)
                                ggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaaaa                                                    <  1:2075544/93‑1 (MQ=255)
                                 gCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAc                                                   <  1:7093566/93‑1 (MQ=255)
                                 gCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAc                                                   <  2:1079180/93‑1 (MQ=255)
                                  gaTGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACg                                                  <  1:1914323/92‑1 (MQ=255)
                                   aTGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACgg                                                 >  2:5841704/1‑93 (MQ=255)
                                     gATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGtt                                               >  2:6915757/1‑93 (MQ=255)
                                      aTCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTc                                              >  1:7341279/1‑93 (MQ=255)
                                      aTCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTc                                              >  1:1738669/1‑93 (MQ=255)
                                       tCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCa                                             >  1:3099725/1‑93 (MQ=255)
                                       tCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCa                                             >  2:3591032/1‑93 (MQ=255)
                                       tCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCa                                             >  2:2652679/1‑93 (MQ=255)
                                        cTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAg                                            <  2:1213362/93‑1 (MQ=255)
                                         tGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAgg                                           >  1:1159209/1‑92 (MQ=255)
                                          gCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTg                                          <  1:8203511/93‑1 (MQ=255)
                                          gCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTg                                          <  1:3229889/93‑1 (MQ=255)
                                           cGTTATGCCATTGATGCTTCGAAAATTGCCCGTGATCTGGGCGGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGATTCAGTgg                                         <  1:5766197/93‑1 (MQ=255)
                                            gTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGGAAAACGGTTCAGTgta                                        >  2:8790504/1‑91 (MQ=255)
                                                tGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCt                                    <  1:953351/93‑1 (MQ=255)
                                                tGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCt                                    <  2:8195828/93‑1 (MQ=255)
                                                tGCCATTGATGCTTCGAAAAGTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCt                                    <  1:6273345/93‑1 (MQ=255)
                                                 gCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTg                                   >  1:225812/1‑93 (MQ=255)
                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 <  2:7004068/93‑1 (MQ=255)
                                                    aTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCt                                >  2:1872220/1‑93 (MQ=255)
                                                     ttGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTa                               >  2:6245898/1‑93 (MQ=255)
                                                     ttGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTa                               >  1:7774391/1‑93 (MQ=255)
                                                      tGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTaa                              >  2:2593218/1‑93 (MQ=255)
                                                       gATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAAt                             >  1:2544680/1‑93 (MQ=255)
                                                       gATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAAt                             >  2:3476965/1‑93 (MQ=255)
                                                          gCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGaa                          <  2:4162525/93‑1 (MQ=255)
                                                           cTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGaaa                         >  2:8086488/1‑93 (MQ=255)
                                                           cTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGaaa                         <  1:2683623/93‑1 (MQ=255)
                                                           cTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGaaa                         <  2:6744253/93‑1 (MQ=255)
                                                            ttCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAg                        >  2:7781923/1‑93 (MQ=255)
                                                             tCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGc                       <  2:7428997/93‑1 (MQ=255)
                                                             tCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGc                       >  2:4284062/1‑93 (MQ=255)
                                                               gAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtg                     >  1:444860/1‑93 (MQ=255)
                                                                aaaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtgg                    <  2:4906894/93‑1 (MQ=255)
                                                                 aaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggt                   <  1:7001081/93‑1 (MQ=255)
                                                                   aTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtgg                 >  2:3170413/1‑93 (MQ=255)
                                                                   aTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtgg                 >  2:4448501/1‑93 (MQ=255)
                                                                   aTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtgg                 >  1:4000722/1‑93 (MQ=255)
                                                                    ttGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGa                >  1:5049026/1‑93 (MQ=255)
                                                                    ttGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGa                >  2:6343277/1‑93 (MQ=255)
                                                                     tGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGaa               <  1:875617/93‑1 (MQ=255)
                                                                     tGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGaa               <  2:5929847/93‑1 (MQ=255)
                                                                     tGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGaa               <  2:7966421/93‑1 (MQ=255)
                                                                       cccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGc             <  2:8657584/93‑1 (MQ=255)
                                                                       cccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGc             <  1:5602736/93‑1 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            <  1:3834319/93‑1 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            <  2:2206073/93‑1 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            >  1:5366045/1‑93 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTCGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            <  2:3303129/93‑1 (MQ=255)
                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          >  2:8323050/1‑93 (MQ=255)
                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          >  1:4796860/1‑93 (MQ=255)
                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          <  1:3362756/93‑1 (MQ=255)
                                                                          ctGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          <  1:1268292/92‑1 (MQ=255)
                                                                           tGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGt         <  1:6593353/93‑1 (MQ=255)
                                                                             atcTTGGCGGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGc       <  2:7509324/91‑1 (MQ=255)
                                                                             aaCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGc       >  2:7263303/1‑93 (MQ=255)
                                                                             aaCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGc       >  1:7992336/1‑93 (MQ=255)
                                                                             aaCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGc       <  2:2317038/93‑1 (MQ=255)
                                                                              aCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATTTGGCTAATGAAAGCTGGTGGAAGCAGGTGCa      <  2:383102/93‑1 (MQ=255)
                                                                              aCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCa      >  1:4260682/1‑93 (MQ=255)
                                                                               cTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAg     >  2:5033776/1‑93 (MQ=255)
                                                                                ttGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAgg    >  1:3566038/1‑93 (MQ=255)
                                                                                 tggctggcTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGa   <  1:2889845/93‑1 (MQ=255)
                                                                                  ggctggctGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAc  <  2:5155741/93‑1 (MQ=255)
                                                                                       ||                                                                                      
TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑