Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,504,800 T→A 6.1% K250* (AAG→TAG) ‡ ypdE ← aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,504,8000TA 6.1% 175.3 / 2.2 66K250* (AAG→TAG) ‡ypdEaminopeptidase
Reads supporting (aligned to +/- strand):  ref base T (28/33);  new base A (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CCACTGTGGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATC  >  NC_000913/2504653‑2504946
                                                                                                                                                   |                                                                                                                                                  
ccacTGTGGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCGCGGGGGTAAGTTTGGGCGGCGCAATCAGCGACTTg                                                                                                                                                   >  2:654070/1‑149 (MQ=255)
    tgtgGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGc                                                                                                                                               <  2:679004/149‑1 (MQ=255)
      tgGGCACGCCGGTGGCGGTTAAGTGCACCGCCACGCCGTCCGTGCCGCCGTTGCTGAACTTATCTGCCTGCAACGGCACGCCTATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTTTCGCtt                                                                                                                                             <  2:634108/149‑1 (MQ=255)
      ggggCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCtt                                                                                                                                             <  2:31727/148‑1 (MQ=255)
       gggCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGATCATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACTGTTTCGACCCAGGCGGGAAGTTTTGGCGGCGCAATCAGCGACAAGTCGCtga                                                                                                                                            <  1:645366/149‑3 (MQ=255)
          cACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAaca                                                                                                                                         <  2:283497/149‑1 (MQ=255)
           aCGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAacac                                                                                                                                        <  2:397926/149‑1 (MQ=255)
                     cggtTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGc                                                                                                                              <  1:678465/149‑1 (MQ=255)
                      ggtTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGcc                                                                                                                             >  2:1028934/1‑149 (MQ=255)
                          aaGTGCACCGCCACGCCGTCAGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGACCGtt                                                                                                                         >  1:469460/1‑149 (MQ=255)
                           aGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTa                                                                                                                        <  2:422992/149‑1 (MQ=255)
                           aGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTa                                                                                                                        <  1:353028/149‑1 (MQ=255)
                               cACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCaa                                                                                                                    >  1:984098/1‑149 (MQ=255)
                                   gccccgccGTCCGTGCCGCCGTTGCTGATCATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCGCGGTTTCGACCCCGGCGGTAAGTTTTGGCGGCGCAATCAGCGCCTTGTCGCTTAACCCCAGCATCGGCCCGTTCCCAATCTg                                                                                                                >  1:662330/1‑149 (MQ=255)
                                       cgccgTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGGTTCGACCCAGGCGGGATGTTTTGGCGGCGCAATCAGCGACTGGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGa                                                                                                            >  1:941606/1‑149 (MQ=255)
                                            tccgtGCCGCCGTTGCTGAACACATCTGCCTGCAACGGCACGCCAATTTCTGCCGACACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGtt                                                                                                       <  2:93528/149‑1 (MQ=255)
                                               gtGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGc                                                                                                    >  2:1034754/1‑149 (MQ=255)
                                                tGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGcc                                                                                                   <  2:377093/149‑1 (MQ=255)
                                                  ccgccgTGGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCcgc                                                                                                 <  2:596594/149‑1 (MQ=255)
                                                   cgccgTTGCTGAACATATCTGCCTGCAACGGCACGCCAAATTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGACGCTTAACCCCAGCATCGGCCCGTTACCAATCAGGCGATGGTAAGCcgcg                                                                                                >  2:590612/1‑149 (MQ=255)
                                                    gccgTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCgcgc                                                                                               <  2:442800/149‑1 (MQ=255)
                                                           cTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATc                                                                                        >  1:1025594/1‑149 (MQ=255)
                                                            tGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCTTTACCAATCTGGCTATGGTTAGCCGCTCCATAATCa                                                                                       >  1:773735/1‑149 (MQ=255)
                                                            tGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCa                                                                                       <  2:1019984/149‑1 (MQ=255)
                                                              aaCATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGTGATGGTTAGCCGCGCCATAATCaaa                                                                                     <  2:585880/149‑1 (MQ=255)
                                                              aaCATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCaaa                                                                                     >  1:986942/1‑149 (MQ=255)
                                                                     ctgcctgcAAAGGCACGCCAATGTCCGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTc                                                                              <  2:1140900/149‑1 (MQ=255)
                                                                      tgcccgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCg                                                                             <  1:640297/149‑1 (MQ=255)
                                                                       gcctgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGttt                                                                                 >  1:516978/1‑144 (MQ=255)
                                                                       gcctgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGttt                                                                                 <  2:516978/144‑1 (MQ=255)
                                                                          tgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCca                                                                         >  1:426278/1‑149 (MQ=255)
                                                                            cAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCcagc                                                                       >  1:93773/1‑149 (MQ=255)
                                                                                 gCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGTCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCg                                                                  <  2:266515/149‑1 (MQ=255)
                                                                                  cACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGATACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCgg                                                                 >  1:289472/1‑149 (MQ=255)
                                                                                     gCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCTTCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTa                                                               <  2:355740/148‑1 (MQ=255)
                                                                                     gCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCTTCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTa                                                               >  1:355740/1‑148 (MQ=255)
                                                                                      ccAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATc                                                             <  1:413247/149‑1 (MQ=255)
                                                                                      ccAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATc                                                             >  2:722929/1‑149 (MQ=255)
                                                                                       cAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTATCACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCa                                                            >  1:331218/1‑149 (MQ=255)
                                                                                        aaTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCTCAATCAGCTAATTGTCGCTTAACACCAGCATCGGCCCGATACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCaa                                                           <  2:611448/149‑1 (MQ=255)
                                                                                        aaTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCAAAAACAAAGTTTTTCGCCCAGCAGGCGGTATCaa                                                           >  2:465442/1‑149 (MQ=255)
                                                                                         aTTTGTGCCGCCACGGTTTCGACCCAGGCGGTAAGGTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCGCGTTACCAATCTGGCGAGGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAg                                                          <  1:692268/149‑1 (MQ=255)
                                                                                         aTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAg                                                          <  2:482039/149‑1 (MQ=255)
                                                                                         aTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAg                                                          <  1:1004397/149‑1 (MQ=255)
                                                                                         aTTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAg                                                          <  2:146611/149‑1 (MQ=255)
                                                                                              tGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACaa                                                     >  2:791603/1‑149 (MQ=255)
                                                                                               gccgccACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAAt                                                    <  1:595223/149‑1 (MQ=255)
                                                                                                ccgccACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATg                                                   <  2:241104/149‑1 (MQ=255)
                                                                                                   ccACGGTTTCGCCCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGCCTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCg                                                >  1:959835/1‑149 (MQ=255)
                                                                                                     aCGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGAc                                              >  1:436915/1‑149 (MQ=255)
                                                                                                       ggTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGCCGATGGTTAGCCGCGCCATAATCAAAGGTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGt                                            >  2:566643/1‑149 (MQ=255)
                                                                                                          ttCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCg                                         >  2:440947/1‑149 (MQ=255)
                                                                                                               cccAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGac                                    <  1:917958/149‑1 (MQ=255)
                                                                                                                    gCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACcgc                               >  1:1006055/1‑149 (MQ=255)
                                                                                                                       gTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACcgcgcg                            >  1:768305/1‑149 (MQ=255)
                                                                                                                        tacgTTGTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCgg                           <  1:452522/146‑1 (MQ=255)
                                                                                                                        tAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCgg                           <  1:457756/149‑1 (MQ=255)
                                                                                                                            ttttGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACCCCGCGCGGGTgg                       >  1:303094/1‑149 (MQ=255)
                                                                                                                                  cggcgCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGttt                 <  1:791603/149‑1 (MQ=255)
                                                                                                                                   ggcgCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTATATCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTg                <  1:598130/149‑1 (MQ=255)
                                                                                                                                    gcgcAATCAGCGACTTGTCGCTTAACACCAGCATCGGCGCGTTACCAATCTGGCGCTGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCTGTTTGc               >  2:651180/1‑149 (MQ=255)
                                                                                                                                     cgcAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGcc              <  2:1119114/149‑1 (MQ=255)
                                                                                                                                      gtgttaaGCGACAAGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGccc             >  1:797750/7‑149 (MQ=255)
                                                                                                                                       tgttaaGCGACAAGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCccg            >  2:701725/6‑149 (MQ=255)
                                                                                                                                           aaGCGACAAGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCcgccgc        <  2:797750/148‑1 (MQ=255)
                                                                                                                                                 cTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCACCCGCGTAATc  >  1:109874/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCACTGTGGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATC  >  NC_000913/2504653‑2504946

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: