Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,125,415 Δ1 bp 6.4% coding (45/1224 nt) glcF ← glycolate oxidase 4Fe‑4S iron‑sulfur cluster subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,125,4150C. 6.4% 264.8 / 10.8 78coding (45/1224 nt)glcFglycolate oxidase 4Fe‑4S iron‑sulfur cluster subunit
Reads supporting (aligned to +/- strand):  ref base C (47/26);  new base . (2/3);  total (49/29)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

GCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGA  >  NC_000913/3125277‑3125562
                                                                                                                                          |                                                                                                                                                   
gCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGt                                                                                                                                           >  2:201975/1‑149 (MQ=255)
  cTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTc                                                                                                                                         >  8:330887/1‑149 (MQ=255)
      cAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGc                                                                                                                                      <  7:330887/148‑1 (MQ=255)
             tGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCAtctc                                                                                                                              >  7:371184/1‑149 (MQ=255)
               tttAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCtt                                                                                                                            <  7:68787/149‑1 (MQ=255)
                  aaTCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCCGCGCGCGCGCGTTCTGCCGCATCTCTTCAg                                                                                                                         >  1:111492/1‑149 (MQ=255)
                  aaTCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAg                                                                                                                         >  3:22570/1‑149 (MQ=255)
                  aaTCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCA‑‑GCGCGCGCGTTCTGCCGCATCTCTTCAGt                                                                                                                        <  5:127848/148‑1 (MQ=255)
                    tCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGt                                                                                                                        >  1:24850/1‑148 (MQ=255)
                       gatagatGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAAt                                                                                                                    >  5:312268/1‑149 (MQ=255)
                         tagatGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCA‑‑GCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTggg                                                                                                                <  2:359686/149‑1 (MQ=255)
                           gatGCGCCCGCGCGGCCCGTCCAGGTCAACGACCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACCGGCGCGCAGGATGCTGGCGGCTTCCAGCGCGCGCGCCGTCTGCCGCATCTCTTCAGTTAATTg                                                                                                                  <  8:315138/147‑1 (MQ=255)
                            atGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGt                                                                                                               <  6:100567/149‑1 (MQ=255)
                                 cccGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCa                                                                                                          >  8:52058/1‑149 (MQ=255)
                                 cccGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCa                                                                                                          <  2:306091/149‑1 (MQ=255)
                                  ccGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCa                                                                                                          >  3:252251/1‑148 (MQ=255)
                                        ggCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTg                                                                                                   <  2:24850/149‑1 (MQ=255)
                                         gCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGc                                                                                                  >  2:250427/1‑149 (MQ=255)
                                         gCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGc                                                                                                  >  5:254521/1‑149 (MQ=255)
                                                   ttCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCTGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGt                                                                                        >  5:57817/1‑149 (MQ=255)
                                                   ttCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGc                                                                                                  >  4:153610/1‑139 (MQ=255)
                                                   ttCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGc                                                                                                  <  3:153610/139‑1 (MQ=255)
                                                      aTCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTAAATTGGGTTTGCATAGCCTGCTACTCAAAGTTcc                                                                                     >  7:4158/1‑149 (MQ=255)
                                                       tCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCcg                                                                                    >  2:367318/1‑149 (MQ=255)
                                                       tCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCcg                                                                                    >  5:100382/1‑149 (MQ=255)
                                                         gCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCcgcg                                                                                  <  7:228465/149‑1 (MQ=255)
                                                          cccAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCcgcg                                                                                  >  5:141236/1‑148 (MQ=255)
                                                             aGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTTGGTTTGCATAGCCTGCGACTCAAAGTTCCGCGTAcc                                                                              >  8:269912/1‑148 (MQ=255)
                                                              gAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACAt                                                                             <  6:141236/149‑1 (MQ=255)
                                                              gAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACAt                                                                             <  1:201975/149‑1 (MQ=255)
                                                               aaGCTGATAGGTTGGGCAGGTTTCGGTACAAAATCCGCAGTGAACAAAGGCGCGCAGGATGCTTTCGGCTTACAGCGCGCGCGAGTTCTGCCGCATATCTTCAGTTAATTGGGTTTGCATAGACTGATAATCAAAGTTCCTCGTACATg                                                                            >  8:101651/1‑149 (MQ=255)
                                                                   tGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGAc                                                                        >  7:283320/1‑149 (MQ=255)
                                                                    gATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCATCGCGCGCGCGTTCTGCCTCATCTCATCAGTTAATTGTGTTTACATAGCCTGCAACTCAAAGTACCGAGTACATgata                                                                         >  3:81470/1‑144 (MQ=255)
                                                                    gATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCg                                                                      >  5:39922/1‑149 (MQ=255)
                                                                     aTAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCg                                                                      >  7:302320/1‑149 (MQ=255)
                                                                       aGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCggg                                                                    >  1:311430/1‑149 (MQ=255)
                                                                        ggTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCgggg                                                                   <  3:116191/149‑1 (MQ=255)
                                                                             ggCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTaaa                                                              >  6:103502/1‑149 (MQ=255)
                                                                             ggCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTaaa                                                              <  5:103502/149‑1 (MQ=255)
                                                                              gCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCA‑‑GCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCa                                                                                                          <  6:326617/102‑1 (MQ=255)
                                                                              gCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCA‑‑GCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCa                                                                                                          >  5:326617/1‑102 (MQ=255)
                                                                               cAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCA‑‑GCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGGTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACAcg                                                          >  7:343464/1‑149 (MQ=255)
                                                                                  gTTGCGGTACAAAATCCGCAGTGAACCCAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCCTCTCTGCAGGTAAGTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGCCCGGGGTTAAACAcgc                                                         >  5:123289/1‑149 (MQ=255)
                                                                                   ttGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCCGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTTATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACAcgcc                                                        >  4:248432/1‑149 (MQ=255)
                                                                                           cAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAggg                                                 <  4:20853/148‑1 (MQ=255)
                                                                                           cAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGt                                                >  5:256716/1‑149 (MQ=255)
                                                                                            aaaaTTCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGc                                                                                                      >  3:55615/1‑94 (MQ=255)
                                                                                            aaaaTTCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGc                                                                                                      <  4:55615/94‑1 (MQ=255)
                                                                                            aaaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCaaa                                                                                          <  1:83330/106‑1 (MQ=255)
                                                                                            aaaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCaaa                                                                                          >  2:83330/1‑106 (MQ=255)
                                                                                            aaaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTc                                               >  3:100622/1‑149 (MQ=255)
                                                                                             ataTCCGCAGTGAACACAGGCGCGCAGGATGATGTCGGCTTCCAGCGCGCGCGAGTTCAGACGCATCTCGTCAGTTAATAGGGTTTACATACAATGCTCGTCAAAGGTCCGCGTACATTCGTCAGGCTTAAAAAACGCCGCAACGGTCg                                              >  2:337860/3‑149 (MQ=255)
                                                                                             aaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTc                                                                                      >  7:96011/1‑109 (MQ=255)
                                                                                             aaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTc                                                                                      <  8:96011/109‑1 (MQ=255)
                                                                                             aaaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTc                                               >  8:253060/1‑148 (MQ=255)
                                                                                              aaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTc                                               >  6:80884/1‑147 (MQ=255)
                                                                                              aaTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTc                                               >  4:210183/1‑147 (MQ=255)
                                                                                               aTCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAg                                            <  2:311430/149‑1 (MQ=255)
                                                                                                tCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTACTTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAgc                                           <  1:250427/149‑1 (MQ=255)
                                                                                                   gCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAgctgc                                        >  1:331342/1‑149 (MQ=255)
                                                                                                          acacAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTaa                                 <  8:270072/149‑1 (MQ=255)
                                                                                                              aGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAgctg                             >  1:97005/1‑149 (MQ=255)
                                                                                                               ggCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAgctgc                            <  8:283320/149‑1 (MQ=255)
                                                                                                                gcgcgcAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAgctgct                           <  4:344279/149‑1 (MQ=255)
                                                                                                                   cgcAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGt                        >  3:190132/1‑149 (MQ=255)
                                                                                                                    gcAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTg                       >  5:157845/1‑149 (MQ=255)
                                                                                                                        gATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAg                   >  1:203813/1‑149 (MQ=255)
                                                                                                                          tGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCg                 >  4:154303/1‑149 (MQ=255)
                                                                                                                          tGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCg                 >  3:55384/1‑149 (MQ=255)
                                                                                                                            cTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGa               <  7:52058/149‑1 (MQ=255)
                                                                                                                             tGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACCCGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGaa              >  2:103661/1‑149 (MQ=255)
                                                                                                                              gTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAg                                >  2:213411/1‑130 (MQ=255)
                                                                                                                                 ggCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTGTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGGTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATaaa          <  8:371184/149‑1 (MQ=255)
                                                                                                                                 ggCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTAAATTGGGTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAg         <  6:39922/149‑1 (MQ=255)
                                                                                                                                 ggCTTCCAGCGCGCGCGCGTTCTGCCGCAGCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCGGGTGATAGCGGAATaaa          <  4:252251/149‑1 (MQ=255)
                                                                                                                                  gCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAg         <  6:100382/149‑1 (MQ=255)
                                                                                                                                  gCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAg         >  4:160984/1‑149 (MQ=255)
                                                                                                                                   cTTCCCGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAg         >  6:320542/1‑148 (MQ=255)
                                                                                                                                         gcgcgcgcgcGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGa  <  3:154303/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
GCCTTCCAGCACCTGTTTAATCAGATAGATGCGCCCGCGCGGCCCGTCCAGTTCATCGCCCAGAAGCTGATAGGTTGGGCAGGTTGCGGTACAAAATCCGCAGTGAACACAGGCGCGCAGGATGCTGTCGGCTTCCAGCGCGCGCGCGTTCTGCCGCATCTCTTCAGTTAATTGGGTTTGCATAGCCTGCTCCTCAAAGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGA  >  NC_000913/3125277‑3125562

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