Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 649,874 G→T S206Y (TCC→TAC)  citE ← citrate lyase beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913649,8740GT94.5% 242.1 / 8.1 73S206Y (TCC→TAC) citEcitrate lyase beta subunit
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (33/36);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

ACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACC  >  NC_000913/649739‑650004
                                                                                                                                       |                                                                                                                                  
aCTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGCTGTTAGCGTCGTAATa                                                                                                                                >  2:497982/1‑140 (MQ=255)
 cttcttTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAg                                                                                                                               >  1:10318/1‑140 (MQ=255)
  ttctttctGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGa                                                                                                                              <  1:446650/140‑1 (MQ=255)
  ttctttctGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGa                                                                                                                              <  2:1075047/140‑1 (MQ=255)
   tctttctGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGAc                                                                                                                             >  2:830950/1‑140 (MQ=255)
       tctGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTa                                                                                                                         >  2:64749/1‑140 (MQ=255)
            gTCGGTGCGTAGAGGTTGTGTAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGaa                                                                                                                    >  1:591742/1‑140 (MQ=255)
            gTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGaa                                                                                                                    >  2:363339/1‑140 (MQ=255)
              cGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAAc                                                                                                                   >  1:1341005/1‑139 (MQ=255)
              cGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAAc                                                                                                                   <  2:1341005/139‑1 (MQ=255)
               ggTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGc                                                                                                                 >  1:76491/1‑140 (MQ=255)
                             gtgCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCag                                                                                                   >  1:1167191/1‑140 (MQ=255)
                             gtgCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCag                                                                                                   >  1:763620/1‑140 (MQ=255)
                             gtgCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCag                                                                                                   >  2:202153/1‑140 (MQ=255)
                              tgCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCaga                                                                                                  >  2:641656/1‑140 (MQ=255)
                              tgCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCaga                                                                                                  >  2:750081/1‑140 (MQ=255)
                                  gcagATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTAGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAgcgc                                                                                              >  2:64409/1‑140 (MQ=255)
                                   cagATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAgcgcg                                                                                             <  2:76491/140‑1 (MQ=255)
                                      aTCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCgg                                                                                          >  2:906563/1‑140 (MQ=255)
                                       tCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGc                                                                                         <  2:1325046/140‑1 (MQ=255)
                                       tCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGc                                                                                         <  1:647336/140‑1 (MQ=255)
                                       tCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGc                                                                                         <  1:351079/140‑1 (MQ=255)
                                       tCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGc                                                                                         <  2:511724/140‑1 (MQ=255)
                                          aTCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTg                                                                                      <  1:830950/140‑1 (MQ=255)
                                          aTCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTg                                                                                      <  2:1167191/140‑1 (MQ=255)
                                          aTCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTg                                                                                      <  2:763620/140‑1 (MQ=255)
                                            cTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCa                                                                                    <  2:27556/140‑1 (MQ=255)
                                             tGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCaa                                                                                   >  2:806126/1‑140 (MQ=255)
                                             tGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCa                                                                                    >  2:1365285/1‑139 (MQ=255)
                                             tGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCa                                                                                    <  1:1365285/139‑1 (MQ=255)
                                               aCGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCaaaa                                                                                 >  1:1198499/1‑140 (MQ=255)
                                                cgcgGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCaaaaa                                                                                >  1:689473/1‑139 (MQ=255)
                                                    ggTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGaa                                                                            >  1:779306/1‑140 (MQ=255)
                                                    ggTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGaa                                                                            >  2:778657/1‑140 (MQ=255)
                                                     gTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAAc                                                                           >  2:1293772/1‑140 (MQ=255)
                                                         aTCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCg                                                                       >  1:1404248/1‑140 (MQ=255)
                                                           cAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTg                                                                     <  2:482647/140‑1 (MQ=255)
                                                             gCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCg                                                                   <  1:939769/140‑1 (MQ=255)
                                                              cGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGa                                                                  >  1:894869/1‑140 (MQ=255)
                                                                  ttGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAAcag                                                              <  1:677944/140‑1 (MQ=255)
                                                                  ttGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAAcag                                                              <  1:812198/140‑1 (MQ=255)
                                                                  ttGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGCTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAAcag                                                              <  1:497982/140‑1 (MQ=255)
                                                                   tGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTg                                                                     <  2:1453104/132‑1 (MQ=255)
                                                                   tGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTg                                                                     >  1:1453104/1‑132 (MQ=255)
                                                                   tGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAAcagc                                                             >  2:973831/1‑140 (MQ=255)
                                                                       gTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGcagtt                                                         <  1:105582/140‑1 (MQ=255)
                                                                           aaGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGcagttcagt                                                     <  1:363339/140‑1 (MQ=255)
                                                                           aaGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGcagttcagt                                                     >  2:140030/1‑140 (MQ=255)
                                                                           aaGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTAGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGcagttcagt                                                     <  1:64409/140‑1 (MQ=255)
                                                                                 aGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGttccttc                                               >  2:301864/1‑140 (MQ=255)
                                                                                       ttGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGa                                         <  1:778657/140‑1 (MQ=255)
                                                                                       ttGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGa                                         <  2:591742/140‑1 (MQ=255)
                                                                                         gATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGc                                       <  1:504161/140‑1 (MQ=255)
                                                                                            gtgGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGtt                                    <  1:984891/140‑1 (MQ=255)
                                                                                            gtgGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGtt                                    <  2:1210351/140‑1 (MQ=255)
                                                                                            gtgGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGtt                                    <  1:1245862/140‑1 (MQ=255)
                                                                                            gtgGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGtt                                    <  1:203018/140‑1 (MQ=255)
                                                                                             tgGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTc                                   >  2:23365/1‑140 (MQ=255)
                                                                                                gcggcTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGt                                >  2:875460/1‑140 (MQ=255)
                                                                                                         tGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCg                                      >  2:485991/1‑125 (MQ=255)
                                                                                                         tGCAGAAATCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCg                                      <  1:485991/125‑1 (MQ=255)
                                                                                                                aTCCGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGcaca                <  1:1117892/140‑1 (MQ=255)
                                                                                                                  ccGGCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATa              <  2:1404248/140‑1 (MQ=255)
                                                                                                                   cGGCTTCGTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGc                            >  1:1023796/1‑125 (MQ=255)
                                                                                                                   cGGCTTCGTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGc                            >  1:1003864/1‑125 (MQ=255)
                                                                                                                   cGGCTTCGTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGc                            <  2:1003864/125‑1 (MQ=255)
                                                                                                                   cGGCTTCGTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGc                            <  2:1023796/125‑1 (MQ=255)
                                                                                                                    ggCTTCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGt            <  1:1028004/140‑1 (MQ=255)
                                                                                                                        tCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc        >  1:1304360/1‑140 (MQ=255)
                                                                                                                        tCGTTGTTAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc        <  2:10318/140‑1 (MQ=255)
                                                                                                                           ttgttAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc     >  1:55598/1‑140 (MQ=255)
                                                                                                                            tgttAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCa    <  1:641656/140‑1 (MQ=255)
                                                                                                                              ttAGCGTCGTAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCAcc  <  2:1220498/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
ACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGACGCGGGTTGATCAGCGATTTGCCGTCAAAGCCCAGCTGTTTGATGTGGGCGGCTTCTTGCAGAAATCCGGCTTCGTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACC  >  NC_000913/649739‑650004

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