Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 551,570 T→G 7.6% E342A (GAA→GCA)  purK ← 5‑(carboxyamino)imidazole ribonucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913551,5700TG 7.6% 168.7 / 4.1 70E342A (GAA→GCA) purK5‑(carboxyamino)imidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base T (28/34);  new base G (2/4);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTC  >  NC_000913/551433‑551710
                                                                                                                                         |                                                                                                                                             
gatgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCcggc                                                                                                                                          >  1:558267/1‑143 (MQ=255)
 atgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGAGCGGTAGGCCGGATAAGGCGTTTACGCGGCATCCGGCAAGGATAGAGCACGAGTTAACCGAACTTCCTGTGCTCTCAACCCACGCAGCTGTGATATTGccgcg                                                                                                                                         >  2:9845/1‑143 (MQ=25)
  ttcTTGACGCTTCTTTGCAGGCCTACACTGGGTACCGGATCGGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGCGCACCCGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCcggcgg                                                                                                                                        <  2:353018/141‑1 (MQ=38)
  tgcTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCcggcgg                                                                                                                                        <  1:560104/143‑1 (MQ=255)
    cTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGgca                                                                                                                                      <  2:841112/143‑1 (MQ=255)
    cTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCCGGAATATGCCagcggcg                                                                                                                                      <  2:124611/143‑2 (MQ=38)
    cTTGACGCATCTTATCAGGCCTACAATGGGGACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCCGGAATATGCCagcggcg                                                                                                                                      <  2:124610/143‑2 (MQ=38)
      tGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAAACACGCCGCTGGCATATTCCGGCGgcagc                                                                                                                                    >  2:684963/1‑143 (MQ=255)
      tGACGCAGCTTAGCAGGCCTACAATGGGTACCGGAGCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCCGGAATATGCCagcggcgtg                                                                                                                                    <  2:813908/143‑4 (MQ=38)
           cATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGAATAGAGCACCAGGTAACCGAACTTACTCTGCGCCCAAATCACGAAGATGGCATCTCCCGACGTGAGCA‑GCgg                                                                                                                              >  2:1025456/1‑143 (MQ=37)
              cTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCagcc‑gccgtct                                                                                                                           >  1:739508/1‑138 (MQ=255)
                 agcAGGGCTGCAAGGGGTACCGGATCGGTAGGCCGAATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTaa                                                                                                                        <  2:952381/141‑1 (MQ=255)
                   cAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTACCCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAgg                                                                                                                      <  1:575833/143‑1 (MQ=255)
                     ggCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGTCCTCA‑ACGGGATTAAGGCt                                                                                                                    >  2:971970/1‑143 (MQ=255)
                      gCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGGTAACCGAACTGACTCTGCGCCCAAGTGACGCCGCTGGCATATTGCGGCGGCAGCA‑GAGGGGTTAAggatg                                                                                                                   >  2:31489/1‑140 (MQ=38)
                       ccTACAAGGGGTCCCGGAGCGGTAGGCCGGATAAGGCGTTTACGCCGCATGCTGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCTCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTc                                                                                                                  <  2:165019/143‑1 (MQ=255)
                        cTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGCACACCAGATAACCGAACTTACTCTGCGCCCAAATCACGCCACTGGCATATTCCGGCGGCCGCA‑GCGGGATCAAGGCTTcc                                                                                                                 >  2:705194/1‑143 (MQ=38)
                        cTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTACCCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATCCCGCCGGCAGCA‑GCGGGATCAAGGCTAcc                                                                                                                 >  2:199931/1‑143 (MQ=39)
                          aCAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAg                                                                                                               <  2:146138/143‑1 (MQ=255)
                            aaTGGGGACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCAGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCg                                                                                                             <  1:1051890/143‑1 (MQ=255)
                              tGGGTACCTGACCGCTAGTCCGGATACGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCTCCCATCTCACGCCGCTGGCATATTCCGGCGGCATCA‑GCGGGATTAAGGCTTCCAGCGTc                                                                                                           <  1:99236/143‑1 (MQ=255)
                              tGGGTACCGGATCGGGAGGCCGGATAAGGTGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GTGGGATTAAGGCTTCCAGCGTc                                                                                                           <  2:676151/143‑1 (MQ=255)
                              tGGGCACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTc                                                                                                           <  1:551660/143‑1 (MQ=255)
                               gggTACCGGATCGGTAGGCCGGATAAGGCGTTGACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCg                                                                                                          <  2:200370/143‑1 (MQ=255)
                                 gTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTACGGCTTCCAGCGTCGCa                                                                                                        >  1:96437/1‑143 (MQ=255)
                                 gTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCa                                                                                                        <  2:376769/143‑1 (MQ=255)
                                        aTCGGTAGAGCGGATAACGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGAc                                                                                                 <  2:410787/143‑1 (MQ=255)
                                         tCGGGAGGCCGGATCAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGAcg                                                                                                <  2:631797/143‑1 (MQ=255)
                                             taggttggATGAGGTGTTTGCGCCGCTGCCTGCAAGAATATAGCTCCTGTTTACCGACAATGCCGGGGGTCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGAt                                                                                            <  2:797654/137‑1 (MQ=255)
                                                gCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGAtgtg                                                                                         <  1:324484/143‑1 (MQ=255)
                                                         ggCGTTTACGCCGCATCCGGCAGGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAGTCACGCGGCTCGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCg                                                                                <  1:148183/143‑1 (MQ=255)
                                                         ggCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCg                                                                                <  2:602227/143‑1 (MQ=255)
                                                           cGTTTCCGACACATCCTTCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGt                                                                              <  2:682821/143‑1 (MQ=255)
                                                              ttACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCaa                                                                           >  2:42449/1‑143 (MQ=255)
                                                               tACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAGTCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCCGACGCGCTGTGTCGCTGTCGGTCaca                                                                          >  1:812075/1‑141 (MQ=255)
                                                                 cgccgtAGCCGGCAAGACTAGAGCACCAGTTGACCGGACGTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCGGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAAtt                                                                        <  2:316786/143‑1 (MQ=255)
                                                                   ccgcATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCa                                                                      <  1:469504/143‑1 (MQ=255)
                                                                   ccgcATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCa                                                                      <  2:244372/143‑1 (MQ=255)
                                                                   ccgcATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCa                                                                      >  2:258142/1‑143 (MQ=255)
                                                                    cgcATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAg                                                                     >  2:458983/1‑143 (MQ=255)
                                                                       aTCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGCTGTGTCGCTGTCGGTCAAAGTCCGacg                                                                  >  2:673694/1‑141 (MQ=255)
                                                                            gCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGCCAAATTCAGATGCCCCa                                                             >  1:98342/1‑143 (MQ=255)
                                                                              aaGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTGCGGCGGCAGCA‑GCGGGATTAAGTCTTCGATCGTCGCGGGCCGACGGGATGGGTCCCTGTCGGTCAGCTTCTGATGCCCCACt                                                           >  1:1016094/1‑143 (MQ=255)
                                                                                 aaTAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCACATTCAGATGCCCCACTTTa                                                        >  2:352676/1‑143 (MQ=255)
                                                                                          ccAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGCCAAATGCAGTTGCCCCCCTCTTCGCCCcgga                                               >  2:227025/1‑143 (MQ=255)
                                                                                           cAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAAATCAGATGCCCCACTTTACGCCCcggac                                              >  1:87306/1‑143 (MQ=255)
                                                                                              tgaaCCGAACTTCCTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCcggacgga                                           <  1:7507/141‑1 (MQ=255)
                                                                                                aaCCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGAct                                         >  2:12728/1‑143 (MQ=255)
                                                                                                   cGAACTTACTCTGCGCCCAAATCACGCCGCTGGTATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGCTGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCGCCGGACGGAcgtct                                      >  1:343856/1‑143 (MQ=255)
                                                                                                   cGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGActtct                                      >  2:338497/1‑143 (MQ=255)
                                                                                                    gAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTGTACGCCCCGGACGGActtctt                                     >  1:717459/1‑143 (MQ=255)
                                                                                                        ttACTCTGCGCCCAAGTCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTc                                 <  1:66775/143‑1 (MQ=255)
                                                                                                        ttACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTc                                 <  2:96437/143‑1 (MQ=255)
                                                                                                            tctGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTAc                             <  1:707936/143‑1 (MQ=255)
                                                                                                            ggtGCGCCCAACGCCCGCGGCTGGCATATTCCGGCGGCAGCA‑GGGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTAc                             <  1:428881/141‑1 (MQ=255)
                                                                                                              tGCGCCCAAATCACGCCGCTGGCATATGCCACCAACAGCA‑GCTCATTTCAGGCACCCAGCGTCGCAGTCAGACGCGAGGTGTCTCTTTCGTTCATACTCAGATGCCCCCCTTTACGCCCAGGACGGCATTCTTTGTCGTACCa                           >  2:703895/1‑143 (MQ=255)
                                                                                                              cgcgcCCAGAGCACGCCGCTGGCGTATGCCCG‑TGCGCCAGGGGGGATTAAGGCTTCCGGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCa                           <  2:465606/142‑1 (MQ=255)
                                                                                                                cgcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGt                         <  1:465276/143‑1 (MQ=255)
                                                                                                                cgcCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCGTTGTCGTACCAGt                         >  1:809064/1‑143 (MQ=255)
                                                                                                                cgcCCAAATCACGCCGCGGGCATATTCCGGCGGCAGCA‑GCGAGGATCAGCCTTCCAGCGTCGCAGGCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGt                         <  1:197935/143‑1 (MQ=255)
                                                                                                                   ccAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATCCCCCACTTTACCCCCCGGCCCCCCTTCTTTCTCCTACCAGTGCa                      >  1:810670/1‑143 (MQ=255)
                                                                                                                      gctgctgccgccGGAATATGCCAGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGAt                   >  1:540927/7‑143 (MQ=255)
                                                                                                                      aaTCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGCCTTCTTTGTCGTACCAGTGCAGAt                   >  1:266228/1‑143 (MQ=255)
                                                                                                                      aaTCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGAt                   <  1:971970/143‑1 (MQ=255)
                                                                                                                      aaTCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGAt                   >  1:155731/1‑143 (MQ=255)
                                                                                                                        tCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATATAGGCGTGCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGc                 <  2:749616/143‑1 (MQ=255)
                                                                                                                                cTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCgg         <  1:338497/143‑1 (MQ=255)
                                                                                                                                   gCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAg      <  1:501101/143‑1 (MQ=255)
                                                                                                                                     atatTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTAAAATTCAGATGCCCCACTTTACGCCCCGGCCGGACTTCTTTGTCGTACCAGGGCAGATGCACCAGCGCCAGtt    >  2:823615/1‑143 (MQ=255)
                                                                                                                                       atTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTc  >  1:91675/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                             
GATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGCAGCA‑GCGGGATTAAGGCTTCCAGCGTCGCAGTCAGACGCGATGTGTCGCTGTCGGTCAAATTCAGATGCCCCACTTTACGCCCCGGACGGACTTCTTTGTCGTACCAGTGCAGATGCACCAGCGGCAGTTTC  >  NC_000913/551433‑551710

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑