Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,593,090 T→C A478A (GCT→GCC tsr → methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,593,0900TC99.7% 1268.3 / ‑5.4 307A478A (GCT→GCCtsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (173/133);  total (173/134)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

GTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAG  >  NC_000913/4592992‑4593187
                                                                                                  |                                                                                                 
 tCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCgag                                                                                                <  1:8499535/101‑1 (MQ=255)
  cAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCgaga                                                                                               >  2:8499945/1‑101 (MQ=255)
   aGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGAt                                                                                              >  1:8499948/1‑101 (MQ=255)
      gcCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGa                                                                                           >  2:8499952/1‑101 (MQ=255)
       ccGTGACCCGCGTGACGGACATTATGGTCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGAc                                                                                          >  1:8499953/1‑101 (MQ=255)
       ccGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGAc                                                                                          >  1:8499956/1‑101 (MQ=255)
         gTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCg                                                                                        >  1:8499961/1‑101 (MQ=255)
         gTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCg                                                                                        >  2:8499962/1‑101 (MQ=255)
          tGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCgg                                                                                       <  2:8499547/101‑1 (MQ=255)
             cccGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCggg                                                                                      >  2:8499966/1‑99 (MQ=255)
             cccGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTa                                                                                    <  1:8499385/101‑1 (MQ=255)
             cccGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTa                                                                                    <  2:8499763/101‑1 (MQ=255)
              ccGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTaa                                                                                   >  1:8499968/1‑101 (MQ=255)
               cgcgTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAAc                                                                                  >  1:8499974/1‑101 (MQ=255)
               cgcgTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAAc                                                                                  >  2:8499973/1‑101 (MQ=255)
                 cgTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTc                                                                                <  2:8499581/101‑1 (MQ=255)
                  gTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCa                                                                               <  2:8499694/101‑1 (MQ=255)
                  gTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCa                                                                               <  2:8499350/101‑1 (MQ=255)
                   tGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCaa                                                                              >  2:8499980/1‑101 (MQ=255)
                   tGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCaa                                                                              >  2:8499983/1‑101 (MQ=255)
                    gACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAAc                                                                             <  1:8499495/101‑1 (MQ=255)
                     aCGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACa                                                                            >  2:8499988/1‑101 (MQ=255)
                      cGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAg                                                                           <  1:8499725/101‑1 (MQ=255)
                      cGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAg                                                                           >  1:8499989/1‑101 (MQ=255)
                       ggACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGAt                                                                                              <  1:8499655/81‑1 (MQ=255)
                       ggACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGa                                                                          <  1:8499406/101‑1 (MQ=255)
                        gACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGaa                                                                         >  2:8499996/1‑101 (MQ=255)
                        gACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGaa                                                                         <  1:8499234/101‑1 (MQ=255)
                        gACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGaa                                                                         >  1:8499994/1‑101 (MQ=255)
                         aCATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAAc                                                                        <  2:8499291/101‑1 (MQ=255)
                         aCATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAAc                                                                        >  2:8499999/1‑101 (MQ=255)
                         aCATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAAc                                                                        <  1:8499331/101‑1 (MQ=255)
                           aTTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAAtgc                                                                      >  1:8500003/1‑101 (MQ=255)
                           aTTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAAcgc                                                                      <  2:8499643/101‑1 (MQ=255)
                             tATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAAcgccg                                                                    >  2:8500006/1‑101 (MQ=255)
                               tGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCcgcg                                                                  <  2:8499671/101‑1 (MQ=255)
                               tGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCcgcg                                                                  <  2:8499471/101‑1 (MQ=255)
                                gggCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCgcgc                                                                 >  2:8500010/1‑101 (MQ=255)
                                gggCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCgcgc                                                                 >  1:8500012/1‑101 (MQ=255)
                                 ggCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCt                                                                >  1:8500016/1‑101 (MQ=255)
                                  gCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCtg                                                               >  2:8500017/1‑101 (MQ=255)
                                    gAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCtggt                                                             >  2:8500022/1‑101 (MQ=255)
                                      aaTTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCtggtgg                                                           >  2:8500025/1‑101 (MQ=255)
                                       aTTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGa                                                          >  1:8500027/1‑101 (MQ=255)
                                        ttGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGaa                                                         <  2:8499409/101‑1 (MQ=255)
                                         tgcttctgcttcTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGaa                                                                         <  2:8499275/84‑1 (MQ=255)
                                         tgcttctgcttcTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAag                                                        <  2:8499584/101‑1 (MQ=255)
                                         tgcttctgcttcTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAag                                                        >  1:8500030/1‑101 (MQ=255)
                                          gcttctgcttctGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAaga                                                       <  1:8499672/101‑1 (MQ=255)
                                            ttctgcttctgATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGt                                                     >  2:8500033/1‑101 (MQ=255)
                                                gcttctgATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgc                                                 >  1:8500037/1‑101 (MQ=255)
                                                gcttctgATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgc                                                 >  2:8500038/1‑101 (MQ=255)
                                                 cttctgATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgcc                                                >  2:8500041/1‑101 (MQ=255)
                                                  ttctgATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgccg                                               <  1:8499682/101‑1 (MQ=255)
                                                   tctgATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgccgc                                              >  2:8500046/1‑101 (MQ=255)
                                                   tctgATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgccgc                                              <  1:8499521/101‑1 (MQ=255)
                                                    ctgATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgccgct                                             >  1:8500049/1‑101 (MQ=255)
                                                        tgaGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgccgctgccg                                         >  1:8500053/1‑101 (MQ=255)
                                                        tgaGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgccgctgc                                           >  2:8500055/1‑99 (MQ=255)
                                                         gaGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCtgccgctgccgc                                        <  2:8499597/101‑1 (MQ=255)
                                                          aGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCga                                                                                                 <  2:8499514/43‑1 (MQ=255)
                                                          aGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTgccgcc                                       >  2:8500059/1‑101 (MQ=255)
                                                            cAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTgccgccgc                                     >  2:8500062/1‑101 (MQ=255)
                                                            cAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTgccgccgc                                     <  2:8499533/101‑1 (MQ=255)
                                                            cAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTgccgccgc                                     >  1:8500065/1‑101 (MQ=255)
                                                             agagCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTgccgccgcc                                    >  2:8500069/1‑101 (MQ=255)
                                                              gagCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGccgccgccg                                   <  2:8499718/101‑1 (MQ=255)
                                                               agCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCtg                                                               >  1:8500072/1‑72 (MQ=255)
                                                                gCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCcgcg                                 >  2:8500074/1‑101 (MQ=255)
                                                                 ccGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCgcgc                                >  1:8500076/1‑101 (MQ=255)
                                                                  cGTGGTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCt                               >  2:8500079/1‑101 (MQ=255)
                                                                     ggTATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGa                            <  2:8499569/101‑1 (MQ=255)
                                                                       tATCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAag                          <  2:8499442/101‑1 (MQ=255)
                                                                        aTCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAaga                         <  1:8499772/101‑1 (MQ=255)
                                                                        aTCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAaga                         >  2:8500084/1‑101 (MQ=255)
                                                                         tCGATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAagag                        >  1:8500087/1‑101 (MQ=255)
                                                                           gATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCa                      >  2:8500088/1‑101 (MQ=255)
                                                                           gATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCa                      >  2:8500091/1‑101 (MQ=255)
                                                                           gATCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCa                      <  1:8499451/101‑1 (MQ=255)
                                                                            aTCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTctgcc                                                >  2:8500098/1‑74 (MQ=255)
                                                                            aTCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAg                     >  2:8500095/1‑101 (MQ=255)
                                                                            aTCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAg                     <  1:8499637/101‑1 (MQ=255)
                                                                             tCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCt                               <  1:8499508/90‑1 (MQ=255)
                                                                             tCAGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAgg                    <  2:8499558/101‑1 (MQ=255)
                                                                               aGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGcc                  <  1:8499457/101‑1 (MQ=255)
                                                                               aGGTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGcc                  <  1:8499435/101‑1 (MQ=255)
                                                                                ggTTGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCa                 >  1:8500100/1‑101 (MQ=255)
                                                                                  ttGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGt               >  2:8500105/1‑101 (MQ=255)
                                                                                   tGGCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTc              <  1:8499550/101‑1 (MQ=255)
                                                                                     gCTTAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTg                              >  1:8500108/1‑83 (MQ=255)
                                                                                        tAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTg         >  2:8500110/1‑101 (MQ=255)
                                                                                        tAGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTg         <  2:8499634/101‑1 (MQ=255)
                                                                                         aGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGa        <  2:8499605/101‑1 (MQ=255)
                                                                                         aGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGa        <  2:8499661/101‑1 (MQ=255)
                                                                                         aGCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGa        <  1:8499620/101‑1 (MQ=255)
                                                                                          gCGGTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGAc       <  2:8499432/101‑1 (MQ=255)
                                                                                            ggTTGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCg     >  2:8500114/1‑101 (MQ=255)
                                                                                              ttGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGaa   >  1:8500118/1‑101 (MQ=255)
                                                                                              ttGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGaa   <  1:8499602/101‑1 (MQ=255)
                                                                                              ttGCCGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGaa   <  1:8499559/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
GTCAGCGCCGTGACCCGCGTGACGGACATTATGGGCGAAATTGCTTCTGCTTCTGATGAGCAGAGCCGTGGTATCGATCAGGTTGGCTTAGCGGTTGCTGAGATGGACCGGGTAACTCAACAGAACGCCGCGCTGGTGGAAGAGTCTGCCGCTGCCGCCGCCGCGCTGGAAGAGCAGGCCAGTCGCCTGACCGAAG  >  NC_000913/4592992‑4593187

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