Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,592,522 A→T 5.2% N289I (AAC→ATC)  tsr → methyl‑accepting chemotaxis protein I, serine sensor receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,5220AT 5.2% 227.4 / 2.7 77N289I (AAC→ATC) tsrmethyl‑accepting chemotaxis protein I, serine sensor receptor
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base T (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

TTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTAAGTGCTTCTGAAACGGCG  >  NC_000913/4592378‑4592667
                                                                                                                                                |                                                                                                                                                 
ttGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGAt                                                                                                                                               <  1:246227/149‑1 (MQ=255)
      tGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATctctct                                                                                                                                         >  6:101617/1‑149 (MQ=255)
       ggATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCtt                                                                                                                                        >  3:155679/1‑149 (MQ=255)
        gATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCtt                                                                                                                                        <  7:19463/148‑1 (MQ=255)
        gATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTc                                                                                                                                       >  2:43821/1‑149 (MQ=255)
            gCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTcgc                                                                                                                                     >  5:62818/1‑147 (MQ=255)
            gCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTcgc                                                                                                                                     <  6:62818/147‑1 (MQ=255)
              tctAATGAGATGGGGCAACTGGCCGAGCGTTTGCGGCAGATGCTGGGAGAGCTGATGCGTACCGTCTGTGATGTGCGCGCGGGGGCGAATGCCATCGATAGCGGTGCCAGCGAAATCGACCCCGGCAATAACGGTCTCGGTTCGAGCAc                                                                                                                                 >  4:35962/1‑149 (MQ=255)
              tctAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGCGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCAGCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCgcgccc                                                                                                                                 >  2:85750/1‑147 (MQ=255)
              tctAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCAc                                                                                                                                 >  2:311999/1‑149 (MQ=255)
                  aTGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCAcc                                                                                                                                <  2:57075/146‑1 (MQ=255)
                  aTGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAg                                                                                                                             <  3:196196/149‑1 (MQ=255)
                    gagaTGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGc                                                                                                                            >  3:2664/1‑148 (MQ=255)
                     agaTGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCaa                                                                                                                          <  4:164935/149‑1 (MQ=255)
                      gaTGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCCACGGGGCCCATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCaa                                                                                                                          >  1:240785/1‑148 (MQ=255)
                            gCAACTGGCAGAGAGTTTGCGCCATATGCAGTGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGAGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAgtc                                                                                                                     >  5:73962/1‑145 (MQ=255)
                             cAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGc                                                                                                                  <  7:184114/149‑1 (MQ=255)
                                cTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCTGTGATGTGCGCAACGGGGCCAATGCCATCGATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTc                                                                                                               <  5:24399/149‑1 (MQ=255)
                                 tGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCg                                                                                                              >  5:68978/1‑149 (MQ=255)
                                 tGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCg                                                                                                              >  3:288134/1‑149 (MQ=255)
                                  ggCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCACCGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGc                                                                                                             >  4:15791/1‑149 (MQ=255)
                                    cAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTg                                                                                                           <  1:43821/149‑1 (MQ=255)
                                     agagagTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATATCGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTg                                                                                                           >  7:30360/1‑148 (MQ=255)
                                     agagagTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATATCGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTg                                                                                                           >  7:35419/1‑148 (MQ=255)
                                       agagTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGaa                                                                                                        <  2:221727/149‑1 (MQ=255)
                                         agTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGc                                                                                                                     >  4:8714/1‑134 (MQ=255)
                                         agTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGc                                                                                                                     <  3:8714/134‑1 (MQ=255)
                                          gTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAaga                                                                                                      <  3:16693/148‑1 (MQ=255)
                                           tttGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAgaga                                                                                                    >  6:295887/1‑149 (MQ=255)
                                             tGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGAc                                                                                                   <  6:309228/148‑1 (MQ=255)
                                                gcCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGca                                                                                               <  4:155679/149‑1 (MQ=255)
                                                 ccATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGcag                                                                                              >  3:56229/1‑149 (MQ=255)
                                                  cATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGcagc                                                                                             <  4:288134/149‑1 (MQ=255)
                                                    tatGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGcagcca                                                                                           <  7:52044/149‑1 (MQ=255)
                                                         aGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATg                                                                                      <  4:38819/149‑1 (MQ=255)
                                                         aGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATg                                                                                      <  6:164873/149‑1 (MQ=255)
                                                         aGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATg                                                                                      <  6:68978/149‑1 (MQ=255)
                                                            gagaGCTGATGCGTACCGTCGGTGATGGGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGTCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAg                                                                                   <  5:101617/149‑1 (MQ=255)
                                                             agagCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAg                                                                                   >  7:318450/1‑148 (MQ=255)
                                                                   gATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGTGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGcagcca                                                                                           <  4:151637/134‑1 (MQ=255)
                                                                   gATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGTGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGcagcca                                                                                           >  3:151637/1‑134 (MQ=255)
                                                                    aTGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGTCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGAc                                                                           <  5:138093/149‑1 (MQ=255)
                                                                    aTGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGAc                                                                           >  4:283110/1‑149 (MQ=255)
                                                                       cGTACCGTCGCTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGAGCCAGCGAAATCGCCCCCGGCAATCACGATCTCTCTGCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGCCGGCAGCCAGCATGGAGCAACTGACCg                                                                         <  3:35962/148‑1 (MQ=255)
                                                                             gTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGt                                                                  <  4:56229/149‑1 (MQ=255)
                                                                             gTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGAGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGt                                                                  <  4:128304/149‑1 (MQ=255)
                                                                                    aTGTGCGCAACGGGGCCAATGCCATCTATAGCGGGGCCAGCGACATCGCCACCGCCTATAACGATCGCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAg                                                           <  2:306881/149‑1 (MQ=255)
                                                                                       tgCGCAACGGGGCCAAATAGATCTTAAGAGGTGCCAGCTAAATCGACACAGTCAATAAAGATCTCTCTTCGGGCACCGAGCAACAGGCCGATTAGCTGGAAGAGACGGCAGCCAGCATGGAGGAACTGCGCGGAACGGTGAAGCAGAAc                                                        >  4:226517/1‑149 (MQ=255)
                                                                                          gcAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCCCCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGcc                                                     >  4:321527/1‑149 (MQ=255)
                                                                                                    cAATGCCATCTATAGCGGTGCCAGCGAAATCGGCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATgcgc                                           <  7:42206/149‑1 (MQ=255)
                                                                                                    cAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGAAGACCGCCGAAAAtcc                                             >  6:83641/1‑145 (MQ=255)
                                                                                                    cAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCAGGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATgc                                             <  5:83641/147‑1 (MQ=255)
                                                                                                        gCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATgcgcgc                                         >  6:90432/1‑147 (MQ=255)
                                                                                                         ccATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCgccag                                      >  5:68842/1‑149 (MQ=255)
                                                                                                         ccATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCgccag                                      <  1:311999/149‑1 (MQ=255)
                                                                                                          cATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCTCAACGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCATCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCgccagg                                     <  6:148266/149‑1 (MQ=255)
                                                                                                            tCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCgccaggc                                    >  1:12374/1‑148 (MQ=255)
                                                                                                                taGCGGTGCCAGCGAAATCGCCACCGGCAATATCGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCgccaggcca                                  <  8:35419/146‑1 (MQ=255)
                                                                                                                taGCGGTGCCAGCGAAATCGCCACCGGCAATATCGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCgccaggcca                                  <  8:30360/146‑1 (MQ=255)
                                                                                                                taGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCgccaggcca                                  <  2:12374/146‑1 (MQ=255)
                                                                                                                      tGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTgg                         >  1:821/1‑149 (MQ=255)
                                                                                                                       gCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGc                        <  1:85750/149‑1 (MQ=255)
                                                                                                                        ccAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGGCATCTGGCg                       >  2:318958/1‑149 (MQ=255)
                                                                                                                        ccAGCGAAATCGCCACCGGCAATAACGAGCTCTCTTCGCGCACCGAGCAACa                                                                                                                        <  2:128269/52‑1 (MQ=255)
                                                                                                                         cAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGt                      >  7:80723/1‑149 (MQ=255)
                                                                                                                         cAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGt                      >  7:71139/1‑149 (MQ=255)
                                                                                                                          aGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGt                      >  2:196864/1‑148 (MQ=255)
                                                                                                                           gCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGtt                     >  6:16094/1‑148 (MQ=255)
                                                                                                                            cGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTaa                   <  5:295887/149‑1 (MQ=255)
                                                                                                                             gAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTAAg                  <  3:132664/149‑1 (MQ=255)
                                                                                                                             gAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTAAg                  <  2:821/149‑1 (MQ=255)
                                                                                                                                 tCGCCCCCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTAAGTGCt              >  8:2082/1‑149 (MQ=255)
                                                                                                                                 tCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAg                                                              >  8:129964/1‑101 (MQ=255)
                                                                                                                                      aCCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGCCGGCAGCCAGGATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTAAGTGCTTCTGa         <  2:240785/149‑1 (MQ=255)
                                                                                                                                            aaTAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTAAGTGCTTCTGAAACGgc   >  2:19889/1‑149 (MQ=255)
                                                                                                                                            aaTAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTAAGTGCTTCTGAAACGgc   <  3:283110/149‑1 (MQ=255)
                                                                                                                                             aTAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTAAGTGCTTCTGAAACGgcg  <  3:15791/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCGCACCGAGCAACAGGCCGCTTCGCTGGAAGAGACGGCAGCCAGCATGGAGCAACTGACCGCAACGGTGAAGCAGAACGCCGAAAATGCGCGCCAGGCCAGCCATCTGGCGTTAAGTGCTTCTGAAACGGCG  >  NC_000913/4592378‑4592667

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