Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,706,606 G→T 5.4% L173L (CTC→CTA) ‡ 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,6060GT 5.4% 576.0 / 15.4 162L173L (CTC→CTA) ‡lepAback‑translocating elongation factor EF4, GTPase
Reads supporting (aligned to +/- strand):  ref base G (71/82);  new base T (6/3);  total (77/85)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 174 total aligned reads displayed.

CTTTGTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCGA  >  NC_000913/2706458‑2706754
                                                                                                                                                    |                                                                                                                                                    
cTTTGTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCg                                                                                                                                                      <  1:88414/149‑1 (MQ=255)
cTTTGTCGCCCTTACGCAGGGTGCCGTTTTTAAAACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCg                                                                                                                                                      <  3:1216399/149‑1 (MQ=255)
 tttGTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCga                                                                                                                                                     >  3:475682/1‑149 (MQ=255)
 tttGTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCga                                                                                                                                                     >  1:548827/1‑149 (MQ=255)
     tCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAAc                                                                                                                                                 <  1:1536547/149‑1 (MQ=255)
       gCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGt                                                                                                                                               >  2:629674/1‑149 (MQ=255)
       gCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGt                                                                                                                                               >  2:1651608/1‑149 (MQ=255)
       gCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGt                                                                                                                                               <  3:1475504/149‑1 (MQ=255)
        cccTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTc                                                                                                                                              <  3:406691/149‑1 (MQ=255)
         ccTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTcc                                                                                                                                             >  2:1571438/1‑149 (MQ=255)
         ccTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTcc                                                                                                                                             <  2:1293381/149‑1 (MQ=255)
          cTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCt                                                                                                                                            <  3:1591918/149‑1 (MQ=255)
           ttACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTg                                                                                                                                           >  2:134507/1‑149 (MQ=255)
             aCGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGca                                                                                                                                         >  1:719927/1‑149 (MQ=255)
                 aGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACAcc                                                                                                                                     <  1:71178/149‑1 (MQ=255)
                       ccGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGcc                                                                                                                               <  4:535172/149‑1 (MQ=255)
                       ccGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGcc                                                                                                                               >  2:1707270/1‑149 (MQ=255)
                          tttttAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctg                                                                                                                                            >  4:885411/1‑131 (MQ=255)
                          tttttAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGt                                                                                                                            >  1:339353/1‑149 (MQ=255)
                           ttttAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGtt                                                                                                                           <  1:144311/149‑1 (MQ=255)
                           ttttAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGtt                                                                                                                           <  1:614984/149‑1 (MQ=255)
                                aTACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTc                                                                                                                                              >  3:522928/1‑125 (MQ=255)
                                aTACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTc                                                                                                                                              <  4:522928/125‑1 (MQ=255)
                                aTACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGc                                                                                                                      <  4:345217/149‑1 (MQ=255)
                                    ggATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGaa                                                                                                                  <  4:307575/149‑1 (MQ=255)
                                     gATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAAc                                                                                                                 <  2:533721/149‑1 (MQ=255)
                                       tAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAg                                                                                                               <  4:1251695/149‑1 (MQ=255)
                                       tAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAg                                                                                                               <  3:1217633/149‑1 (MQ=255)
                                           tGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTc                                                                                                                        <  3:324838/136‑1 (MQ=255)
                                            gAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCAc                                                                                                          >  2:990837/1‑149 (MQ=255)
                                             aaaCAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCTCGAACGTctggt                                                                                                                                          <  2:143452/116‑5 (MQ=255)
                                                 aacGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGt                                                                                                     <  2:568551/149‑1 (MQ=255)
                                                  acGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTc                                                                                                    <  2:1298230/149‑1 (MQ=255)
                                                  acGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTc                                                                                                    <  2:684662/149‑1 (MQ=255)
                                                   cGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCg                                                                                                   <  4:444368/149‑1 (MQ=255)
                                                      ccAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTg                                                                                                >  1:391294/1‑149 (MQ=255)
                                                      ccAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTg                                                                                                <  2:548827/149‑1 (MQ=255)
                                                      ccAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTg                                                                                                <  4:1277053/149‑1 (MQ=255)
                                                             gTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCt                                                                                                                     >  4:1102231/1‑121 (MQ=255)
                                                             gTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCt                                                                                                                     <  3:1102231/121‑1 (MQ=255)
                                                              ttGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGc                                                                                                                                <  1:129535/109‑1 (MQ=255)
                                                                gTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGc                                                                                      >  3:1670199/1‑149 (MQ=255)
                                                                   gAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcga                                                                                   >  1:372122/1‑149 (MQ=255)
                                                                    aaCCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcgac                                                                                  >  3:785638/1‑149 (MQ=255)
                                                                      ccATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcgccga                                                                                <  4:785638/149‑1 (MQ=255)
                                                                      ccATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcgacga                                                                                <  2:20673/149‑1 (MQ=255)
                                                                       cATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGat                                                                               <  1:1396572/149‑1 (MQ=255)
                                                                         tGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCcgac                                                                                  >  3:220870/1‑144 (MQ=255)
                                                                          gAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATc                                                                            >  1:1362040/1‑149 (MQ=255)
                                                                           aGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCt                                                                           <  2:353020/149‑1 (MQ=255)
                                                                            gTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCtt                                                                          <  4:475682/149‑1 (MQ=255)
                                                                               gATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTcc                                                                                                                                             <  3:345568/79‑1 (MQ=255)
                                                                               gATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTcc                                                                                                                                             >  4:345568/1‑79 (MQ=255)
                                                                               gATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCaa                                                                       <  1:1135824/149‑1 (MQ=255)
                                                                               gATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCaa                                                                       >  4:497922/1‑149 (MQ=255)
                                                                                aTAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGAt                                                                                        <  2:1158189/131‑1 (MQ=255)
                                                                                aTAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAAt                                                                      >  1:681524/1‑149 (MQ=255)
                                                                                aTAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAAt                                                                      <  4:1670199/149‑1 (MQ=255)
                                                                                  aaTTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATg                                                                                       >  3:1513625/1‑130 (MQ=255)
                                                                                    ttAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATttct                                                                  >  2:1631862/1‑149 (MQ=255)
                                                                                      aGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGt                                                                                            <  3:1404102/121‑1 (MQ=255)
                                                                                      aGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGt                                                                                            >  4:1404102/1‑121 (MQ=255)
                                                                                      aGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATttcttc                                                                <  2:181166/149‑1 (MQ=255)
                                                                                       gTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTcc                                                               >  1:1211952/1‑149 (MQ=255)
                                                                                        tGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCga                                                                                                                                                     >  1:1042995/1‑62 (MQ=255)
                                                                                        tGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCg                                                                                                                              <  3:456061/85‑1 (MQ=255)
                                                                                           cTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCa                                                           >  2:1117111/1‑149 (MQ=255)
                                                                                              cAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTACGCCAcgc                                                        <  4:143813/149‑1 (MQ=255)
                                                                                                   ggCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCa                                                   <  1:629674/149‑1 (MQ=255)
                                                                                                     ccTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAgg                                                 <  2:1362040/149‑1 (MQ=255)
                                                                                                      cTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGa                                                <  1:990837/149‑1 (MQ=255)
                                                                                                       ttCCGGATCGCCTTCTGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGAt                                               >  4:1065345/1‑149 (MQ=255)
                                                                                                       ttCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAAt                                                                      <  1:694766/126‑1 (MQ=255)
                                                                                                       ttCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGAt                                               >  3:930885/1‑149 (MQ=255)
                                                                                                            gATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCt                                          <  1:1651608/149‑1 (MQ=255)
                                                                                                              tCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGc                                        >  2:108613/1‑149 (MQ=255)
                                                                                                                   ttCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCa                                   >  3:1445939/1‑149 (MQ=255)
                                                                                                                    tCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAg                                  <  3:368464/149‑1 (MQ=255)
                                                                                                                    tCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAg                                  <  4:1598903/149‑1 (MQ=255)
                                                                                                                    gggcgcggcggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAg                                  >  3:679699/5‑149 (MQ=255)
                                                                                                                       ggcggcggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAg                                  >  1:19720/1‑146 (MQ=255)
                                                                                                                       ggcggcggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTc                               >  1:1591207/1‑149 (MQ=255)
                                                                                                                       ggcggcggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTc                               >  4:280529/1‑149 (MQ=255)
                                                                                                                          ggcggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAAt                            >  4:591823/1‑149 (MQ=255)
                                                                                                                             ggAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCtt                         >  1:1647137/1‑149 (MQ=255)
                                                                                                                              gAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTg                        >  3:473535/1‑149 (MQ=255)
                                                                                                                                aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCATGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  4:1595220/149‑1 (MQ=255)
                                                                                                                                aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  2:717215/149‑1 (MQ=255)
                                                                                                                                aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  4:1506598/149‑1 (MQ=255)
                                                                                                                                aTGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  3:52722/149‑1 (MQ=255)
                                                                                                                                  gTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCa                    <  3:730532/149‑1 (MQ=255)
                                                                                                                                    cgcgCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGt                  <  4:901977/149‑1 (MQ=255)
                                                                                                                                     gcgcACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTa                 >  4:740300/1‑149 (MQ=255)
                                                                                                                                        tgcaggACGTTCTCGAACGTCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCgg             >  2:1375312/6‑149 (MQ=255)
                                                                                                                                          ccAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  3:566815/139‑1 (MQ=255)
                                                                                                                                          ccAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGc            >  1:717444/1‑149 (MQ=255)
                                                                                                                                          ccAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGc            <  1:1452587/149‑1 (MQ=255)
                                                                                                                                                gTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCCGGATCGGCTGCCGGCAGGTCCATCTTGtt                      >  1:1001496/1‑133 (MQ=255)
                                                                                                                                                gTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGtt                      <  2:1001496/133‑1 (MQ=255)
                                                                                                                                                  tCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTc    <  2:207163/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTTGTCGCCCTTACGCAGGGTGCCGTTTTTAATACGGATAAGTGAAACAACGCCCAGGTAGTTGTCGAACCATGAGTCGATAATTAGTGCCTGCAACGGGCCTTCCGGATCGCCTTCCGGCGGCGGAATGTCGCGCACCAGACGTTCGAGAACGTCCTGCACACCAACGCCGGTTTTCGCTGAACAGCGCACCGCGTCGGTGGCGTCGATGCCGACGATATCTTCAATTTCTTCCGCCACGCGTTCAGGATCGGCTGCCGGCAGGTCAATCTTGTTCAGTACCGGCACAACTTCGA  >  NC_000913/2706458‑2706754

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