Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,783,479 A→G D198G (GAT→GGT)  yibO → putative 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,783,4790AG86.1% 668.5 / 99.8 288D198G (GAT→GGT) yibOputative 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (18/22);  new base G (113/134);  total (132/156)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 292 total aligned reads displayed.

AAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATG  >  NC_000913/3783387‑3783568
                                                                                            |                                                                                         
  aaTTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTc                                                                                         <  1:6099243/93‑1 (MQ=255)
  aaTTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTc                                                                                         <  2:5750936/93‑1 (MQ=255)
  aaTTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTc                                                                                         >  2:2508579/1‑93 (MQ=255)
    ttCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGTGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATcgc                                                                                       <  1:1877468/93‑1 (MQ=255)
     tCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTcgcg                                                                                      >  1:358596/1‑93 (MQ=255)
       gaagaaAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTa                                                                                    <  2:1567952/93‑1 (MQ=255)
       gaagaaAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTa                                                                                    >  2:1663425/1‑93 (MQ=255)
         agaaAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGa                                                                                  >  2:4422748/1‑93 (MQ=255)
           aaaaaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGaaa                                                                                >  1:5588323/1‑93 (MQ=255)
             aaaTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGaaaaa                                                                              <  2:6100094/93‑1 (MQ=255)
               aTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGc                                                                            <  1:3380306/93‑1 (MQ=255)
                tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCt                                                                           <  2:5581990/93‑1 (MQ=255)
                  tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTa                                                                         <  1:4115735/93‑1 (MQ=255)
                  tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTa                                                                         <  2:4663174/93‑1 (MQ=255)
                  tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTa                                                                         <  1:3441194/93‑1 (MQ=255)
                   gCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTAt                                                                        <  1:3993346/93‑1 (MQ=255)
                   gCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTAt                                                                        <  2:5138668/93‑1 (MQ=255)
                      gcgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGAc                                                                     >  2:703057/1‑93 (MQ=255)
                      gcgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGAc                                                                     <  1:4320922/93‑1 (MQ=255)
                      gcgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTACAAAAAGCTTATGAc                                                                     <  2:6445265/93‑1 (MQ=255)
                       cgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGTGGGCGCGTAGAAAAAGCTTATGAcc                                                                    >  2:4258619/1‑93 (MQ=255)
                        gcTGGGCAAATGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACct                                                                   >  1:4852748/1‑93 (MQ=255)
                        gcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACct                                                                   <  1:2022330/93‑1 (MQ=255)
                        gcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACct                                                                   <  2:1850543/93‑1 (MQ=255)
                        gcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACct                                                                   <  1:5756742/93‑1 (MQ=255)
                         cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACctg                                                                  >  2:1438089/1‑93 (MQ=255)
                         cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACctg                                                                  >  2:2240293/1‑93 (MQ=255)
                          tGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACctgc                                                                 <  2:6764638/93‑1 (MQ=255)
                          tGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACctgc                                                                 <  2:359750/93‑1 (MQ=255)
                           gggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACctgct                                                                >  2:3939629/1‑93 (MQ=255)
                            ggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACctgctg                                                               >  2:4554718/1‑93 (MQ=255)
                             gCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGa                                                              <  2:4203740/93‑1 (MQ=255)
                              cAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGAc                                                             <  2:901580/93‑1 (MQ=255)
                                aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActc                                                           >  1:6200191/1‑93 (MQ=255)
                                aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActc                                                           <  1:5407131/93‑1 (MQ=255)
                                aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActc                                                           >  2:2974453/1‑93 (MQ=255)
                                aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActc                                                           <  1:4952375/93‑1 (MQ=255)
                                 aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActct                                                          >  1:3951758/1‑93 (MQ=255)
                                 aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGActct                                                          <  2:7112692/93‑1 (MQ=255)
                                  ggCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTg                                                         <  1:5964236/93‑1 (MQ=255)
                                   gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTgg                                                        <  2:6950147/93‑1 (MQ=255)
                                   gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTgg                                                        <  2:3751891/93‑1 (MQ=255)
                                     cgcgTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGgcg                                                      <  1:691702/93‑1 (MQ=255)
                                       cgTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCa                                                    <  2:730522/93‑1 (MQ=255)
                                        gTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAg                                                   <  1:1748397/93‑1 (MQ=255)
                                        gTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAg                                                   <  1:3076137/93‑1 (MQ=255)
                                          aGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAggg                                                 <  1:6005185/93‑1 (MQ=255)
                                          aGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAggg                                                 <  1:3695516/93‑1 (MQ=255)
                                           gCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGc                                                <  1:5883131/93‑1 (MQ=255)
                                           gCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGc                                                >  2:3210344/1‑93 (MQ=255)
                                           gCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGc                                                >  1:2183058/1‑93 (MQ=255)
                                             gTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGa                                              <  2:2009247/93‑1 (MQ=255)
                                              tCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAg                                             <  2:2842541/93‑1 (MQ=255)
                                                 atcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTc                                          <  2:6214966/93‑1 (MQ=255)
                                                 atcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTc                                          >  1:2738123/1‑93 (MQ=255)
                                                  tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTcc                                         >  1:2618918/1‑93 (MQ=255)
                                                  tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTcc                                         <  1:4423878/93‑1 (MQ=255)
                                                  tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTcc                                         >  1:5131250/1‑93 (MQ=255)
                                                    atTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAg                                       >  2:5888523/1‑93 (MQ=255)
                                                     ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAgg                                      <  2:6473932/93‑1 (MQ=255)
                                                     ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAgg                                      >  1:7236029/1‑93 (MQ=255)
                                                      tGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGc                                     <  1:7352560/93‑1 (MQ=255)
                                                        gTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCg                                   >  1:141170/1‑93 (MQ=255)
                                                         tCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGa                                  <  1:5691751/93‑1 (MQ=255)
                                                            ttACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAc                               <  2:3057435/93‑1 (MQ=255)
                                                               ctacGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAccgc                            <  2:2350495/93‑1 (MQ=255)
                                                                tacGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAccgcc                           >  2:2076538/1‑93 (MQ=255)
                                                                tacGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAccgcc                           <  1:6989782/93‑1 (MQ=255)
                                                                tacGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAccgcc                           >  2:4030470/1‑93 (MQ=255)
                                                                tacGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAccgcc                           >  1:1291182/1‑93 (MQ=255)
                                                                  cgcgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGt                         <  2:1089708/93‑1 (MQ=255)
                                                                  cgcgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGt                         <  1:1648815/93‑1 (MQ=255)
                                                                   gcgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGtt                        >  2:6559926/1‑93 (MQ=255)
                                                                   gcgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGtt                        >  1:5827841/1‑93 (MQ=255)
                                                                   gcgATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGtt                        <  2:1191529/93‑1 (MQ=255)
                                                                     gATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGc                      >  1:802179/1‑93 (MQ=255)
                                                                     gATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGc                      <  1:6717652/93‑1 (MQ=255)
                                                                      aTGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                     >  2:764304/1‑93 (MQ=255)
                                                                      aTGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                     <  1:4556197/93‑1 (MQ=255)
                                                                      aTGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCt                     <  1:3476865/93‑1 (MQ=255)
                                                                       tGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTg                    <  2:506383/93‑1 (MQ=255)
                                                                       tGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTg                    <  2:4348881/93‑1 (MQ=255)
                                                                        ggACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTgg                   >  1:1074334/1‑93 (MQ=255)
                                                                           ccGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                <  2:537184/93‑1 (MQ=255)
                                                                            cgcgATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTg               <  1:2086072/93‑1 (MQ=255)
                                                                              cgATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGGCCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCa             <  1:5508539/93‑1 (MQ=255)
                                                                              cgATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCa             <  1:2923215/93‑1 (MQ=255)
                                                                               gATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAg            <  1:6629381/93‑1 (MQ=255)
                                                                               gATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAg            >  2:2691039/1‑93 (MQ=255)
                                                                                 tAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgc          >  2:1485750/1‑93 (MQ=255)
                                                                                  aaCCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct         <  1:363722/93‑1 (MQ=255)
                                                                                    ccGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgc       >  2:6900852/1‑93 (MQ=255)
                                                                                     cGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgct      <  1:6489223/93‑1 (MQ=255)
                                                                                       ttGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa    <  2:1759804/93‑1 (MQ=255)
                                                                                       ttGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa    <  1:2404392/93‑1 (MQ=255)
                                                                                       ttGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa    <  1:6369451/93‑1 (MQ=255)
                                                                                       ttGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa    >  1:7458581/1‑93 (MQ=255)
                                                                                        tGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt   <  1:3193355/93‑1 (MQ=255)
                                                                                        tGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt   <  1:7200723/93‑1 (MQ=255)
                                                                                         ggggTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAACGCCGTTGCTGGTTTGCAGGCTGCTTATg  >  1:1155980/1‑93 (MQ=255)
                                                                                            |                                                                                         
AAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATG  >  NC_000913/3783387‑3783568

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