Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,342,683 T→C G137G (GGT→GGC)  kdsC → 3‑deoxy‑D‑manno‑octulosonate 8‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,342,6830TC84.2% 199.5 / 23.7 76G137G (GGT→GGC) kdsC3‑deoxy‑D‑manno‑octulosonate 8‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (6/6);  new base C (34/30);  total (40/36)
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.57e-01

ACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAACTGGATGAAG  >  NC_000913/3342549‑3342819
                                                                                                                                      |                                                                                                                                        
actcactTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTaa                                                                                                                                      >  2:706215/1‑139 (MQ=255)
   cactTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCg                                                                                                                                   <  1:158939/139‑1 (MQ=255)
        gTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCg                                                                                                                                   >  1:825950/1‑134 (MQ=255)
        gTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCGGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCg                                                                                                                                   <  2:825950/134‑1 (MQ=255)
         tATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCg                                                                                                                                   >  2:1608976/1‑133 (MQ=255)
         tATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCg                                                                                                                                   <  1:1608976/133‑1 (MQ=255)
           tCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTg                                                                                                                           >  2:906002/1‑139 (MQ=255)
             aGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGc                                                                                                                         <  2:22505/139‑1 (MQ=255)
              ggggCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGcc                                                                                                                        >  2:67561/1‑139 (MQ=255)
              ggggCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGcc                                                                                                                        >  1:817568/1‑139 (MQ=255)
                ggCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGa                                                                                                                      >  2:36858/1‑139 (MQ=255)
                   aGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATgc                                                                                                                   >  2:106422/1‑139 (MQ=255)
                     tCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATgcgc                                                                                                                 <  1:106422/139‑1 (MQ=255)
                      cagacaaacTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCa                                                                                                                <  2:446837/136‑1 (MQ=255)
                        aacaaacTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATc                                                                                                              <  2:1374996/139‑1 (MQ=255)
                         acaaacTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATcc                                                                                                             >  1:1024503/1‑139 (MQ=255)
                                 aTCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGa                                                                                                     >  2:1476089/1‑139 (MQ=255)
                                   cGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATc                                                                                                   >  1:464130/1‑139 (MQ=255)
                                   cGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATc                                                                                                   >  1:182708/1‑139 (MQ=255)
                                     ccTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATccc                                                                                                 >  2:308502/1‑139 (MQ=255)
                                       tttAGCGCTCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGCCGGCGACGATCCCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCcgc                                                                                               >  1:1453272/1‑139 (MQ=255)
                                        ttAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGAATGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCcgcg                                                                                              <  2:1024503/139‑1 (MQ=255)
                                           gCGATCTGCTGGAAAAACTGGCGATTGCCCGGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCgcgcgc                                                                                           <  1:36858/139‑1 (MQ=255)
                                           gCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCAATGCGCATCCACTGTTGATCCCgcgcgc                                                                                           >  1:1368269/1‑139 (MQ=255)
                                                ctgctgGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATTCACTGTTg                                                                                                      <  2:1070652/123‑1 (MQ=255)
                                                ctgctgGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTg                                                                                                      >  1:1070652/1‑123 (MQ=255)
                                                ctgctgGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGAtt                                                                                      >  2:17367/1‑139 (MQ=255)
                                                ctgctgGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGAtt                                                                                      >  2:1140937/1‑139 (MQ=255)
                                                ctgctgGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGAtt                                                                                      >  2:626254/1‑139 (MQ=255)
                                                    tgGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGt                                                                                  >  1:58389/1‑139 (MQ=255)
                                                             tGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCATTGTTGATCCCGCGCGCCGATTACGTGACGCGCAt                                                                         >  1:688241/1‑139 (MQ=255)
                                                             tGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGa                                                                                                                      >  2:1531629/1‑94 (MQ=255)
                                                             tGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGa                                                                                                                      <  1:1531629/94‑1 (MQ=255)
                                                             tGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCAt                                                                                                               <  1:449133/101‑1 (MQ=255)
                                                             tGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCAt                                                                                                               >  2:449133/1‑101 (MQ=255)
                                                                    tGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGc                                                                  <  1:557335/139‑1 (MQ=255)
                                                                        ccGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCATTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTc                                                              <  2:688241/139‑1 (MQ=255)
                                                                               aTGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGgcgc                                                       <  2:464130/139‑1 (MQ=255)
                                                                               aTGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATACCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGgcgc                                                       <  1:1162833/139‑1 (MQ=255)
                                                                               aTGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCAATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGgcgc                                                       <  2:1368269/139‑1 (MQ=255)
                                                                               aTGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGAGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGtcgc                                                       <  1:1476089/139‑1 (MQ=255)
                                                                               aTGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGAGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGgcgc                                                       <  1:645625/139‑1 (MQ=255)
                                                                                        aTGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGa                                                                                                                      >  2:735840/1‑67 (MQ=255)
                                                                                        aTGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGa                                                                                                                      <  1:735840/67‑1 (MQ=255)
                                                                                          gTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAg                                            >  2:1461265/1‑139 (MQ=255)
                                                                                              gCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTg                                        <  1:1159212/139‑1 (MQ=255)
                                                                                              gCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTg                                        >  2:1346247/1‑139 (MQ=255)
                                                                                                 atgatCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGa                                     >  1:345305/1‑139 (MQ=255)
                                                                                                  tgatCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCAt                                                                                                               >  2:1033899/1‑64 (MQ=255)
                                                                                                  tgatCTCATCGACTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCAt                                                                                                               <  1:1033899/64‑1 (MQ=255)
                                                                                                   gatCTCATCGACTGGCCGGTAATGG‑AAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGGTCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACtt                                  >  1:946806/1‑139 (MQ=255)
                                                                                                              cTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATgcgc                                                                                                                 >  1:1332019/1‑50 (MQ=255)
                                                                                                              cTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATgcgc                                                                                                                 <  2:1332019/50‑1 (MQ=255)
                                                                                                              cTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTg                                                                                 >  1:1448721/1‑82 (MQ=255)
                                                                                                              cTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTg                                                                                 <  2:1448721/82‑1 (MQ=255)
                                                                                                              cTGGCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTg                        >  2:1177228/1‑139 (MQ=255)
                                                                                                               tGGCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCg                                                                                                                             >  2:931417/1‑37 (MQ=255)
                                                                                                               tGGCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCg                                                                                                                             <  1:931417/37‑1 (MQ=255)
                                                                                                                 gCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCgg                                                                <  2:343002/96‑1 (MQ=255)
                                                                                                                 gCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCgg                                                                >  1:343002/1‑96 (MQ=255)
                                                                                                                 gCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATcc                                                                                                             <  2:1381568/51‑1 (MQ=255)
                                                                                                                 gCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATcc                                                                                                             >  1:1381568/1‑51 (MQ=255)
                                                                                                                 gCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGgcg                     <  2:345305/139‑1 (MQ=255)
                                                                                                                 gCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGgcg                     >  1:268763/1‑139 (MQ=255)
                                                                                                                 gCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGgcg                     <  2:1453272/139‑1 (MQ=255)
                                                                                                                 gCCGGTAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGgcg                     <  1:67561/139‑1 (MQ=255)
                                                                                                                      tAATGGAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAggg                <  2:817568/139‑1 (MQ=255)
                                                                                                                           gAAAAAGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAAc           <  1:308502/139‑1 (MQ=255)
                                                                                                                            aaaaaGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGa                                     >  1:1024148/1‑112 (MQ=255)
                                                                                                                            aaaaaGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTCTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGa                                     <  2:1024148/112‑1 (MQ=255)
                                                                                                                               aaGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAACTGGa       <  2:946806/139‑1 (MQ=255)
                                                                                                                                aGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGttt                                         >  1:365151/1‑104 (MQ=255)
                                                                                                                                aGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGttt                                         <  2:365151/104‑1 (MQ=255)
                                                                                                                                aGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGa                                     >  1:847329/1‑108 (MQ=255)
                                                                                                                                aGTGGGCTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGa                                     <  2:847329/108‑1 (MQ=255)
                                                                                                                                    ggTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAACTGGATGAAg  >  1:985407/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                        
ACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTATGTCGGCGATGATCTCATCGACTGGCCGGTAATGGAAAAAGTGGGTTTAAGCGTCGCCGTGGCCGATGCGCATCCACTGTTGATCCCGCGCGCCGATTACGTGACGCGCATTGCTGGCGGTCGTGGCGCAGTGCGCGAAGTTTGCGACTTATTACTCCTGGCGCAGGGCAAACTGGATGAAG  >  NC_000913/3342549‑3342819

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