Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 711,488 C→G 5.4% T429S (ACC→AGC)  ileS → isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947711,4880CG 5.4% 300.0 / 4.8 92T429S (ACC→AGC) ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (40/47);  new base G (3/2);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATCCC  >  NC_002947/711340‑711635
                                                                                                                                                     |                                                                                                                                                          
gCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCAc                                                                                                                                                            <  5:99934/149‑1 (MQ=255)
    tcgtcgAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGc                                                                                                                                                        <  5:4881/149‑1 (MQ=255)
      gtcgAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGa                                                                                                                                                      <  2:118038/149‑1 (MQ=255)
           gaAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGGGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑c                                                                                                                                               >  2:142864/1‑149 (MQ=255)
            aaGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACTTGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑t                                                                                                                                           <  3:196975/149‑1 (MQ=255)
                 gAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGt                                                                                                                                       >  6:86615/1‑148 (MQ=255)
                  aGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGa                                                                                                                                     <  8:74408/149‑1 (MQ=255)
                   gCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCCCTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCAGCGGCTCGCCG‑‑G‑‑‑TGctcggc                                                                                                                                    <  7:341127/149‑6 (MQ=255)
                    cGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGc                                                                                                                                   >  6:12564/1‑149 (MQ=255)
                      aaGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGt                                                                                                                                 <  2:201602/149‑1 (MQ=255)
                           ggTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCt                                                                                                                            >  5:39808/1‑149 (MQ=255)
                               cgcTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGaaa                                                                                                                        <  2:273831/149‑1 (MQ=255)
                               cgcTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGaaa                                                                                                                        <  7:109421/149‑1 (MQ=255)
                                        acacCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCg                                                                                                                 >  8:145839/1‑147 (MQ=255)
                                        acacCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCg                                                                                                                 >  4:199479/1‑147 (MQ=255)
                                        acacCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCg                                                                                                                 >  3:175623/1‑147 (MQ=255)
                                        acacCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑aa                                                                                                             <  3:15385/149‑1 (MQ=255)
                                        acacCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGAGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCg                                                                                                                 >  3:169652/1‑147 (MQ=255)
                                          acCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGa                                                                                                           >  7:89246/1‑149 (MQ=255)
                                            cGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTGGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                         <  5:309230/149‑1 (MQ=255)
                                            cGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                         >  8:340850/1‑149 (MQ=255)
                                            cGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                         <  5:12564/149‑1 (MQ=255)
                                            cGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                         >  7:95046/1‑149 (MQ=255)
                                             gAAACCATCAGCAGACCGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                         <  1:131527/149‑1 (MQ=255)
                                                   aTCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCATTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGGGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGt                                                                                                   >  5:200659/1‑148 (MQ=255)
                                                       gC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTc                                                                                              <  1:167792/149‑1 (MQ=255)
                                                           acaGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCg                                                                                           <  4:175838/149‑1 (MQ=255)
                                                           acaGATACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCg                                                                                           >  4:293817/1‑149 (MQ=255)
                                                            caGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCTCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCgg                                                                                          <  7:145839/149‑1 (MQ=255)
                                                            caGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGaaa                                                                                                                        <  6:211369/121‑1 (MQ=255)
                                                            caGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGaaa                                                                                                                        >  5:211369/1‑121 (MQ=255)
                                                               cTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCt                                                                                       >  4:203908/1‑149 (MQ=255)
                                                                tACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCTTCCCGGCCTg                                                                                      >  8:238925/1‑149 (MQ=255)
                                                                  cATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTggg                                                                                    >  2:36199/1‑149 (MQ=255)
                                                                  cATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTgg                                                                                     <  6:96112/148‑1 (MQ=255)
                                                                  cATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTg                                                                                      <  5:433/147‑1 (MQ=255)
                                                                         tgctgGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGgc                                                                             >  3:307633/1‑149 (MQ=255)
                                                                           ctgGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGgc                                                                             <  4:223841/147‑1 (MQ=255)
                                                                               cgcCACAAGACCCCGCTGATCTACCGCGCCACCGCGCATTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGGGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCt                                                                       <  6:200659/149‑1 (MQ=255)
                                                                                gcCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGTTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTg                                                                      >  4:170200/1‑149 (MQ=255)
                                                                                  cacaAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGTTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCa                                                                    <  3:170200/149‑1 (MQ=255)
                                                                                      aGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTggg                                                                                    >  2:183512/1‑129 (MQ=255)
                                                                                      aGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTggg                                                                                    <  1:183512/129‑1 (MQ=255)
                                                                                       gACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCga                                                               <  1:36199/149‑1 (MQ=255)
                                                                                          cccGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGGCCAGGTGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCAcgc                                                                 >  1:118731/1‑149 (MQ=255)
                                                                                          cccGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCgatga                                                            >  6:226719/1‑149 (MQ=255)
                                                                                          cccGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCgatga                                                            >  6:10002/1‑149 (MQ=255)
                                                                                               tGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTcc                                                                                             <  6:292167/111‑1 (MQ=255)
                                                                                               tGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTcc                                                                                             >  5:292167/1‑111 (MQ=255)
                                                                                               tGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGc                                                        >  6:130526/1‑148 (MQ=255)
                                                                                                   cTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTc                                                                                                 <  5:188593/103‑1 (MQ=255)
                                                                                                   cTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTc                                                                                                 >  6:188593/1‑103 (MQ=255)
                                                                                                      ccGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCaa                                                     <  7:290033/144‑1 (MQ=255)
                                                                                                       ctcgcCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGATCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGccc                                               <  5:204963/147‑1 (MQ=255)
                                                                                                               cgcgCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGt                                       <  5:10002/149‑1 (MQ=255)
                                                                                                               cgcgCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGt                                       <  4:292724/149‑1 (MQ=255)
                                                                                                                 cgcAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGc                                     <  5:130526/149‑1 (MQ=255)
                                                                                                                    aGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCAtc                                  <  5:86615/149‑1 (MQ=255)
                                                                                                                    aGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCAtc                                  >  3:72053/1‑149 (MQ=255)
                                                                                                                     gTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCAtct                                 >  7:110881/1‑149 (MQ=255)
                                                                                                                        gTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCTTCCCGGCCTGGGGCCAGGTGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTcgc                              <  7:238925/149‑1 (MQ=255)
                                                                                                                        gTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTcgc                              <  8:110881/149‑1 (MQ=255)
                                                                                                                            gTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGTCCAGGCGCGCCTGCACGCGATGATCGCCACCCGCCCGGACTGGTGCATCTCGCGTCa                          >  3:252907/1‑149 (MQ=255)
                                                                                                                             tGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAg                         >  2:302627/1‑149 (MQ=255)
                                                                                                                              gggCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGccc                                               >  4:210615/1‑126 (MQ=255)
                                                                                                                              gggCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGccc                                               <  3:210615/126‑1 (MQ=255)
                                                                                                                               ggCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCGATAAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAg                         >  1:108482/1‑149 (MQ=255)
                                                                                                                               ggCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAgcg                       <  5:226719/149‑1 (MQ=255)
                                                                                                                                gCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAgcg                       >  8:46423/1‑148 (MQ=255)
                                                                                                                                 cATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTc                                                                                              <  8:203235/76‑1 (MQ=255)
                                                                                                                                  aTGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCaa                    >  6:161127/1‑149 (MQ=255)
                                                                                                                                  aTGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCaa                    <  6:39808/149‑1 (MQ=255)
                                                                                                                                   tGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAAc                   <  3:293817/149‑1 (MQ=255)
                                                                                                                                     gACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTg                 >  6:197428/1‑149 (MQ=255)
                                                                                                                                     gACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTg                 >  3:78986/1‑149 (MQ=255)
                                                                                                                                      aCAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTgg                <  2:3448/149‑1 (MQ=255)
                                                                                                                                         agcagcCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgt             <  1:142864/149‑1 (MQ=255)
                                                                                                                                         agcagcCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgt             <  3:199479/149‑1 (MQ=255)
                                                                                                                                         agcagcCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgt             <  7:340850/149‑1 (MQ=255)
                                                                                                                                         agcagcCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgt             <  5:197428/149‑1 (MQ=255)
                                                                                                                                          gcagcCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgtg            >  3:72387/1‑149 (MQ=255)
                                                                                                                                            agcCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAgcgc                      <  3:30636/137‑1 (MQ=255)
                                                                                                                                            agcCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAgcgc                      >  4:30636/1‑137 (MQ=255)
                                                                                                                                              tcacgCAGCGGCTCGCCG‑‑G‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTggg               >  3:237437/5‑142 (MQ=255)
                                                                                                                                              tcacgCAGCGGCTCGCCG‑‑G‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTggg               <  4:237437/138‑1 (MQ=255)
                                                                                                                                              tcacgCAGCGGCTCGCCG‑‑G‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGcc        >  2:367818/5‑149 (MQ=255)
                                                                                                                                              ccGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGcc        >  8:218018/1‑149 (MQ=255)
                                                                                                                                                gAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGa      >  2:224765/1‑149 (MQ=255)
                                                                                                                                                acgCAGCGGCTCGCCG‑‑G‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGa      >  1:198321/3‑149 (MQ=255)
                                                                                                                                                  gCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATc    <  8:95046/149‑1 (MQ=255)
                                                                                                                                                  gCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATc    <  4:78986/149‑1 (MQ=255)
                                                                                                                                                    aCCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATccc  <  7:218018/149‑1 (MQ=255)
                                                                                                                                                    aCCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATccc  <  3:203908/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                          
GCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGC‑CACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATCCC  >  NC_002947/711340‑711635

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