Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,259,839 G→A L183L (CTG→TTG)  fruK ← fructose‑1‑phosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,8390GA99.6% 825.8 / ‑5.5 229L183L (CTG→TTG) fruKfructose‑1‑phosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (106/122);  total (107/122)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

Only 100 of 250 total aligned reads displayed.

GTAGCGCATGTGCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCC  >  NC_000913/2259749‑2259928
                                                                                          |                                                                                         
  aGCGCATGTGCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCAc                                                                                       <  2:1389149/93‑1 (MQ=255)
    cgcATGTGCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACgg                                                                                     <  1:728806/93‑1 (MQ=255)
     gcATGTGCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGt                                                                                    >  1:2779793/1‑93 (MQ=255)
     gcATGTGCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGt                                                                                    >  2:2331839/1‑93 (MQ=255)
     gcATGTGCAGCTTCAATCACATCTTTCACTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGt                                                                                    <  2:6447116/93‑1 (MQ=255)
       aTGTGCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGc                                                                                  >  2:2700030/1‑93 (MQ=255)
       aTGTGCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGc                                                                                  >  2:559336/1‑93 (MQ=255)
         gtgCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCg                                                                                >  2:5906263/1‑93 (MQ=255)
           gCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCt                                                                              >  2:4711027/1‑93 (MQ=255)
           gCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCt                                                                              >  2:3721076/1‑93 (MQ=255)
             aGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCttt                                                                            >  2:4522289/1‑93 (MQ=255)
             aGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCttt                                                                            <  2:448075/93‑1 (MQ=255)
             aGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACCGTGCCGCttt                                                                            >  1:2010568/1‑93 (MQ=255)
              gCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGGTAGGGTTCACCAACCACGGTGCCGCTTTc                                                                           >  2:1078354/1‑93 (MQ=255)
                  cAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAAc                                                                       <  2:5859971/93‑1 (MQ=255)
                   aaTCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCaaagc                                                                      <  1:4783936/93‑3 (MQ=255)
                   aaTCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGATAGGCTTCACCAACCACGGTGCCGCTTTCACAcc                                                                      >  1:1045070/1‑93 (MQ=255)
                    aTCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCt                                                                     >  1:2729408/1‑93 (MQ=255)
                    aTCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCt                                                                     >  1:4082828/1‑93 (MQ=255)
                     tCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTg                                                                    >  1:6549161/1‑93 (MQ=255)
                      cacaTCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGc                                                                   <  1:3018578/93‑1 (MQ=255)
                       acaTCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGct                                                                  >  2:3179817/1‑93 (MQ=255)
                        caTCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGcta                                                                 >  1:1804150/1‑93 (MQ=255)
                         aTCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGctac                                                                >  1:6407149/1‑93 (MQ=255)
                         aTCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGctac                                                                <  2:5474472/93‑1 (MQ=255)
                         aTCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGctac                                                                <  2:2256455/93‑1 (MQ=255)
                          tCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGctact                                                               >  2:3584574/1‑93 (MQ=255)
                            tttcatttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTaa                                                             <  1:2967319/93‑1 (MQ=255)
                            tttcatttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTaa                                                             >  2:3478179/1‑93 (MQ=255)
                              tcatttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGGGGCGGGTAGGATTCACCAACCACGGTGCCGCTTTCAAACCAGCTACTAAca                                                           >  1:3287161/1‑92 (MQ=255)
                              tcatttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACg                                                           <  1:5076629/93‑1 (MQ=255)
                               catttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGc                                                          <  1:5584349/93‑1 (MQ=255)
                                atttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCt                                                         >  1:4179449/1‑93 (MQ=255)
                                 tttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCtt                                                        >  2:4924679/1‑93 (MQ=255)
                                  ttcaGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTc                                                       <  1:1353364/93‑1 (MQ=255)
                                    caGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCAc                                                     <  1:914211/93‑1 (MQ=255)
                                     aGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACg                                                    >  2:159182/1‑93 (MQ=255)
                                     aGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACg                                                    >  2:1541195/1‑93 (MQ=255)
                                     aGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACg                                                    >  1:2341239/1‑93 (MQ=255)
                                     aGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGCTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACg                                                    >  1:2152485/1‑93 (MQ=255)
                                      ggCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACgg                                                   >  2:4175344/1‑93 (MQ=255)
                                       gCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGc                                                  <  1:1574632/93‑1 (MQ=255)
                                       gCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGc                                                  <  2:2311709/93‑1 (MQ=255)
                                       gCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGATTCACCTACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACCGc                                                  >  2:3695050/1‑93 (MQ=255)
                                        cAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGct                                                 <  1:4465180/93‑1 (MQ=255)
                                        cAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGGTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGct                                                 >  2:3313892/1‑93 (MQ=255)
                                         aGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGcta                                                >  1:3463189/1‑93 (MQ=255)
                                          gTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACTGTGCCGCCTTCAAACCTGCTACTAACGCTTCACGGctac                                               >  2:505508/1‑93 (MQ=255)
                                          gTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGctac                                               <  2:6022386/93‑1 (MQ=255)
                                          gTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGctac                                               <  1:6221116/93‑1 (MQ=255)
                                           tttACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGctact                                              >  2:1681877/1‑93 (MQ=255)
                                             tACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTCCTAACGCTTCACGGCTACTAt                                            <  1:6130654/93‑1 (MQ=255)
                                             tACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTAt                                            <  2:2869721/93‑1 (MQ=255)
                                              aCGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATc                                           <  1:1315153/93‑1 (MQ=255)
                                              aCGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATc                                           >  1:1005820/1‑93 (MQ=255)
                                                ggCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCaa                                         >  1:2064922/1‑93 (MQ=255)
                                                ggCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCaa                                         <  2:2964343/93‑1 (MQ=255)
                                                ggCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCaa                                         >  2:1324610/1‑93 (MQ=255)
                                                 gCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCcaa                                        >  1:1160218/1‑93 (MQ=255)
                                                 gCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCaaa                                        <  2:922395/93‑1 (MQ=255)
                                                 gCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCaaa                                        <  1:2509112/93‑1 (MQ=255)
                                                 gCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCaaa                                        >  2:771444/1‑93 (MQ=255)
                                                    tGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGAt                                     <  2:5203774/93‑1 (MQ=255)
                                                     gCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATa                                    <  1:6004620/93‑1 (MQ=255)
                                                      cccAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATaa                                   <  1:5410997/93‑1 (MQ=255)
                                                      cccAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATaa                                   <  1:4968752/93‑1 (MQ=255)
                                                       ccAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAAt                                  >  2:2791301/1‑93 (MQ=255)
                                                       ccAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAAt                                  >  2:1309869/1‑93 (MQ=255)
                                                         aGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGc                                <  2:1736114/93‑1 (MQ=255)
                                                            tttCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAg                             <  2:4328287/93‑1 (MQ=255)
                                                            tttCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAg                             <  1:1210956/93‑1 (MQ=255)
                                                            tttCCAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAg                             <  2:6250519/93‑1 (MQ=255)
                                                             ttCCAGCTCGCGGCGGTTATGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACTCTTCACGGCTACTATCAAAGATAATGCaatg                            <  2:4043698/93‑3 (MQ=255)
                                                                cAGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGaca                         >  2:2401289/1‑93 (MQ=255)
                                                                 aGCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGacac                        >  2:3389818/1‑93 (MQ=255)
                                                                  gCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACt                       >  2:3081756/1‑93 (MQ=255)
                                                                  gCTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACt                       <  2:5998973/93‑1 (MQ=255)
                                                                   cTCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTg                      >  2:6559216/1‑93 (MQ=255)
                                                                    tCGCGGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGa                     >  2:2754208/1‑93 (MQ=255)
                                                                    ccgcgGCGGTTAGGTTTCACCAACGACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGa                     <  1:3546711/92‑1 (MQ=255)
                                                                     cgcgGCGGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGAc                    <  2:5271598/93‑1 (MQ=255)
                                                                      gcggcgGTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACt                   <  1:520632/93‑1 (MQ=255)
                                                                      gccgcggTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACt                   <  2:2143349/90‑1 (MQ=255)
                                                                       cggcggTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTa                  <  2:3531428/93‑1 (MQ=255)
                                                                       cggcggTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTa                  >  2:3919666/1‑93 (MQ=255)
                                                                          cggTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGc               >  2:1121850/1‑93 (MQ=255)
                                                                           gtttAGGTTTCACCACCCACGGTGCCGCTGTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCa              <  1:1454638/91‑1 (MQ=255)
                                                                           ggTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCa              >  2:2455922/1‑93 (MQ=255)
                                                                           ggTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCa              >  2:172623/1‑93 (MQ=255)
                                                                            gTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGGTACTATCAAAGATAATGCAAGGACACTGACTACGCAg             >  1:4420112/1‑93 (MQ=255)
                                                                            gTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAg             >  1:4887156/1‑93 (MQ=255)
                                                                            gTTAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAg             >  1:6693915/1‑93 (MQ=255)
                                                                             ttAGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAgg            <  2:5226391/93‑1 (MQ=255)
                                                                               aGGTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAggtg          <  1:137790/93‑3 (MQ=255)
                                                                                ggTTTCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTCACTACGCAGGTGa         >  1:5515367/1‑93 (MQ=255)
                                                                                  tttCACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGt       >  2:6675914/1‑93 (MQ=255)
                                                                                   ttCACCAACCACGGTGCCGGTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTc      >  1:6475947/1‑93 (MQ=255)
                                                                                     cACCAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCAt    <  1:6345780/93‑1 (MQ=255)
                                                                                      accaaccaCGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATc   <  1:3690938/93‑1 (MQ=255)
                                                                                       ccAACCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATcc  >  1:1694702/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTAGCGCATGTGCAGCTTCAATCACATCTTTCATTTCAGGCAGTTTACGGCCTGCCCAGATTTCCAGCTCGCGGCGGTTAGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCC  >  NC_000913/2259749‑2259928

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