Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,286 G→T 6.1% E66D (GAG→GAT) ‡ 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,2860GT 6.1% 290.0 / 4.9 82E66D (GAG→GAT) ‡tonBmembrane spanning protein in TonB‑ExbB‑ExbD transport complex
Reads supporting (aligned to +/- strand):  ref base G (36/41);  new base T (3/2);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAAAC  >  NC_000913/1311138‑1311434
                                                                                                                                                    |                                                                                                                                                    
tCTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCCGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAg                                                                                                                                                      >  2:503066/1‑149 (MQ=255)
 cTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGGTATTTAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGGTGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGc                                                                                                                                                     >  1:439577/1‑149 (MQ=255)
   gCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCg                                                                                                                                                   <  2:64985/149‑1 (MQ=255)
     aTTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggt                                                                                                                                                 <  3:154579/149‑1 (MQ=255)
     aTTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggt                                                                                                                                                 <  2:439577/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                                                                                                                   >  2:657135/1‑149 (MQ=255)
                                   tataCCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGAt                                                                                                                   >  3:713099/1‑149 (MQ=255)
                                      aCCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATccc                                                                                                                >  1:330986/1‑149 (MQ=255)
                                        cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                              >  4:232101/1‑149 (MQ=255)
                                        cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                              >  2:620690/1‑149 (MQ=255)
                                        cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                              >  4:324044/1‑149 (MQ=255)
                                        cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                              >  4:593444/1‑149 (MQ=255)
                                        cTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                              >  3:217051/1‑149 (MQ=255)
                                         tCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGa                                                                                                             >  4:735525/1‑149 (MQ=255)
                                         tCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGa                                                                                                             >  2:195261/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                                                                             <  4:797887/141‑1 (MQ=255)
                                                                                 cGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGaaa                                                                             >  3:797887/1‑141 (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                                                         >  1:57949/1‑128 (MQ=255)
                                                                                                                  tcGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAg                                                         <  2:57949/128‑1 (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                                                               <  2:155812/105‑1 (MQ=255)
                                                                                                                                   gTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGCGGTCATTGAAAAGCCGAAGCCGaaa                                                               >  1:155812/1‑105 (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          <  1:495063/149‑1 (MQ=255)
                                                                                                                                            ccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGa          <  2:294174/149‑1 (MQ=255)
                                                                                                                                               ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCggtgg                                                                                      <  1:673364/70‑1 (MQ=255)
                                                                                                                                               ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCggtgg                                                                                      >  2:673364/1‑70 (MQ=255)
                                                                                                                                               ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGGAGCCTACGCCAAAGCCGACGCCGGTGAAAAAGGTACACGAGCAGCCAAATCGCGATGt       <  2:748855/149‑1 (MQ=255)
                                                                                                                                               ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGt       <  4:217051/149‑1 (MQ=255)
                                                                                                                                                  gAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCAGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGTGCAGCCAAAACGCGATGTCaa    <  2:47420/149‑1 (MQ=255)
                                                                                                                                                    gCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAAAc  <  3:593444/149‑1 (MQ=255)
                                                                                                                                                    gCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAAAc  <  3:324044/149‑1 (MQ=255)
                                                                                                                                                    gCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAAAc  <  3:232101/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCAAAC  >  NC_000913/1311138‑1311434

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