Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,345 A→G 6.9% intergenic (‑141/+167) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,177,3450AG 6.9% 271.9 / 6.8 84intergenic (‑141/+167)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base A (46/32);  new base G (3/3);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGT  >  NC_000913/2177197‑2177492
                                                                                                                                                     |                                                                                                                                                   
aCGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTCGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGaa                                                                                                                                                     <  4:590096/149‑1 (MQ=255)
 cGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAg                                                                                                                                                    <  4:582725/149‑1 (MQ=255)
  gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGc                                                                                                                                                   >  2:1113262/1‑149 (MQ=255)
  gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGc                                                                                                                                                   <  3:63232/149‑1 (MQ=255)
  gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGc                                                                                                                                                   >  2:695267/1‑149 (MQ=255)
   tACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGcc                                                                                                                                                  <  1:1113262/149‑1 (MQ=255)
    aCATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCg                                                                                                                                                 >  2:970896/1‑149 (MQ=255)
      aTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATAGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                                                                               <  4:190589/149‑1 (MQ=255)
       ttttCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                                                                                                                              <  4:1124422/149‑1 (MQ=255)
             tatCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGt                                                                                                                                        <  1:357037/149‑1 (MQ=255)
                      gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                                                                                                                               <  1:507141/149‑1 (MQ=255)
                          gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAATTTCGCGTGTTTTTa                                                                                                                           <  3:85260/149‑1 (MQ=255)
                           ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt                                                                                                                          >  2:964939/1‑149 (MQ=255)
                             ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                                                                                                                        <  2:340231/149‑1 (MQ=255)
                                    tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc                                                                                                                 >  3:638579/1‑149 (MQ=255)
                                             aaTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                                                                        >  2:504654/1‑149 (MQ=255)
                                               tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACAt                                                                                                      >  3:874079/1‑149 (MQ=255)
                                                gAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa                                                                                                     >  2:325600/1‑149 (MQ=255)
                                                 aaaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                                    <  1:695267/149‑1 (MQ=255)
                                                 aaaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                                    >  2:111001/1‑149 (MQ=255)
                                                  aaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACg                                                                                                   >  4:922274/1‑149 (MQ=255)
                                                   attattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGt                                                                                                  <  3:922274/149‑1 (MQ=255)
                                                    ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                                                                               <  3:948202/103‑1 (MQ=255)
                                                    ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                                                                               >  4:948202/1‑103 (MQ=255)
                                                    ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                                                                                 >  3:409295/1‑149 (MQ=255)
                                                     tattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTgg                                                                                                >  4:1125029/1‑149 (MQ=255)
                                                       tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                                                                                              <  4:874079/149‑1 (MQ=255)
                                                       tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                                                                                              <  3:1012296/149‑1 (MQ=255)
                                                       tttGGTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                                                                                              <  3:1125029/149‑1 (MQ=255)
                                                         tGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGCGGCAGt                                                                                            >  2:238395/1‑149 (MQ=255)
                                                            ttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                                                                         <  2:1145573/149‑1 (MQ=255)
                                                                           tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATaa                                                                          >  1:1013647/1‑149 (MQ=255)
                                                                            aaCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAg                                                                         <  1:970896/149‑1 (MQ=255)
                                                                                   acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCtt                                                                   >  4:1146738/1‑148 (MQ=255)
                                                                                   acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCtt                                                                   <  3:1146738/148‑1 (MQ=255)
                                                                                     aaaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCa                                                                 >  1:899862/1‑149 (MQ=255)
                                                                                     aaaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCtgaag                                                                >  2:734833/1‑145 (MQ=255)
                                                                                      aaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGc                                                               >  1:336995/1‑149 (MQ=255)
                                                                                            aCAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCa                                                         <  1:964939/149‑1 (MQ=255)
                                                                                            aCAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTGAAGCGTctta                                                         <  1:734833/149‑4 (MQ=255)
                                                                                             cAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCAt                                                        >  4:299314/1‑149 (MQ=255)
                                                                                              aaTATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATc                                                       >  2:24326/1‑149 (MQ=255)
                                                                                                tatGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                                                     <  4:409295/149‑1 (MQ=255)
                                                                                                 atGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                                    <  3:918326/101‑1 (MQ=255)
                                                                                                 atGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                                    >  4:918326/1‑101 (MQ=255)
                                                                                                 atGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGCTTTTATGCCGGATGCGGTGTAAACATACGTGGCACTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAgg                                                    >  4:938114/1‑149 (MQ=255)
                                                                                                    aaaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGTCATGATAAGACGCTTCAGCGTCGCATCAGGCAt                                                 >  1:903422/1‑149 (MQ=255)
                                                                                                    aaaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGa                                                                        <  4:306643/126‑1 (MQ=255)
                                                                                                    aaaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGa                                                                        >  3:306643/1‑126 (MQ=255)
                                                                                                    aaaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCAt                                                 >  1:1068838/1‑149 (MQ=255)
                                                                                                     aaaTTATTCGAGGAGTGAAAGGCAAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCAt                                                 <  2:899862/149‑1 (MQ=255)
                                                                                                     aaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATa                                                >  2:1094412/1‑149 (MQ=255)
                                                                                                           ttCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGc                                          <  2:1068838/149‑1 (MQ=255)
                                                                                                              gaggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                                       >  3:817613/1‑149 (MQ=255)
                                                                                                              gaggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                                       >  1:68304/1‑149 (MQ=255)
                                                                                                               aggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGa                                      >  3:974308/1‑149 (MQ=255)
                                                                                                                    tGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGGGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAg                                 >  2:982065/1‑149 (MQ=255)
                                                                                                                     gAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGc                                >  2:875334/1‑149 (MQ=255)
                                                                                                                       aaGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGt                              >  4:534236/1‑149 (MQ=255)
                                                                                                                       aaGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGt                              >  1:917727/1‑149 (MQ=255)
                                                                                                                        aGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                                                     >  1:1084391/1‑125 (MQ=255)
                                                                                                                        aGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                                                     <  2:1084391/125‑1 (MQ=255)
                                                                                                                          gC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTaaa                           >  2:583418/1‑149 (MQ=255)
                                                                                                                           c‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc                          <  2:903422/149‑1 (MQ=255)
                                                                                                                           c‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGCGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc                          <  1:238395/149‑1 (MQ=255)
                                                                                                                             aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg                         >  4:1037349/1‑149 (MQ=255)
                                                                                                                             aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg                         >  3:775196/1‑149 (MQ=255)
                                                                                                                             aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg                         >  4:235748/1‑149 (MQ=255)
                                                                                                                                aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCGGGCATAGGTTGCCGGATGCAGCGTAAACGCCt                      >  4:947729/1‑149 (MQ=255)
                                                                                                                                aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCt                      >  2:199098/1‑149 (MQ=255)
                                                                                                                                 aaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCtt                     <  2:68304/149‑1 (MQ=255)
                                                                                                                                  taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGc                        <  4:1139069/144‑1 (MQ=255)
                                                                                                                                  taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGc                        >  3:1139069/2‑145 (MQ=39)
                                                                                                                                   aCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTAt                   <  1:504654/149‑1 (MQ=255)
                                                                                                                                     ggcttCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTaaa                           <  4:875519/135‑1 (MQ=39)
                                                                                                                                     ggCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTaaa                           >  3:875519/1‑139 (MQ=255)
                                                                                                                                        cTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTc              >  1:836384/1‑149 (MQ=255)
                                                                                                                                         tCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTcc             >  2:647344/1‑149 (MQ=255)
                                                                                                                                         tCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTcc             >  3:95539/1‑149 (MQ=255)
                                                                                                                                         tCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTcc             <  1:325600/149‑1 (MQ=255)
                                                                                                                                              aTAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACgg        <  4:638579/149‑1 (MQ=255)
                                                                                                                                                  ggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCTGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtg    <  3:938114/149‑1 (MQ=255)
                                                                                                                                                  ggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtg    >  1:922532/1‑149 (MQ=255)
                                                                                                                                                  ggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtg    >  3:230279/1‑149 (MQ=255)
                                                                                                                                                   gAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtgg   >  2:755944/1‑149 (MQ=255)
                                                                                                                                                    aaGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtggt  <  3:299314/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGT  >  NC_000913/2177197‑2177492

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