Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,589 G→C 6.9% intergenic (‑80/+27) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5890GC 6.9% 386.2 / 17.3 116intergenic (‑80/+27)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (65/43);  new base C (4/4);  total (69/47)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

Only 100 of 118 total aligned reads displayed.

GTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGG  >  NC_000913/3439441‑3439733
                                                                                                                                                    |                                                                                                                                                
gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                                                                  >  5:151362/1‑149 (MQ=255)
gTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                                                                  <  3:21064/149‑1 (MQ=255)
      aaaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCggg                                                                                                                                            >  8:247435/1‑149 (MQ=255)
       aaCTCACCTTTCGCTTGCTCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGt                                                                                                                                           <  7:368314/149‑1 (MQ=255)
       aaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGt                                                                                                                                           >  1:116879/1‑149 (MQ=255)
        aCTCACCTTTAGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCAAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtt                                                                                                                                          >  3:270984/1‑149 (MQ=255)
           cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttt                                                                                                                                       <  5:191229/149‑1 (MQ=255)
            aCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                      <  3:155195/149‑1 (MQ=255)
            aCCGTTCGCTTGCGCTGCAGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                       <  2:212774/148‑1 (MQ=255)
             ccTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttttt                                                                                                                                     <  1:227393/149‑1 (MQ=255)
             ccTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                      >  1:89178/1‑148 (MQ=255)
             ccTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                      >  1:45481/1‑148 (MQ=255)
                    cTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGc                                                                                                                              >  1:221126/1‑149 (MQ=255)
                         gcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCtt                                                                                                                            >  8:170272/1‑146 (MQ=255)
                          cTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGa                                                                                                                        <  4:123579/149‑1 (MQ=255)
                           ttCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAATAGAGTACCGCTATTGTAGATAAGTGGTGATACTAGGGGTATAAAAAAACCCGCAGGGGCGGGTTTTTTTACGTTGCTTCAGag                                                                                                                       >  5:91122/1‑148 (MQ=255)
                           ttCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                       >  4:185202/1‑133 (MQ=255)
                           ttCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                       <  3:185202/133‑1 (MQ=255)
                             cTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTa                                                                                                                     >  7:170641/1‑149 (MQ=255)
                                cTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActc                                                                                                                  <  3:258974/149‑1 (MQ=255)
                                  atatGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTg                                                                                                                >  8:361546/1‑149 (MQ=255)
                                               cccACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTc                                                                                                      <  5:191605/146‑1 (MQ=255)
                                               cccACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTc                                                                                                      <  8:237562/146‑1 (MQ=255)
                                               cccACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGc                                                                                                   >  5:116003/1‑149 (MQ=255)
                                                  aCTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                      <  3:325859/111‑1 (MQ=255)
                                                  aCTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGAAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                                                                >  7:119757/1‑149 (MQ=255)
                                                   cTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACg                                                                                                                                  <  1:188635/114‑1 (MQ=255)
                                                   cTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACg                                                                                                                                  >  2:188635/1‑114 (MQ=255)
                                                    tGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTg                                                                                                                               >  5:14050/1‑116 (MQ=255)
                                                      gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtct                                                                                              <  4:315562/147‑1 (MQ=255)
                                                      gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTctct                                                                                            >  3:284893/1‑149 (MQ=255)
                                                         aaGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                         >  5:331046/1‑149 (MQ=255)
                                                          aGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                        <  7:170272/149‑1 (MQ=255)
                                                          aGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                        >  5:54970/1‑149 (MQ=255)
                                                           gTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                        >  8:325968/1‑148 (MQ=255)
                                                              aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATc                                                                                    <  1:173906/149‑1 (MQ=255)
                                                               ccACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATc                                                                                    >  8:316059/1‑148 (MQ=255)
                                                               ccACATAAATAATCCAGGGGATGAACAAAAAGAGTACATCTATTGTAGATCAGTGGGGATACTACGGGTAAAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCa                                                                                   <  6:359265/149‑1 (MQ=255)
                                                                cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCaa                                                                                  >  3:237351/1‑149 (MQ=255)
                                                                cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCaa                                                                                  <  6:221980/149‑1 (MQ=255)
                                                                  caTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAc                                                                                 >  1:116518/1‑148 (MQ=255)
                                                                   aTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCa                                                                               <  6:54970/149‑1 (MQ=255)
                                                                   aTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCa                                                                               <  2:45481/149‑1 (MQ=255)
                                                                       aTACACCAGGGGATGAACAAACAGAGTACAGCTATTGTAGATAAGGGGGGATACTACGGGTATAAACAAACCCGCCGGGGCGGGTTTTTTTACGGTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTc                                                                           <  2:258806/149‑1 (MQ=255)
                                                                        tACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTc                                                                                                      >  6:14486/1‑121 (MQ=255)
                                                                          cTCCAGGGGATGAACAAAAAGAGTACAGCTATTGAAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGa                                                                         <  8:119757/148‑1 (MQ=255)
                                                                           tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATg                                                                       >  7:176288/1‑149 (MQ=255)
                                                                           tCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAACTGGGGATACTACGGGTATAAAAAAACCCGCAGGGTCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATg                                                                       >  1:261517/1‑149 (MQ=255)
                                                                             cAGGGGATGAACAAAAAGAGTACAGCTATTGTAGAGAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGCTTACTATGCAGCAGCTTCTGCTGTCTCTGAACGATCAACCAGCTCGATGt                                                                      >  8:332383/1‑148 (MQ=255)
                                                                               ggggATGCACAAAAAGAGTACAGCTATTGTAGATAAGTGTGTATACTACGGGTATAAAAAAACCCGCCGGTGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                   <  4:270984/149‑1 (MQ=255)
                                                                               ggggATGAACAAAAAGAGTTGAGCTTTTGTTGATAAGTTGGGATACTTCGGGTATAAATAATCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                   <  3:119502/149‑1 (MQ=255)
                                                                                     gaagaaAAAGAGTACATCTATTGTTGATATGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCATCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                    <  2:281065/142‑1 (MQ=255)
                                                                                     gAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTATGAACGATCAACCAGATCGATGTaa                                                                    >  1:281065/1‑142 (MQ=255)
                                                                                       aCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCAGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc                                                           >  1:166407/1‑149 (MQ=255)
                                                                                        cAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGATTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcg                                                          >  6:31482/1‑149 (MQ=255)
                                                                                          aaaaGAGTACAGCTATAGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGATTCAGATTACTCTGACTCAGATTCTCATTTCTCTGAACGAACAAACAGCTCGATGTAACCCATCcgcgcg                                                        >  8:106172/1‑149 (MQ=255)
                                                                                           aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGAt                                                                        >  8:351199/1‑132 (MQ=255)
                                                                                           aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGAt                                                                        <  7:351199/132‑1 (MQ=255)
                                                                                           aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGt                                                       >  5:57548/1‑149 (MQ=255)
                                                                                               agTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTc                                                   <  7:316059/149‑1 (MQ=255)
                                                                                                    aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCt                                               >  3:98528/1‑148 (MQ=255)
                                                                                                    aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTg                                              >  1:265384/1‑149 (MQ=255)
                                                                                                    aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                       >  8:263625/1‑60 (MQ=37)
                                                                                                    aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                       <  7:263625/60‑1 (MQ=37)
                                                                                                     gCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTg                                              >  8:215127/1‑148 (MQ=255)
                                                                                                         ttGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACg                                          >  4:113820/1‑148 (MQ=255)
                                                                                                            tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGaa                                       <  7:361546/148‑1 (MQ=255)
                                                                                                            tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGaa                                       <  7:325968/148‑1 (MQ=255)
                                                                                                            tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAg                                      >  5:16654/1‑149 (MQ=255)
                                                                                                             aGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCATCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACTCTCAACACGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGc                                     >  2:278152/1‑149 (MQ=255)
                                                                                                              gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                        <  3:268159/97‑1 (MQ=255)
                                                                                                              gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                        >  4:268159/1‑97 (MQ=255)
                                                                                                                  aGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCcaca                                 <  3:185854/148‑1 (MQ=255)
                                                                                                                      gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACtt                              >  8:25572/1‑147 (MQ=255)
                                                                                                                       ggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                           <  3:323497/149‑1 (MQ=255)
                                                                                                                       ggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGGTGTCGCCTGCACGGAAGCCACACtt                              >  7:83785/1‑146 (MQ=255)
                                                                                                                        gATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCACCCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                           >  7:98273/1‑148 (MQ=255)
                                                                                                                        gATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAg                           >  5:99994/1‑148 (MQ=255)
                                                                                                                        gATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCAAGGAAGCCACACTTCAg                           >  5:142178/1‑148 (MQ=255)
                                                                                                                          tactacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTTAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAAt                        <  4:98528/149‑1 (MQ=255)
                                                                                                                           actacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                    >  5:183124/1‑104 (MQ=255)
                                                                                                                           actacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                    <  6:183124/104‑1 (MQ=255)
                                                                                                                             tacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCAACCTTCAGAATAc                      >  5:201410/1‑148 (MQ=255)
                                                                                                                             tacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTGTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATAc                      >  8:396856/1‑148 (MQ=255)
                                                                                                                              acGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATAc                      >  5:381559/1‑147 (MQ=255)
                                                                                                                               cGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAg                   <  6:151362/149‑1 (MQ=255)
                                                                                                                                gggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgt                  <  4:237351/149‑1 (MQ=255)
                                                                                                                                  gTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtg                 >  6:348596/1‑148 (MQ=255)
                                                                                                                                   cgtaAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtgt                >  7:104047/3‑148 (MQ=255)
                                                                                                                                        aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTTCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGATTACGAGTGTAaccac          >  1:123975/1‑149 (MQ=255)
                                                                                                                                        aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac          >  6:369464/1‑149 (MQ=255)
                                                                                                                                        aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTGCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTGCAGAATACGCGTGGAacccc          >  1:323802/1‑147 (MQ=255)
                                                                                                                                          aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCCACCAGCTCGATGTACGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCCCACGTCAGAATACGAGTGTAACCCCCg        >  4:196900/1‑149 (MQ=255)
                                                                                                                                          aaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCg        >  6:360034/1‑149 (MQ=255)
                                                                                                                                           aaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGTTTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCgg       <  2:261517/149‑1 (MQ=255)
                                                                                                                                           aaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                                                                             <  4:90527/63‑1 (MQ=255)
                                                                                                                                           aaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                                                                             >  3:90527/1‑63 (MQ=255)
                                                                                                                                             cccGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCa     <  6:16654/149‑1 (MQ=255)
                                                                                                                                              ccgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                         >  4:131932/1‑96 (MQ=255)
                                                                                                                                                gccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACgg  >  1:252539/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGG  >  NC_000913/3439441‑3439733

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