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,4790AG99.7% 1189.8 / ‑5.0 351D198G (GAT→GGT) yibOputative 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (158/192);  total (159/192)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 365 total aligned reads displayed.

AATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATG  >  NC_000913/3783389‑3783568
                                                                                          |                                                                                         
   tCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTcgcg                                                                                      <  1:7579282/93‑1 (MQ=255)
    cGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGt                                                                                     >  2:6671668/1‑93 (MQ=255)
      aagaaAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAg                                                                                   >  1:4561096/1‑93 (MQ=255)
        gaaAAATTTGCCGCGCTGGGCAAATGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGaa                                                                                 <  1:5273461/93‑1 (MQ=255)
         aaaaaGTTGCCGCGCTGGGCAAACGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGaaa                                                                                <  1:8316383/93‑1 (MQ=255)
             aTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGc                                                                            >  2:2481212/1‑93 (MQ=255)
             aTTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGc                                                                            <  1:7179952/93‑1 (MQ=255)
              tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCt                                                                           <  2:8227302/93‑1 (MQ=255)
              tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCt                                                                           <  2:8108681/93‑1 (MQ=255)
              tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCt                                                                           <  1:5829499/93‑1 (MQ=255)
              tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCt                                                                           <  1:7735976/93‑1 (MQ=255)
              tttGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCt                                                                           <  2:379732/93‑1 (MQ=255)
                tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTa                                                                         <  2:3519261/93‑1 (MQ=255)
                tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTa                                                                         <  1:6721797/93‑1 (MQ=255)
                tGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTa                                                                         <  1:4065441/93‑1 (MQ=255)
                  ccGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATg                                                                       >  1:2945862/1‑93 (MQ=255)
                   cgcgCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGa                                                                      >  1:2033212/1‑93 (MQ=255)
                    gcgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATATCCGTTGGGGTCGCGTAGAAAAAGCTTATGAc                                                                     >  2:4858753/1‑93 (MQ=255)
                    gcgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTTTGAc                                                                     <  2:5603534/93‑1 (MQ=255)
                    gcgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGAc                                                                     >  2:5918026/1‑93 (MQ=255)
                    gcgcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGAc                                                                     >  2:2462192/1‑93 (MQ=255)
                      gcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGAccg                                                                   >  1:489396/1‑92 (MQ=255)
                      gcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACct                                                                   >  1:5541520/1‑93 (MQ=255)
                      gcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACct                                                                   <  1:3723118/93‑1 (MQ=255)
                      gcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACct                                                                   >  2:2646006/1‑93 (MQ=255)
                      gcTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACct                                                                   >  2:4410239/1‑93 (MQ=255)
                       cTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACctg                                                                  >  2:4710370/1‑93 (MQ=255)
                          ggCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACctgctg                                                               <  1:133555/93‑1 (MQ=255)
                           gCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGTGGTCGCGTAGAAAAAGCTTATGACCTGCTGa                                                              >  1:5054638/1‑93 (MQ=255)
                           gCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGa                                                              <  2:7363904/93‑1 (MQ=255)
                            cAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGAc                                                             <  1:8705093/93‑1 (MQ=255)
                             aaaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGAct                                                            >  1:3982070/1‑93 (MQ=255)
                              aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActc                                                           >  2:2036604/1‑93 (MQ=255)
                              aaGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActc                                                           >  2:117999/1‑93 (MQ=255)
                               agtccGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATATCCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActct                                                          <  1:4687602/90‑1 (MQ=255)
                               aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActct                                                          <  1:4373937/93‑1 (MQ=255)
                               aGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGActct                                                          <  2:1521529/93‑1 (MQ=255)
                                ggCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTg                                                         >  2:6184178/1‑93 (MQ=255)
                                 gCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTgg                                                        <  1:6888371/93‑1 (MQ=255)
                                       tAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAgg                                                  >  2:3097470/1‑93 (MQ=255)
                                       tAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAgg                                                  <  2:4005441/93‑1 (MQ=255)
                                       tAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAgg                                                  <  2:5510816/93‑1 (MQ=255)
                                         gCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGc                                                >  2:1881591/1‑93 (MQ=255)
                                           gTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGa                                              >  1:175420/1‑93 (MQ=255)
                                            tCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAg                                             <  1:1936551/93‑1 (MQ=255)
                                             ccATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGt                                            <  2:6599716/93‑1 (MQ=255)
                                              catcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGtt                                           <  1:3066349/93‑1 (MQ=255)
                                              catcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGtt                                           <  2:5169617/93‑1 (MQ=255)
                                              catcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGACGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGtt                                           >  2:5697568/1‑93 (MQ=255)
                                               atcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTc                                          <  1:3179041/93‑1 (MQ=255)
                                               atcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTc                                          >  2:7117338/1‑93 (MQ=255)
                                               atcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTc                                          <  1:5180585/93‑1 (MQ=255)
                                                tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTcc                                         <  1:8447024/93‑1 (MQ=255)
                                                tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTcc                                         <  1:1160270/93‑1 (MQ=255)
                                                tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTcc                                         >  1:253330/1‑93 (MQ=255)
                                                tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTcc                                         <  2:2969506/93‑1 (MQ=255)
                                                tcatTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTcc                                         <  1:6503318/93‑1 (MQ=255)
                                                 catTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCa                                        <  2:1638442/93‑1 (MQ=255)
                                                 catTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCa                                        <  2:5287914/93‑1 (MQ=255)
                                                   ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGANTNNGGCGCAGGGCGGNTTCAAgn                                      >  1:2432481/1‑92 (MQ=255)
                                                   ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAgg                                      >  1:100272/1‑93 (MQ=255)
                                                   ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAgg                                      >  2:239034/1‑93 (MQ=255)
                                                   ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAgg                                      <  2:6257560/93‑1 (MQ=255)
                                                   ttGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAgc                                      >  1:76891/1‑92 (MQ=255)
                                                    tGGTCGTGACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGc                                     <  1:5090111/93‑1 (MQ=255)
                                                         gTTACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATa                                <  2:412255/93‑1 (MQ=255)
                                                          ttACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAc                               <  1:621165/93‑1 (MQ=255)
                                                          ttACTACGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAc                               <  1:904936/93‑1 (MQ=255)
                                                               acGCGATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATAccgccg                          <  2:6823871/93‑1 (MQ=255)
                                                                cgcgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGt                         <  2:357919/93‑1 (MQ=255)
                                                                 gcgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGtt                        >  1:485182/1‑93 (MQ=255)
                                                                 gcgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGtt                        <  1:7327092/93‑1 (MQ=255)
                                                                 gcgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGtt                        <  1:2176175/93‑1 (MQ=255)
                                                                 gcgAGGGACCGCGATAGCCGTGGGGGTCGCGTAGAAAAATCTTATGACCTGTTGGCTCTGGCGCAGGGCGATTTCCAGGCCGATACCGCCGtt                        <  2:8606672/93‑1 (MQ=255)
                                                                  cgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTg                       <  1:2692803/93‑1 (MQ=255)
                                                                  cgATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTg                       <  2:1805416/93‑1 (MQ=255)
                                                                   gATGGACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGc                      >  2:3918933/1‑93 (MQ=255)
                                                                      ggACCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTgg                   <  2:3213404/93‑1 (MQ=255)
                                                                        aCCGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGtt                 <  2:1303983/93‑1 (MQ=255)
                                                                         ccGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                >  2:6724368/1‑93 (MQ=255)
                                                                         ccGCGATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGttt                <  1:6743994/93‑1 (MQ=255)
                                                                            cgATAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCa             <  1:3305249/93‑1 (MQ=255)
                                                                               tAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgc          >  1:944147/1‑93 (MQ=255)
                                                                               tAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgc          >  1:8122569/1‑93 (MQ=255)
                                                                               tAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgc          >  2:6911991/1‑93 (MQ=255)
                                                                               tAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgc          <  2:2759847/93‑1 (MQ=255)
                                                                               tAACCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgc          <  2:42844/93‑1 (MQ=255)
                                                                                aaCCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct         <  2:754926/93‑1 (MQ=255)
                                                                                aaCCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgct         <  1:7565356/93‑1 (MQ=255)
                                                                                 aCCGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctg        >  2:7664896/1‑93 (MQ=255)
                                                                                  ccGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGGGCAGGGCGAGTTCCATGCCGATACAGCCGTTGCTGGTTTGCAGcctgc       >  2:6518684/1‑93 (MQ=255)
                                                                                  ccGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgc       >  2:2586546/1‑93 (MQ=255)
                                                                                   cGTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGCGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGgctgtt      >  1:2954738/1‑91 (MQ=255)
                                                                                    gTTGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCtt     >  2:3439319/1‑93 (MQ=255)
                                                                                     ttGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa    <  1:6791281/93‑1 (MQ=255)
                                                                                     ttGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa    <  1:5634863/93‑1 (MQ=255)
                                                                                     ttGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa    <  2:175303/93‑1 (MQ=255)
                                                                                     ttGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTa    >  2:2222141/1‑93 (MQ=255)
                                                                                      tGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt   <  1:7129788/93‑1 (MQ=255)
                                                                                      tGGGGTCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTAt   >  1:1104366/1‑93 (MQ=255)
                                                                                          |                                                                                         
AATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATG  >  NC_000913/3783389‑3783568

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