Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,689,817 C→A 16.5% intergenic (+66/‑35) 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,8170CA16.5% 466.1 / 70.7 181intergenic (+66/‑35)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base C (78/71);  new base A (18/12);  total (96/85)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

CGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCG  >  NC_000913/1689718‑1689915
                                                                                                   |                                                                                                  
cGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAggcggc                                                                                                    <  2:903160/100‑1 (MQ=255)
cGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAggcggc                                                                                                    <  1:3837529/100‑1 (MQ=255)
   ccGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAg                                                                                                 >  1:1589477/1‑100 (MQ=255)
     gTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTc                                                                                               >  1:1908216/1‑100 (MQ=255)
      tttAGCGCGTGTTTACAACAATCTGTAAAAAAAAAAAAAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGcc                                                                                                   <  2:2775275/95‑1 (MQ=255)
        tAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGc                                                                                            <  1:2905087/100‑1 (MQ=255)
          gcgcGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCt                                                                                          >  2:1685225/1‑100 (MQ=255)
           cgcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTg                                                                                         >  2:1507294/1‑100 (MQ=255)
           cgcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTg                                                                                         <  1:464466/100‑1 (MQ=255)
           cgcgTGTTTACAACAAGCTGTAAAAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTg                                                                                         <  2:2007708/100‑1 (MQ=255)
            gcgTGTTTACAACAAGCTGTAAGAGATAACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGa                                                                                        <  2:736702/100‑1 (MQ=255)
             cgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGAt                                                                                       <  2:3791780/100‑1 (MQ=255)
             cgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGAt                                                                                       <  2:2830415/100‑1 (MQ=255)
             cgTGTTTACAACAAGCTGTAAGAAATAAATAAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGcc                                                                                           <  2:3213102/96‑1 (MQ=255)
              gtgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaag                                                                                      <  2:646625/100‑3 (MQ=255)
              gtggTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTAGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGACCGCCTGaag                                                                                      <  2:2876682/100‑3 (MQ=255)
               tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaagg                                                                                     <  1:2407817/100‑4 (MQ=255)
                gtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTc                                                                                    >  1:3626948/1‑100 (MQ=255)
                gtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTc                                                                                    >  2:2582764/1‑100 (MQ=255)
                gtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaaggt                                                                                    <  1:851516/100‑5 (MQ=255)
                gtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaaggt                                                                                    <  1:2012281/100‑5 (MQ=255)
                gtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaaggt                                                                                    <  1:247249/100‑5 (MQ=255)
                 tttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCa                                                                                   >  2:204884/1‑100 (MQ=255)
                  ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAt                                                                                  <  2:913862/100‑1 (MQ=255)
                    acaacaAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAttt                                                                                >  1:923115/1‑100 (MQ=255)
                     caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                               >  2:3152495/1‑100 (MQ=255)
                     caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAAGCTTATAACACCTTCAGGCGGCCAGTACGCCTGATTTCAtttt                                                                               >  2:793293/1‑100 (MQ=255)
                     caacaaGCTGCAAGACCTTACTGAAACAATTAACACCTCCTGCTAAGCGTATTAACGGCTTATAACACCTTCAGTCGGCCAGTCCGCCTGATTTCGttt                                                                                <  2:2216568/99‑1 (MQ=255)
                      aacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                               >  2:1707591/1‑99 (MQ=255)
                       acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTan                                                                             <  2:543111/100‑2 (MQ=255)
                       acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                             >  2:2982647/1‑100 (MQ=255)
                        caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                            <  1:998271/100‑1 (MQ=255)
                           gCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGGTTATTAAAGGATTATAACACCTTCAGGCGGCCAGTCCGCATGATTTCATTTTATGGAt                                                                         >  2:2271491/1‑100 (MQ=255)
                           gCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGAt                                                                         >  1:2092303/1‑100 (MQ=255)
                           gCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGAt                                                                         <  2:276691/100‑1 (MQ=255)
                            cTGTAAGAGCTTACTGAAAAAATTAACCTCTCTTGCTAAGCTGATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATa                                                                        <  2:2169178/100‑1 (MQ=255)
                            cTGTAAGAGCTTACTGAAAAAATTAACAGCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATa                                                                        <  2:1365522/100‑1 (MQ=255)
                                 agagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAt                                                                   >  2:3928493/1‑100 (MQ=255)
                                 agagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAt                                                                   <  2:2087555/100‑1 (MQ=255)
                                  gagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAtt                                                                  <  1:1315284/100‑1 (MQ=255)
                                   agattaaaaaaaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                 <  2:3925913/91‑1 (MQ=255)
                                   agCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                 <  2:2658734/100‑1 (MQ=255)
                                   agCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTa                                                                 >  2:3458118/1‑100 (MQ=255)
                                    gCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                <  1:2340048/100‑1 (MQ=255)
                                       tagtGAAAAAATTAACATCTCTTGCTAAGATTATTAAAGGCTTATTACATCTTCAGGGGGCCAGTCCGCCTGATTTGATTTTATGGATAGTCATTATg                                                               >  2:3259300/4‑98 (MQ=255)
                                        aCTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTGCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAAt                                                            <  2:296394/100‑1 (MQ=255)
                                           gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                         >  2:2355537/1‑100 (MQ=255)
                                           gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                         >  2:2461172/1‑100 (MQ=255)
                                           gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                         >  1:1092212/1‑100 (MQ=255)
                                            aaaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAAt                                                        <  1:968215/100‑1 (MQ=255)
                                             aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATc                                                       <  2:3497214/100‑1 (MQ=255)
                                             aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATc                                                       <  2:3443406/100‑1 (MQ=255)
                                             aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATc                                                       <  1:2297163/100‑1 (MQ=255)
                                             aaaaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAAAACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAACTc                                                       >  2:3752804/1‑100 (MQ=255)
                                                  ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTgg                                                  <  1:1803740/100‑1 (MQ=255)
                                                  ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTgg                                                  >  1:2731526/1‑100 (MQ=255)
                                                  ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTgg                                                  >  1:3253525/1‑100 (MQ=255)
                                                    aaCATCTCTTGCTAAGCTTATTAAAGTCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGTTa                                                <  2:1594666/100‑1 (MQ=255)
                                                    aaCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTa                                                >  1:3456348/1‑100 (MQ=255)
                                                        tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCgg                                            <  2:2276864/100‑1 (MQ=255)
                                                        tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCgg                                            <  1:694171/100‑1 (MQ=255)
                                                           ctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAg                                         >  1:2198907/1‑100 (MQ=255)
                                                           ctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAAGCATTATGAATAAATAGGTGGTAGAGGGAg                                         >  1:3549240/1‑100 (MQ=255)
                                                            ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAgg                                        >  1:634641/1‑100 (MQ=255)
                                                            ttGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAgg                                        <  1:204884/100‑1 (MQ=255)
                                                              gCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCg                                      <  1:292131/100‑1 (MQ=255)
                                                                  aGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                  <  1:116170/100‑1 (MQ=255)
                                                                   gcgccttaAAGGCTGATAACACCTGCTGGGCGCCATTCCTCCTGCTTTCATTTTATGGATAATCATTAGGCATATGCCGCTGGTAGCGGTAGGCg                                      >  2:2702511/6‑95 (MQ=255)
                                                                   gcatattaAAGGCTTATAACACACACAGGCGGTAAGTACACCTCATTTCATTTTATGGATAATCAATAGGAAGAAATCGCTGGTAGCGGTAGGCGTCAtt                                 <  2:1883317/97‑1 (MQ=255)
                                                                    cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttg                                <  2:1818880/100‑1 (MQ=255)
                                                                      tattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgtt                              >  1:3881004/1‑100 (MQ=255)
                                                                       attaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgttg                             <  1:1905934/100‑1 (MQ=255)
                                                                        ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                            <  2:421676/100‑1 (MQ=255)
                                                                          aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                          >  1:3555971/1‑100 (MQ=255)
                                                                           aaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCt                         <  1:154676/100‑1 (MQ=255)
                                                                             ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                       >  2:1127063/1‑100 (MQ=255)
                                                                              gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                      >  1:2569031/1‑100 (MQ=255)
                                                                              gCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                      <  1:141481/100‑1 (MQ=255)
                                                                                ttATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                    >  2:318499/1‑100 (MQ=255)
                                                                                ttATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                    <  2:2482958/100‑1 (MQ=255)
                                                                                 tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGt                   <  1:2827392/100‑1 (MQ=255)
                                                                                   taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                 <  2:3242822/100‑1 (MQ=255)
                                                                                   taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                 <  1:3675522/100‑1 (MQ=255)
                                                                                   taACACCCTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                 >  1:2066569/1‑100 (MQ=255)
                                                                                   aaatgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAg                 >  2:3582823/9‑100 (MQ=37)
                                                                                     acacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTc               <  1:318499/100‑1 (MQ=255)
                                                                                      tgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATACATCGCTGGTAGCGGTAGGCGTCATTGTTGCTCTAGGCGTa                  >  2:1201497/6‑96 (MQ=37)
                                                                                       gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg             >  2:1423687/5‑100 (MQ=255)
                                                                                       acCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg             <  2:1218528/100‑1 (MQ=255)
                                                                                       acCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg             <  1:3045932/100‑1 (MQ=255)
                                                                                        aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTGGTCTgg            >  1:1713736/4‑100 (MQ=255)
                                                                                        aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg            >  2:3626598/4‑100 (MQ=255)
                                                                                        aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTAGGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGATAGGCGTAGTCTgg            >  2:3812465/4‑100 (MQ=255)
                                                                                          atCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGTCGTCATTGTTGCGCTAGGCGTAGTCGGGAc          >  2:3390336/2‑100 (MQ=255)
                                                                                          atCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc          >  2:851174/2‑100 (MQ=255)
                                                                                           tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa         >  2:3382305/1‑100 (MQ=255)
                                                                                           tCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa         >  1:1853909/1‑100 (MQ=255)
                                                                                              ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc      <  2:3089166/100‑1 (MQ=255)
                                                                                              ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc      >  1:659962/1‑100 (MQ=255)
                                                                                                cggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcgg    <  1:1058935/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCG  >  NC_000913/1689718‑1689915

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