Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 711,493 G→T 5.1% E431* (GAG→TAG) ‡ ileS → isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947711,4930GT 5.1% 316.8 / 4.3 97E431* (GAG→TAG) ‡ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (44/48);  new base T (3/2);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

GTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATCCCGT  >  NC_002947/711346‑711637
                                                                                                                                                    |                                                                                                                                                       
gtcgAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGa                                                                                                                                                        <  2:118038/149‑1 (MQ=255)
     gaAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGGGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑c                                                                                                                                                 >  2:142864/1‑149 (MQ=255)
      aaGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACTTGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑t                                                                                                                                             <  3:196975/149‑1 (MQ=255)
           gAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGt                                                                                                                                         >  6:86615/1‑148 (MQ=255)
            aGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGa                                                                                                                                       <  8:74408/149‑1 (MQ=255)
             gCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCCCTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCAGCGGCTCGCCG‑‑G‑‑‑TGctcggc                                                                                                                                      <  7:341127/149‑6 (MQ=255)
              cGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGc                                                                                                                                     >  6:12564/1‑149 (MQ=255)
                aaGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGt                                                                                                                                   <  2:201602/149‑1 (MQ=255)
                     ggTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCt                                                                                                                              >  5:39808/1‑149 (MQ=255)
                         cgcTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGaaa                                                                                                                          <  7:109421/149‑1 (MQ=255)
                         cgcTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGaaa                                                                                                                          <  2:273831/149‑1 (MQ=255)
                                  acacCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCg                                                                                                                   >  4:199479/1‑147 (MQ=255)
                                  acacCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCg                                                                                                                   >  8:145839/1‑147 (MQ=255)
                                  acacCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCg                                                                                                                   >  3:175623/1‑147 (MQ=255)
                                  acacCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑aa                                                                                                               <  3:15385/149‑1 (MQ=255)
                                  acacCGAAACCATCAGCCACAGCTACATGCACTGCTGGAGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCg                                                                                                                   >  3:169652/1‑147 (MQ=255)
                                    acCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGa                                                                                                             >  7:89246/1‑149 (MQ=255)
                                      cGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTGGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                           <  5:309230/149‑1 (MQ=255)
                                      cGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                           <  5:12564/149‑1 (MQ=255)
                                      cGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                           >  7:95046/1‑149 (MQ=255)
                                      cGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                           >  8:340850/1‑149 (MQ=255)
                                       gAAACCATCAGCAGACCGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGaca                                                                                                           <  1:131527/149‑1 (MQ=255)
                                             aTCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCATTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGGGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGt                                                                                                     >  5:200659/1‑148 (MQ=255)
                                                 gCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑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                                                                                                                          >  5:211369/1‑121 (MQ=255)
                                                      caGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGaaa                                                                                                                          <  6:211369/121‑1 (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                                                                                      <  1:183512/129‑1 (MQ=255)
                                                                                aGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTggg                                                                                      >  2:183512/1‑129 (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:10002/1‑149 (MQ=255)
                                                                                    cccGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCgatga                                                              >  6:226719/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                                    >  3:72053/1‑149 (MQ=255)
                                                                                                              aGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCAtc                                    <  5:86615/149‑1 (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                                                 <  3:210615/126‑1 (MQ=255)
                                                                                                                        gggCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGccc                                                 >  4:210615/1‑126 (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:39808/149‑1 (MQ=255)
                                                                                                                            aTGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCaa                      >  6:161127/1‑149 (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               <  3:199479/149‑1 (MQ=255)
                                                                                                                                   agcagcCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgt               <  1:142864/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    <  3:203908/149‑1 (MQ=255)
                                                                                                                                              aCCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATccc    <  7:218018/149‑1 (MQ=255)
                                                                                                                                                cGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgt               >  8:113903/1‑136 (MQ=255)
                                                                                                                                                cGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgt               <  7:113903/136‑1 (MQ=255)
                                                                                                                                                cGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATCCCGt  >  1:295034/1‑149 (MQ=255)
                                                                                                                                                cGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATCCCGt  >  8:263988/1‑149 (MQ=255)
                                                                                                                                                  gCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgt               >  2:105350/1‑134 (MQ=255)
                                                                                                                                                  gCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGgt               <  1:105350/134‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
GTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCG‑‑C‑‑‑TGCGTGAGCGTGCGCTGAAAGCCATCG‑‑AAGACACCAAGTTCGTCCCGGCCTGGGGCCAGGCGCGCCTGCACGCGATGATCGCCAACCGCCCGGACTGGTGCATCTCGCGTCAGCGCAACTGGGGTGTGCCGATCCCGT  >  NC_002947/711346‑711637

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: