Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,689,817 C→A 15.8% 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,8170CA15.8% 353.4 / 48.7 132intergenic (+66/‑35)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base C (58/52);  new base A (11/10);  total (69/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 148 total aligned reads displayed.

CGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGC  >  NC_000913/1689718‑1689916
                                                                                                   |                                                                                                   
cGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAggcggc                                                                                                     >  2:81728/1‑100 (MQ=255)
      tttAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTcc                                                                                               <  1:3296061/100‑1 (MQ=255)
         aGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGcc                                                                                            <  2:3907902/100‑1 (MQ=255)
           cgcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTg                                                                                          >  2:280431/1‑100 (MQ=255)
           cgcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTg                                                                                          <  2:3857609/100‑1 (MQ=255)
           cgcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTg                                                                                          <  1:649054/100‑1 (MQ=255)
           cgcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCCCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTg                                                                                          <  2:1029217/100‑1 (MQ=255)
            gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGa                                                                                         <  2:1751233/100‑1 (MQ=255)
            gcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGa                                                                                         >  1:2218745/1‑100 (MQ=255)
             cgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGAt                                                                                        >  1:2823496/1‑100 (MQ=255)
              gtgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaag                                                                                       <  2:3475090/100‑3 (MQ=255)
               tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGACTGGCCGCCTGaagg                                                                                      <  2:3182312/100‑4 (MQ=255)
                 tttACAACAAGCTGTAAGAGCTTACTGAAAAAGTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCa                                                                                    >  1:4017002/1‑100 (MQ=255)
                 tttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCATATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCa                                                                                    >  1:329366/1‑100 (MQ=255)
                  ttACAACAAGCTTTAAGAGCTTACTGAAAAAATTAACTTCTCTTGGTAAGCTCATTAACGGAGTATAACACCTTCAGGCGGACTGGCCGCCTGaaggtgt                                                                                   <  1:4094994/100‑7 (MQ=19)
                   tACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtt                                                                                  >  2:2038299/1‑100 (MQ=255)
                   tACAACAAGCTGTAAGAGCTTACTGAAAAAATAAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtt                                                                                  >  2:4179075/1‑100 (MQ=255)
                    acaacaAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAttt                                                                                 >  2:736966/1‑100 (MQ=255)
                     caacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAtttt                                                                                >  2:332420/1‑100 (MQ=255)
                      aacaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTa                                                                               <  2:4020754/100‑1 (MQ=255)
                       acaaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTAt                                                                              >  2:1209766/1‑100 (MQ=255)
                        caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGCTTTCATTCTATg                                                                             <  2:2893262/100‑1 (MQ=255)
                        caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                             <  1:2494630/100‑1 (MQ=255)
                         atgCTGTAAGAGATAAATGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATg                                                                             <  1:1358395/97‑1 (MQ=255)
                          aGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACAACTTCAAGCGGCCAGACCGCCTGATTTCATTTTATGGa                                                                           >  1:2963104/1‑100 (MQ=255)
                             tGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATaa                                                                        <  1:1097992/100‑1 (MQ=255)
                              gTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAAt                                                                       >  1:1694377/1‑100 (MQ=255)
                              gTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAAt                                                                       >  2:1094184/1‑100 (MQ=255)
                               tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATc                                                                      >  2:243669/1‑100 (MQ=255)
                                 agagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCAt                                                                    >  2:2447244/1‑100 (MQ=255)
                                   agCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGACTTATAACACCTTCAGGCGTCCAGTCCGCCTGATTTCATATTATGGATAATCATTa                                                                  >  1:2563620/1‑100 (MQ=255)
                                    gCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTAt                                                                 <  1:695334/100‑1 (MQ=255)
                                       taatGACAAAATTAACATCTCTTGCTAAGCTTATTACAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGa                                                               <  2:200643/96‑1 (MQ=255)
                                       tACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGaa                                                              <  1:3786970/100‑1 (MQ=255)
                                        aCTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAAt                                                             <  2:2222520/100‑1 (MQ=255)
                                         cTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATa                                                            >  1:4345394/1‑100 (MQ=255)
                                          tGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGATTATAACACCTTCAGCCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaa                                                           >  1:3990127/1‑100 (MQ=255)
                                           gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                          <  2:1221068/100‑1 (MQ=255)
                                           gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                          >  2:1929447/1‑100 (MQ=255)
                                           gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                          >  1:1844315/1‑100 (MQ=255)
                                           gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                          >  1:1307468/1‑100 (MQ=255)
                                           gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATaaa                                                          >  1:3067470/1‑100 (MQ=255)
                                            aacaaaTTAACATCTCTTGCTAAGCTTATTAACGGCTTATAACACATTCAGACGGATATTCCACCTGATTTCAGTTTATGGATAATCATTATGAATAAAt                                                         <  1:2963443/97‑1 (MQ=255)
                                               aaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGc                                                      <  2:3386888/100‑1 (MQ=255)
                                               aaaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGc                                                      <  2:2580167/100‑1 (MQ=255)
                                                aaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCt                                                     <  2:2301956/100‑1 (MQ=255)
                                                aaTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCt                                                     <  2:1074132/100‑1 (MQ=255)
                                                 aTTAACCTATCTTGCTAAGCTTATTAAAGGCTTATCACAACTCCAGGCGGTCAGTCCGCCCGATTTCATTTTATGGATAATCATTATGAATAAATCGc                                                      <  1:1461227/98‑1 (MQ=255)
                                                       aTCTCTTGCTAAGCTTATTAAAGTCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCGTTTTATGGATAATCATTATGAATAAATCGCTGGTAGCg                                              <  1:876864/100‑1 (MQ=255)
                                                        tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCgg                                             <  2:2635492/100‑1 (MQ=255)
                                                        tctcTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCgg                                             <  1:1176287/100‑1 (MQ=255)
                                                             tGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGc                                        >  1:978160/1‑100 (MQ=255)
                                                              gCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCg                                       <  1:3175069/100‑1 (MQ=255)
                                                              gCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCg                                       <  2:329367/100‑1 (MQ=255)
                                                                  aGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                   >  2:3570850/1‑100 (MQ=255)
                                                                   gCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAtt                                  >  2:4359363/1‑100 (MQ=255)
                                                                    cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttg                                 >  2:3153865/1‑100 (MQ=255)
                                                                    cTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttg                                 <  2:418158/100‑1 (MQ=255)
                                                                     ttattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgt                                >  1:2385823/1‑100 (MQ=255)
                                                                      tattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgtt                               >  2:189949/1‑100 (MQ=255)
                                                                        ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgttg                              >  1:3424723/1‑99 (MQ=255)
                                                                        ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                             <  2:1138288/100‑1 (MQ=255)
                                                                        ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                             <  2:2400069/100‑1 (MQ=255)
                                                                         taAAGGCTTATAACACCTTCAGTCGGCCAGTCCGCCTGATTTCATGTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcg                            <  1:1208077/100‑1 (MQ=255)
                                                                         taAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcg                            <  2:2218754/100‑1 (MQ=255)
                                                                         taAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttgttgag                            >  1:2374942/1‑98 (MQ=255)
                                                                          aaaGGCTTATAACACCTTCAGGCGGCCAGTCCTCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                           <  1:286229/100‑1 (MQ=255)
                                                                          aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                           >  2:1923438/1‑100 (MQ=255)
                                                                          aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                           <  2:276267/100‑1 (MQ=255)
                                                                           aaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgcgc                           >  1:1717761/1‑99 (MQ=255)
                                                                           aaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCt                          >  1:1313958/1‑100 (MQ=255)
                                                                           aaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCt                          >  1:432033/1‑100 (MQ=255)
                                                                             ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                        >  1:3663116/1‑100 (MQ=255)
                                                                             ggCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAg                        >  2:1643977/1‑100 (MQ=255)
                                                                               cgtattACACCTTCAGGCGGCCAGTCCGCCTGATTTTCTTTTATGGCTGATCATTTTGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                       <  1:1864109/97‑1 (MQ=255)
                                                                                  aaaatgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTa                   >  1:1990871/10‑100 (MQ=37)
                                                                                    aatgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGt                 >  2:675838/8‑100 (MQ=38)
                                                                                    aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGt                 >  2:180766/1‑100 (MQ=255)
                                                                                     atgaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTc                >  2:642685/7‑100 (MQ=38)
                                                                                       gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg              >  1:3851106/5‑100 (MQ=255)
                                                                                       gaaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGACATTGTTGCGCTAGGCGTAGTCt               >  2:1191412/5‑99 (MQ=38)
                                                                                        tcTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCACTGGTAGCGGTATGTGTCATTGTTGCGCTAGGCGTAGTCTgg             >  1:755757/2‑100 (MQ=255)
                                                                                        ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg             <  2:1069367/100‑1 (MQ=255)
                                                                                        ccTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg             >  1:1029020/1‑100 (MQ=255)
                                                                                        aaatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg             >  2:1484792/4‑100 (MQ=255)
                                                                                         aatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa            >  1:3627312/3‑100 (MQ=255)
                                                                                         aatCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa            >  1:1075869/3‑100 (MQ=255)
                                                                                          ttCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc           >  2:482977/1‑100 (MQ=255)
                                                                                          atCAGGCGGACTGGCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc           >  1:3427079/2‑100 (MQ=255)
                                                                                           tCAGTCGGCCAGTCCGCCTGATTTCATTTTATGGATATTCATGATGAATAAATCGCTTGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa          <  1:703561/100‑1 (MQ=255)
                                                                                            cAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAg         <  1:3168173/100‑1 (MQ=255)
                                                                                            cAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAg         >  1:2529539/1‑100 (MQ=255)
                                                                                              ggcggcCAGTCCGGCTGATTTCATTTTATGGATAATCATTATAAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc       >  2:2982380/1‑100 (MQ=255)
                                                                                              ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggc       >  1:2495842/1‑100 (MQ=255)
                                                                                               gcggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcg      >  2:731989/1‑100 (MQ=255)
                                                                                                cggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcgg     <  2:3037388/100‑1 (MQ=255)
                                                                                                 ggcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcggc    >  1:57312/1‑100 (MQ=255)
                                                                                                  gcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg   >  1:4079498/1‑100 (MQ=255)
                                                                                                  gcCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg   <  2:765104/100‑1 (MQ=255)
                                                                                                   ccAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGgcgc  <  2:196203/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGC  >  NC_000913/1689718‑1689916

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: