Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,436,200 G→A 7.4% L45L (CTG→CTA) ‡ gpI → tail formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,436,2000GA 7.4% 200.6 / 3.5 67L45L (CTG→CTA) ‡gpItail formation protein
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base A (2/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCG  >  NC_002947/3436054‑3436347
                                                                                                                                                    |                                                                                                                                                   
gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCtggctggc                                                                                                                                                   >  8:321022/1‑149 (MQ=255)
gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCtggctggc                                                                                                                                                   >  3:30684/1‑149 (MQ=255)
        tgAAATGAGCCTGTTACCGCACAACGCCACCCATTTGGAGCGCGCCCTGGAGGCGGCCGCCGTCCTTGGCCTTTACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGTCCCGCCGAACTTCCTGCCCAGGCCAGCCAGGGCa                                                                                                                                            <  1:229480/149‑1 (MQ=255)
          aCATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATgg                                                                                                                                         <  7:321022/149‑1 (MQ=255)
          aCATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATgg                                                                                                                                         <  2:80552/149‑1 (MQ=255)
              gAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAg                                                                                                                                      <  1:334050/148‑1 (MQ=255)
                    ggtACCGCACAACGCCACCCAGTGGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGCCCCGGACATTAGCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAgg                                                                                                                               <  1:182392/147‑1 (MQ=255)
                     ttACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAggg                                                                                                                              >  6:189796/1‑149 (MQ=255)
                       aCCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCt                                                                                                                            >  7:327920/1‑149 (MQ=255)
                       aCCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCt                                                                                                                            >  2:269261/1‑149 (MQ=255)
                        ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTg                                                                                                                           >  6:181116/1‑149 (MQ=255)
                        ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTg                                                                                                                           <  2:184178/149‑1 (MQ=255)
                           cacaACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTggg                                                                                                                                              >  5:211678/1‑127 (MQ=255)
                           cacaACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTggg                                                                                                                                              <  6:211678/127‑1 (MQ=255)
                                cGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGTCATGGAAGGTCGAGGGCTGGGAggcgg                                                                                                                   >  6:93603/1‑149 (MQ=255)
                                cGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcgg                                                                                                                   <  3:247645/149‑1 (MQ=255)
                                 gCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                                  <  7:93351/149‑1 (MQ=255)
                                 gCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                                  >  5:346036/1‑149 (MQ=255)
                                    aCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTa                                                                                                               >  7:240956/1‑149 (MQ=255)
                                      ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGTGCTGGGAGGCGGCTTaca                                                                                                             >  8:235029/1‑149 (MQ=255)
                                       cAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAggg                                                                                                                              >  2:213787/1‑131 (MQ=255)
                                       cAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAggg                                                                                                                              <  1:213787/131‑1 (MQ=255)
                                                 gcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAg                                                                                                  >  1:339101/1‑149 (MQ=255)
                                                 gcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAg                                                                                                  >  5:185657/1‑149 (MQ=255)
                                                              ggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                               >  3:105558/1‑91 (MQ=255)
                                                              ggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                               <  4:105558/91‑1 (MQ=255)
                                                                 cgccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATcc                                                                                  >  1:192239/1‑149 (MQ=255)
                                                                 cgccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTAACACCGACGACCAGCGCCGCGCGCTGATcc                                                                                  >  2:283236/1‑149 (MQ=255)
                                                                  gccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcg                                                                                 >  4:347440/1‑149 (MQ=255)
                                                                     gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                                >  5:160145/1‑147 (MQ=255)
                                                                             ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCggcg                                                                                                                               <  2:106397/92‑4 (MQ=17)
                                                                             ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCggcg                                                                                                                               >  1:106397/1‑89 (MQ=17)
                                                                             ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTc                                                                      >  1:126526/1‑149 (MQ=255)
                                                                              cTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTcc                                                                     <  6:346036/149‑1 (MQ=255)
                                                                                   cccGGATAATGTCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                                 >  1:83437/1‑148 (MQ=255)
                                                                                    ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                               <  5:181116/149‑1 (MQ=255)
                                                                                    ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCCACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                               >  3:240550/1‑149 (MQ=255)
                                                                                     cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                              <  3:165315/149‑1 (MQ=255)
                                                                                      ggACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTCGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                              >  7:3783/1‑148 (MQ=255)
                                                                                       gACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                            <  5:189796/149‑1 (MQ=255)
                                                                                       gACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCACGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                              >  6:212971/1‑147 (MQ=255)
                                                                                                cGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                                    <  3:347440/148‑1 (MQ=255)
                                                                                                cGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                                    <  8:240956/148‑1 (MQ=255)
                                                                                                 ggggCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGACTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAggg                                                   >  4:287458/1‑148 (MQ=255)
                                                                                                  gggCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCa                                                 <  6:211347/149‑1 (MQ=255)
                                                                                                         gCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCg                                           <  2:339101/148‑1 (MQ=255)
                                                                                                         gCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCGTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCgg                                          <  4:30684/149‑1 (MQ=255)
                                                                                                           cGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGTCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                                        <  5:93603/149‑1 (MQ=255)
                                                                                                           cGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                                        >  4:167404/1‑149 (MQ=255)
                                                                                                           cGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                                        >  5:44564/1‑149 (MQ=255)
                                                                                                            gACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                               >  1:170596/1‑125 (MQ=255)
                                                                                                            gACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                               <  2:170596/125‑1 (MQ=255)
                                                                                                            gACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcgc                                       <  2:126526/149‑1 (MQ=255)
                                                                                                                  tCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcga                                                                                                         <  4:151277/77‑1 (MQ=255)
                                                                                                                  tCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcga                                                                                                         >  3:151277/1‑77 (MQ=255)
                                                                                                                    aCGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggc                               <  2:192239/149‑1 (MQ=255)
                                                                                                                      gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                                        >  2:128265/1‑138 (MQ=255)
                                                                                                                      gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                             <  6:160145/149‑1 (MQ=255)
                                                                                                                      gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                             >  4:5491/1‑149 (MQ=255)
                                                                                                                      gCTGCCCGCCGAACTTCCTGCCCTGACTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                             <  3:287458/149‑1 (MQ=255)
                                                                                                                                 gAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCTGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAgg                   <  1:269261/149‑1 (MQ=255)
                                                                                                                                 gAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAgg                   <  8:327920/149‑1 (MQ=255)
                                                                                                                                     ttCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGGGCGATTCGGCGGGTg                          >  3:238205/1‑138 (MQ=255)
                                                                                                                                     ttCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGGGCGATTCGGCGGGTg                          <  4:238205/138‑1 (MQ=255)
                                                                                                                                       ccatGCCCAGGCCAGCCAGGGCAGGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                                        <  6:185954/71‑1 (MQ=18)
                                                                                                                                       ccatGCCCAGGCCAGCCAGGGCAGGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                                        >  5:185954/4‑74 (MQ=18)
                                                                                                                                          tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgc           >  5:117469/1‑148 (MQ=255)
                                                                                                                                            cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg           <  2:83437/147‑1 (MQ=255)
                                                                                                                                            cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg           <  7:235029/147‑1 (MQ=255)
                                                                                                                                                   tggcCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGAGTCGGCGGGTGCTCAAGGCGGTTCGGGTCAAcgcg  >  2:315048/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCG  >  NC_002947/3436054‑3436347

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: