Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,285 A→C 5.4% E66A (GAG→GCG) ‡ tonB → membrane spanning protein in TonB‑ExbB‑ExbD transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,2850AC 5.4% 288.8 / 6.9 83E66A (GAG→GCG) ‡tonBmembrane spanning protein in TonB‑ExbB‑ExbD transport complex
Reads supporting (aligned to +/- strand):  ref base A (37/41);  new base C (3/2);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

GTCTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAA  >  NC_000913/1311137‑1311432
                                                                                                                                                    |                                                                                                                                                   
gTCTGCATTCATGGTGCTGTTGTGGCGGGGCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGGAGCCGATTTTTGTTAAGATGGTTACGCCCGCTGATCTTGAACCCCCACAAGCCGTTCAGCCGCCACCGGa                                                                                                                                                     >  3:712077/1‑149 (MQ=255)
 tCTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCCGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAg                                                                                                                                                    >  2:503066/1‑149 (MQ=255)
  cTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGGTATTTAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGGTGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGc                                                                                                                                                   >  1:439577/1‑149 (MQ=255)
    gCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCg                                                                                                                                                 <  2:64985/149‑1 (MQ=255)
      aTTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggt                                                                                                                                               <  2:439577/149‑1 (MQ=255)
      aTTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggt                                                                                                                                               <  3:154579/149‑1 (MQ=255)
         cATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggtggt                                                                                                                                            <  2:397355/149‑1 (MQ=255)
         cATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggtggt                                                                                                                                            <  2:494457/149‑1 (MQ=255)
                cTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCa                                                                                                                                     <  1:357406/149‑1 (MQ=255)
                        gCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGa                                                                                                                             >  3:384736/1‑149 (MQ=255)
                        gCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGa                                                                                                                             >  1:748855/1‑149 (MQ=255)
                            gTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCt                                                                                                                         >  4:657285/1‑149 (MQ=255)
                              cTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGa                                                                                                                       <  1:516930/149‑1 (MQ=255)
                               tGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAg                                                                                                                      >  2:598555/1‑149 (MQ=255)
                               tGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAg                                                                                                                      >  1:847/1‑149 (MQ=255)
                                gCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGc                                                                                                                     >  4:384337/1‑149 (MQ=255)
                                 ctctATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGcc                                                                                                                    <  2:136214/149‑1 (MQ=255)
                                    tataCCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGAt                                                                                                                 >  3:713099/1‑149 (MQ=255)
                                    tataCCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGAt                                                                                                                 >  2:657135/1‑149 (MQ=255)
                                       aCCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATccc                                                                                                              >  1:330986/1‑149 (MQ=255)
                                         cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                            >  2:620690/1‑149 (MQ=255)
                                         cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                            >  3:217051/1‑149 (MQ=255)
                                         cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                            >  4:593444/1‑149 (MQ=255)
                                         cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                            >  4:232101/1‑149 (MQ=255)
                                         cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                            >  4:324044/1‑149 (MQ=255)
                                          tCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGa                                                                                                           >  2:195261/1‑149 (MQ=255)
                                          tCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGa                                                                                                           >  4:735525/1‑149 (MQ=255)
                                            ggTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAAc                                                                                                         <  4:171872/149‑1 (MQ=255)
                                               aCATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGTTTTCTGTCACGTTGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGc                                                                                                      <  3:384337/149‑1 (MQ=255)
                                                 aTCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCa                                                                                                    >  1:47420/1‑149 (MQ=255)
                                                    aGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCaaaa                                                                                                 >  3:696766/1‑149 (MQ=255)
                                                       ttATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAaagaa                                                                                              <  1:736275/149‑1 (MQ=255)
                                                              aCTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCgg                                                                                       <  1:503066/149‑1 (MQ=255)
                                                                tACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCggtg                                                                                     >  2:702732/1‑149 (MQ=255)
                                                                 aCCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCggtgg                                                                                    <  2:330986/149‑1 (MQ=255)
                                                                      cgcctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCAtt                                                                               <  3:657285/149‑1 (MQ=255)
                                                                        cctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGa                                                                             >  1:294174/1‑149 (MQ=255)
                                                                         ctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGaa                                                                            >  1:177932/1‑149 (MQ=255)
                                                                          tgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGaaa                                                                           <  4:713099/149‑1 (MQ=255)
                                                                           gcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGaaaa                                                                          <  2:847/149‑1 (MQ=255)
                                                                           gcgcAGCCGATTCCTTTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGaaaa                                                                          <  1:428387/149‑1 (MQ=255)
                                                                               aGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCACGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccg                                                                      >  4:759638/1‑149 (MQ=255)
                                                                                gCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggtggt                                                                                                                                            >  1:164762/1‑78 (MQ=255)
                                                                                gCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggtggt                                                                                                                                            <  2:164762/78‑1 (MQ=255)
                                                                                 ccGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccgaa                                                                    <  2:573825/149‑1 (MQ=255)
                                                                                  cGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGaaa                                                                           >  3:797887/1‑141 (MQ=255)
                                                                                  cGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGaaa                                                                           <  4:797887/141‑1 (MQ=255)
                                                                                   gATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccgaagc                                                                  >  4:671555/1‑149 (MQ=255)
                                                                                     tttCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccgaagccg                                                                >  3:646790/1‑149 (MQ=255)
                                                                                        cTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGaaa                                                             <  4:712077/149‑1 (MQ=255)
                                                                                          gTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAAcc                                                           <  1:598555/149‑1 (MQ=255)
                                                                                              cGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAg                                                       <  1:480579/149‑1 (MQ=255)
                                                                                                 tGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCa                                                    <  1:620690/149‑1 (MQ=255)
                                                                                                              tgATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCgg                                       >  4:149124/1‑149 (MQ=255)
                                                                                                                 tctcGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGa                                    <  3:213695/149‑1 (MQ=255)
                                                                                                                  ctcGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaa                                   >  4:251070/1‑149 (MQ=255)
                                                                                                                   tcGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAg                                                       <  2:57949/128‑1 (MQ=255)
                                                                                                                   tcGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAg                                                       >  1:57949/1‑128 (MQ=255)
                                                                                                                     taaCCGCCACAAGGCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACGGCCAAAAGAAGCACTGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTCaaaaa                                >  4:517721/2‑149 (MQ=255)
                                                                                                                     gAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaaaaa                                <  4:384736/149‑1 (MQ=255)
                                                                                                                      aaCCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAg                               <  3:735525/149‑1 (MQ=255)
                                                                                                                      aaCCGCCACAAGCCGTTCAGCCGCAACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAAACGCCAAAAGAAGCACCGGTGATCATTGAAATGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAg                               <  1:195261/149‑1 (MQ=255)
                                                                                                                        ccGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAg                               >  1:623872/1‑147 (MQ=255)
                                                                                                                        ccGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAg                               <  2:623872/147‑1 (MQ=255)
                                                                                                                         cGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTa                            <  1:657135/149‑1 (MQ=255)
                                                                                                                         cGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTa                            <  3:149124/149‑1 (MQ=255)
                                                                                                                           ccACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACa                          >  3:309125/1‑149 (MQ=255)
                                                                                                                              caAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGa                       <  4:456475/149‑1 (MQ=255)
                                                                                                                                    gTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAgg                              >  1:791396/1‑136 (MQ=255)
                                                                                                                                    gTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAgg                              <  2:791396/136‑1 (MQ=255)
                                                                                                                                    gTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGCGGTCATTGAAAAGCCGAAGCCGaaa                                                             >  1:155812/1‑105 (MQ=255)
                                                                                                                                    gTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGCGGTCATTGAAAAGCCGAAGCCGaaa                                                             <  2:155812/105‑1 (MQ=255)
                                                                                                                                         tctaccACCGGCTCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAAc            >  4:608400/5‑149 (MQ=255)
                                                                                                                                          cttccACCGGCTCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccga                                                                     >  3:242899/4‑91 (MQ=255)
                                                                                                                                          ctaccACCGGCTCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccga                                                                     <  4:242899/88‑1 (MQ=255)
                                                                                                                                           taccACCGGCTCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTa                            <  4:96060/129‑1 (MQ=255)
                                                                                                                                           taccACCGGCTCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTa                            >  3:96060/3‑131 (MQ=255)
                                                                                                                                             ccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGa        <  2:294174/149‑1 (MQ=255)
                                                                                                                                             ccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGa        <  1:495063/149‑1 (MQ=255)
                                                                                                                                                ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCggtgg                                                                                    >  2:673364/1‑70 (MQ=255)
                                                                                                                                                ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCggtgg                                                                                    <  1:673364/70‑1 (MQ=255)
                                                                                                                                                ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGGAGCCTACGCCAAAGCCGACGCCGGTGAAAAAGGTACACGAGCAGCCAAATCGCGATGt     <  2:748855/149‑1 (MQ=255)
                                                                                                                                                ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGt     <  4:217051/149‑1 (MQ=255)
                                                                                                                                                   gAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCAGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGTGCAGCCAAAACGCGATGTCaa  <  2:47420/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAA  >  NC_000913/1311137‑1311432

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