Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,689,819 A→T 15.0% intergenic (+68/‑33) manA → / → ydgA mannose‑6‑phosphate isomerase/DUF945 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,689,8190AT15.0% 487.5 / 68.5 186intergenic (+68/‑33)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base A (81/77);  new base T (17/11);  total (98/88)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

Only 100 of 187 total aligned reads displayed.

CCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAT  >  NC_000913/1689721‑1689918
                                                                                                  |                                                                                                   
ccGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAg                                                                                                    >  2:1786595/1‑100 (MQ=255)
ccGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAg                                                                                                    >  1:1589477/1‑100 (MQ=255)
  gTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTc                                                                                                  >  1:1908216/1‑100 (MQ=255)
    ttAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCg                                                                                                <  1:2270103/100‑1 (MQ=255)
     tAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCCAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGc                                                                                               <  1:70213/100‑1 (MQ=255)
       gcgcGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCt                                                                                             >  2:210761/1‑100 (MQ=255)
        cgcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTg                                                                                            <  1:464466/100‑1 (MQ=255)
         gcgTGTTTACAACAAGCTGTAAGAGCTTTCTGAAAAAATTATCATCTCTTGCTAAGCTTATTCAAGGCTTATAACACCTTCATGCGGCCAGTCCGCCTg                                                                                            >  2:772209/1‑99 (MQ=255)
           gtggTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTAGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGACCGCCTGaag                                                                                         <  2:2876682/100‑3 (MQ=255)
            tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGAttt                                                                                        >  1:2192905/1‑100 (MQ=255)
            tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaagg                                                                                        <  1:665141/100‑4 (MQ=255)
             gtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTc                                                                                       >  2:87975/1‑100 (MQ=255)
             gtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaaggt                                                                                       <  1:2012281/100‑5 (MQ=255)
              tttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaa                                                                                          <  2:3613283/96‑2 (MQ=255)
               ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAt                                                                                     <  2:913862/100‑1 (MQ=255)
                tACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtt                                                                                    >  2:2286173/1‑100 (MQ=255)
                tACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtt                                                                                    <  1:2454029/100‑1 (MQ=255)
                 acaacaAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAttt                                                                                   >  1:923115/1‑100 (MQ=255)
                  caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                                  <  2:376957/100‑1 (MQ=255)
                  caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                                  >  2:3152495/1‑100 (MQ=255)
                  caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                                  >  2:248200/1‑100 (MQ=255)
                    acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                                >  1:2390598/1‑100 (MQ=255)
                    acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                                >  2:2047238/1‑100 (MQ=255)
                    acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                                >  1:3497239/1‑100 (MQ=255)
                    acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                                >  1:3119620/1‑100 (MQ=255)
                    acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                                >  2:2982647/1‑100 (MQ=255)
                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                               <  1:998271/100‑1 (MQ=255)
                     caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                               <  2:2207560/100‑1 (MQ=255)
                       aGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGa                                                                             <  1:1799363/100‑1 (MQ=255)
                        gCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGAt                                                                            >  1:2092303/1‑100 (MQ=255)
                         cTGTAAGAGCTTACCGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCCCCAGATTTCACTTTATGCATa                                                                           >  2:2655442/1‑100 (MQ=255)
                           gTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAAt                                                                         >  2:472744/1‑100 (MQ=255)
                           gTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCAGCCAGTCCGCCTGATTTCATTTTATGGATAAt                                                                         >  2:3803994/1‑100 (MQ=255)
                            tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATc                                                                        <  2:369071/100‑1 (MQ=255)
                              agagCTTACTGAAAAAATTAACATCTCTTGCTACGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAt                                                                      >  2:221907/1‑100 (MQ=255)
                              agagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAt                                                                      >  2:3928493/1‑100 (MQ=255)
                               gagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAtt                                                                     <  1:1315284/100‑1 (MQ=255)
                                agCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                    >  2:3458118/1‑100 (MQ=255)
                                agCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                    >  2:3351353/1‑100 (MQ=255)
                                     aCTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTGCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAAt                                                               <  2:296394/100‑1 (MQ=255)
                                      cTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATa                                                              >  1:3754048/1‑100 (MQ=255)
                                        gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                            >  1:1092212/1‑100 (MQ=255)
                                         aaaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAAt                                                           <  1:968215/100‑1 (MQ=255)
                                         aaaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAAt                                                           >  1:2842174/1‑100 (MQ=255)
                                         aaaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAAt                                                           >  2:1406865/1‑100 (MQ=255)
                                          aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATc                                                          <  2:3497214/100‑1 (MQ=255)
                                          aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATc                                                          <  2:3443406/100‑1 (MQ=255)
                                              aTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTg                                                      <  2:3192461/100‑1 (MQ=255)
                                               ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTgg                                                     >  1:2731526/1‑100 (MQ=255)
                                                   tcttccgatcTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGc                                                 <  2:1503176/91‑1 (MQ=255)
                                                   cATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGc                                                 <  1:2334068/100‑1 (MQ=255)
                                                     tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCgg                                               <  2:2276864/100‑1 (MQ=255)
                                                        ctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAAGCATTATGAATAAATAGGTGGTAGAGGGAg                                            >  1:3549240/1‑100 (MQ=255)
                                                         ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAgg                                           >  1:634641/1‑100 (MQ=255)
                                                         ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAgg                                           <  1:204884/100‑1 (MQ=255)
                                                         ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAgg                                           <  2:2379247/100‑1 (MQ=255)
                                                         ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAgg                                           <  1:3302519/100‑1 (MQ=255)
                                                         ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCCGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAgg                                           <  2:2667852/100‑1 (MQ=255)
                                                             tAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTc                                       <  1:2925127/100‑1 (MQ=255)
                                                               aGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                     <  1:116170/100‑1 (MQ=255)
                                                               aGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                     <  1:2912900/100‑1 (MQ=255)
                                                               aGCTTATTAAAGGCTTATAACACCTGCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                     <  2:685790/100‑1 (MQ=255)
                                                                gcatattaAAGGCTTATAACACACACAGGCGGTAAGTACACCTCATTTCATTTTATGGATAATCAATAGGAAGAAATCGCTGGTAGCGGTAGGCGTCAtt                                    <  2:1883317/97‑1 (MQ=255)
                                                                 cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttg                                   >  2:2827374/1‑100 (MQ=255)
                                                                 cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttg                                   <  2:1818880/100‑1 (MQ=255)
                                                                   tattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgtt                                 >  1:3881004/1‑100 (MQ=255)
                                                                   tattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgt                                  <  2:2053508/99‑1 (MQ=255)
                                                                   tattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCATGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgtt                                 >  2:794723/1‑100 (MQ=255)
                                                                     ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                               <  2:421676/100‑1 (MQ=255)
                                                                        aaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCt                            <  1:154676/100‑1 (MQ=255)
                                                                         aGGCTTATAACACCGTCAGGCGGTCAGTCCGCCTAATTTCATTTTCTGGATAATCATTATGAATCAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCt                            <  2:3513858/99‑1 (MQ=255)
                                                                          ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                          >  2:1127063/1‑100 (MQ=255)
                                                                          ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                          >  2:2773269/1‑100 (MQ=255)
                                                                           gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                         >  1:2569031/1‑100 (MQ=255)
                                                                           gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                         >  2:3226151/1‑100 (MQ=255)
                                                                           gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                         <  1:141481/100‑1 (MQ=255)
                                                                           gCTTATAACACCTTCAGACGGACAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                         <  2:1704520/100‑1 (MQ=255)
                                                                             ttATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                       >  2:318499/1‑100 (MQ=255)
                                                                             ttATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                       >  2:2705614/1‑100 (MQ=255)
                                                                              tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGt                      <  1:2827392/100‑1 (MQ=255)
                                                                                taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                    <  1:1840252/100‑1 (MQ=255)
                                                                                taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                    >  1:3556213/1‑100 (MQ=255)
                                                                                taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                    <  2:159345/100‑1 (MQ=255)
                                                                                taACACCCTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                    >  1:2066569/1‑100 (MQ=255)
                                                                                  acacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGtn                  <  1:3654270/100‑2 (MQ=255)
                                                                                   tgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATACATCGCTGGTAGCGGTAGGCGTCATTGTTGCTCTAGGCGTa                     >  2:1201497/6‑96 (MQ=37)
                                                                                   tgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCt                 >  2:2505313/6‑100 (MQ=38)
                                                                                    acCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                <  2:1218528/100‑1 (MQ=255)
                                                                                     ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg               >  2:3902742/1‑100 (MQ=255)
                                                                                     ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg               >  2:3258105/1‑100 (MQ=255)
                                                                                     aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTGGTCTgg               >  1:1713736/4‑100 (MQ=255)
                                                                                     aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg               >  1:3230819/4‑100 (MQ=255)
                                                                                       atCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGTCGTCATTGTTGCGCTAGGCGTAGTCGGGAc             >  2:3390336/2‑100 (MQ=255)
                                                                                        tCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGCCa            >  2:1259104/1‑100 (MQ=255)
                                                                                        tCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCCGGACa            >  2:552967/1‑100 (MQ=255)
                                                                                           ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc         >  1:659962/1‑100 (MQ=255)
                                                                                           ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc         <  2:3089166/100‑1 (MQ=255)
                                                                                             cggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcgg       <  1:1058935/100‑1 (MQ=255)
                                                                                               gcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg     >  1:1856978/1‑100 (MQ=255)
                                                                                                  aGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAt  <  2:1625274/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAT  >  NC_000913/1689721‑1689918

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