Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,706,607 A→C 8.9% L173R (CTC→CGC) ‡ lepA ← back‑translocating elongation factor EF4, GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,706,6070AC 8.9% 283.7 / 20.7 90L173R (CTC→CGC) ‡lepAback‑translocating elongation factor EF4, GTPase
Reads supporting (aligned to +/- strand):  ref base A (41/41);  new base C (4/4);  total (45/45)
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.54e-01

GTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCGA  >  NC_000913/2706462‑2706754
                                                                                                                                                 |                                                                                                                                                   
gTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGaa                                                                                                                                                  <  2:641552/149‑1 (MQ=255)
     ccTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTcc                                                                                                                                             <  1:589188/149‑1 (MQ=255)
      cTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctg                                                                                                                                            <  4:559365/149‑3 (MQ=255)
       ttACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctgg                                                                                                                                           <  3:494897/149‑4 (MQ=255)
           gCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGcaca                                                                                                                                       >  2:489897/1‑149 (MQ=255)
                  gCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGc                                                                                                                                <  2:671167/149‑1 (MQ=255)
                      tttttAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGt                                                                                                                            >  3:335962/1‑149 (MQ=255)
                        tttAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGttt                                                                                                                          <  2:896543/149‑1 (MQ=255)
                          tAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTc                                                                                                                        <  4:219198/149‑1 (MQ=255)
                            aTACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGc                                                                                                                      <  3:789887/149‑1 (MQ=255)
                            aTACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGc                                                                                                                      >  4:79490/1‑149 (MQ=255)
                               cGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGa                                                                                                                   <  3:318655/149‑1 (MQ=255)
                                      gTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGcc                                                                                                                               >  1:769040/1‑130 (MQ=255)
                                      gTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGcc                                                                                                                               <  2:769040/130‑1 (MQ=255)
                                       tGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGaa                                                                                                                                                  <  4:39370/110‑1 (MQ=255)
                                       tGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGaa                                                                                                                                                  >  3:39370/1‑110 (MQ=255)
                                       tGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCa                                                                                                           >  3:886541/1‑149 (MQ=255)
                                       tGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCa                                                                                                           >  1:32185/1‑149 (MQ=255)
                                       tGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCa                                                                                                           >  4:448306/1‑149 (MQ=255)
                                         aaaCAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCTGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCAcc                                                                                                         >  2:255954/1‑149 (MQ=255)
                                         aaaCAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCAcc                                                                                                         >  4:884436/1‑149 (MQ=255)
                                              acGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctggt                                                                                                                                          >  3:419389/1‑107 (MQ=255)
                                              acGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctggt                                                                                                                                          <  4:419389/111‑5 (MQ=255)
                                              acGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTc                                                                                                    <  4:886541/149‑1 (MQ=255)
                                                gCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCgg                                                                                                  >  4:392195/1‑149 (MQ=255)
                                                  ccAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTg                                                                                                >  2:809123/1‑149 (MQ=255)
                                                   cAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCAGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTgg                                                                                               <  4:609998/149‑1 (MQ=255)
                                                         gTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGa                                                                                         >  2:349792/1‑149 (MQ=255)
                                                         gTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGa                                                                                         <  1:827696/149‑1 (MQ=255)
                                                          ttGTCGAACCATGAGTCGATTATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGAt                                                                                        >  2:654543/1‑149 (MQ=255)
                                                          ttGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGAt                                                                                        <  4:422239/149‑1 (MQ=255)
                                                          ttGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGAt                                                                                        >  3:282108/1‑149 (MQ=255)
                                                                aaCCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcgac                                                                                  >  4:259025/1‑149 (MQ=255)
                                                                  ccATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcgacga                                                                                <  3:392195/149‑1 (MQ=255)
                                                                    aTGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGata                                                                              >  4:56578/1‑149 (MQ=255)
                                                                    aTGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGata                                                                              <  2:32185/149‑1 (MQ=255)
                                                                     tGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGatat                                                                             <  1:17103/149‑1 (MQ=255)
                                                                      gAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATc                                                                            <  2:505824/149‑1 (MQ=255)
                                                                       aGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCt                                                                           <  2:600969/149‑1 (MQ=255)
                                                                         tCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTc                                                                                                                                              <  4:270140/80‑1 (MQ=255)
                                                                         tCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTc                                                                                                                                              >  3:270140/1‑80 (MQ=255)
                                                                           gATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGc                                                                                      <  1:377182/134‑1 (MQ=255)
                                                                           gATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGc                                                                                      >  2:377182/1‑134 (MQ=255)
                                                                           gATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCaa                                                                       <  1:489897/149‑1 (MQ=255)
                                                                           gATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCaa                                                                       >  2:2952/1‑149 (MQ=255)
                                                                            aTAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAAt                                                                      <  3:48548/149‑1 (MQ=255)
                                                                                ttAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATttct                                                                  <  2:318953/149‑1 (MQ=255)
                                                                                  aGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATttcttc                                                                <  1:474591/149‑1 (MQ=255)
                                                                                      ccTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGcc                                                            >  2:57453/1‑149 (MQ=255)
                                                                                       cTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCa                                                           >  3:750766/1‑149 (MQ=255)
                                                                                       cTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCa                                                           >  1:757376/1‑149 (MQ=255)
                                                                                            aCGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGt                                                      >  2:468202/1‑149 (MQ=255)
                                                                                            aCGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGt                                                      >  3:293544/1‑149 (MQ=255)
                                                                                              gggCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTc                                                    <  1:468202/149‑1 (MQ=255)
                                                                                               ggCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTc                                                                                                                                              <  3:504276/58‑1 (MQ=255)
                                                                                               ggCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTc                                                                                                                                              >  4:504276/1‑58 (MQ=255)
                                                                                                 ccTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAgg                                                 <  2:757376/149‑1 (MQ=255)
                                                                                                  cTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGa                                                >  3:813987/1‑149 (MQ=255)
                                                                                                   ttCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTc                                                                                                    <  2:15514/96‑1 (MQ=255)
                                                                                                   ttCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTc                                                                                                    >  1:15514/1‑96 (MQ=255)
                                                                                                   ttCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGAt                                               >  4:804985/1‑149 (MQ=255)
                                                                                                    tCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATc                                              <  4:263373/149‑1 (MQ=255)
                                                                                                         aTCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTg                                         >  4:798128/1‑149 (MQ=255)
                                                                                                           cGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGcc                                       <  1:809123/149‑1 (MQ=255)
                                                                                                              cTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGc                                    >  4:286457/1‑149 (MQ=255)
                                                                                                              cTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGc                                    <  1:874681/149‑1 (MQ=255)
                                                                                                               ttCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCa                                   >  4:495938/1‑149 (MQ=255)
                                                                                                               ttCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCa                                   >  2:214831/1‑149 (MQ=255)
                                                                                                                tCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAg                                  <  3:266051/149‑1 (MQ=255)
                                                                                                                  cggcggcggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGt                                <  1:57453/149‑1 (MQ=255)
                                                                                                                   ggcggcggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTc                               >  2:463184/1‑149 (MQ=255)
                                                                                                                         ggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGtt                                                     <  1:415684/121‑1 (MQ=255)
                                                                                                                         ggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGtt                                                     >  2:415684/1‑121 (MQ=255)
                                                                                                                           aaTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGt                       <  2:700864/149‑1 (MQ=255)
                                                                                                                            aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  1:100913/149‑1 (MQ=255)
                                                                                                                            aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  4:783336/149‑1 (MQ=255)
                                                                                                                                 gcgcACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTTGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTa                 >  2:600173/1‑149 (MQ=255)
                                                                                                                                   gcACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTAcc               >  2:657988/1‑149 (MQ=255)
                                                                                                                                     aCCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      >  1:611755/1‑140 (MQ=255)
                                                                                                                                     aCCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  2:611755/140‑1 (MQ=255)
                                                                                                                                       cAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGca           <  1:121827/149‑1 (MQ=255)
                                                                                                                                       aggACGTTCTCGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGca           >  4:260948/3‑149 (MQ=255)
                                                                                                                                        ggACGTTCTCGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGcac          >  4:12790/2‑149 (MQ=255)
                                                                                                                                        aGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGcac          <  2:109609/149‑1 (MQ=255)
                                                                                                                                        aGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGAGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGcac          >  3:613096/1‑149 (MQ=255)
                                                                                                                                         gACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGcaca         >  3:110623/1‑149 (MQ=255)
                                                                                                                                            gTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTTAGTACCGGCACAACt      >  1:858583/1‑149 (MQ=255)
                                                                                                                                             ttCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACtt     >  3:343925/1‑149 (MQ=255)
                                                                                                                                               cGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCg   <  2:475933/149‑1 (MQ=255)
                                                                                                                                               cGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCg   <  1:214831/149‑1 (MQ=255)
                                                                                                                                                gagaACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCga  <  1:654543/149‑1 (MQ=255)
                                                                                                                                                gagaACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCga  >  3:100843/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCGA  >  NC_000913/2706462‑2706754

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