Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,436,195 T→A 6.5% W44R (TGG→AGG)  gpI → tail formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,436,1950TA 6.5% 183.8 / 2.6 61W44R (TGG→AGG) gpItail formation protein
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base A (2/2);  total (30/31)
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.64e-01

ACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTC  >  NC_002947/3436049‑3436341
                                                                                                                                                   |                                                                                                                                                   
aCGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGCGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCct                                                                                                                                                     <  1:14810/147‑1 (MQ=255)
aCGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCctgg                                                                                                                                                   <  6:163954/149‑1 (MQ=255)
     gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCctggctgg                                                                                                                                               <  4:11314/148‑1 (MQ=255)
     gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCtggctggc                                                                                                                                              <  3:298356/149‑1 (MQ=255)
               aCATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCAGCCCTCGACCTTCCATGCCCAGGCCAGCCAGGGCagga                                                                                                                                    <  5:319179/149‑4 (MQ=37)
                 aTGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATgg                                                                                                                                    >  7:132172/1‑147 (MQ=255)
                         gTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCctggctgg                                                                                                                                               <  8:319054/128‑1 (MQ=255)
                         gTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCctggctgg                                                                                                                                               >  7:319054/1‑128 (MQ=255)
                            aCCGCACAACACCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCt                                                                                                                       >  6:355225/1‑149 (MQ=255)
                              cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                     <  2:256691/149‑1 (MQ=255)
                              cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTCGGCATGGAAGGTCGAGGGCTgg                                                                                                                     >  7:152688/1‑149 (MQ=255)
                                   aaCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                                >  8:92122/1‑149 (MQ=255)
                                   aaCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggc                                                                                                                <  7:201739/149‑1 (MQ=255)
                                      gCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                             <  2:331140/149‑1 (MQ=255)
                                        cACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCt                                                                                                            >  6:92390/1‑148 (MQ=255)
                                          cccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTa                                                                                                          >  6:197777/1‑148 (MQ=255)
                                           ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                        <  8:132172/149‑1 (MQ=255)
                                           ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                        <  5:92390/149‑1 (MQ=255)
                                           ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                        >  1:338598/1‑149 (MQ=255)
                                                  gAGCGCGCCCTGGAGGTGGCCGACGACCTTGGGCTTGAGCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGGGCTGGCTGGGGTGGGCATGGAAGGTCGTGGGCTGGGAGGCGGCt                                                                                                            <  8:184519/138‑1 (MQ=255)
                                                  gAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATACGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCAAGGAAGGTCGAGGGCTGGGAGGCGGCt                                                                                                            >  7:184519/1‑138 (MQ=255)
                                                    gcgcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGAcc                                                                                               >  2:133000/1‑149 (MQ=255)
                                                          ccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgcc                                                                                         <  5:355225/149‑1 (MQ=255)
                                                                 gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCt                                                                                  >  6:70162/1‑149 (MQ=255)
                                                                    gccgccGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTg                                                                                                                                           <  4:289747/89‑1 (MQ=255)
                                                                    gccgccGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTg                                                                                                                                           >  3:289747/1‑89 (MQ=255)
                                                                      cgccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATcc                                                                             >  6:121314/1‑149 (MQ=255)
                                                                         cgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                           >  6:232090/1‑148 (MQ=255)
                                                                          gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGa                                                                         <  1:136321/149‑1 (MQ=255)
                                                                           aCCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                                        <  7:92122/149‑1 (MQ=255)
                                                                                  ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTc                                                                 >  5:246400/1‑149 (MQ=255)
                                                                                      gACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                               <  5:232090/147‑1 (MQ=255)
                                                                                       aCCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                            <  1:133000/149‑1 (MQ=255)
                                                                                       aCCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGt                                                             <  5:70162/148‑1 (MQ=255)
                                                                                          cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                         <  5:197777/149‑1 (MQ=255)
                                                                                          cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                         <  6:246400/149‑1 (MQ=255)
                                                                                          cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                         <  5:196016/149‑1 (MQ=255)
                                                                                          cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                          >  6:337837/1‑148 (MQ=255)
                                                                                               aTTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                                    >  8:49438/1‑149 (MQ=255)
                                                                                               aTTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                                    >  7:105605/1‑149 (MQ=255)
                                                                                                   tCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                                   <  1:82390/114‑1 (MQ=255)
                                                                                                   tCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                                   >  2:82390/1‑114 (MQ=255)
                                                                                                          cGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACgg                                         <  5:121314/149‑1 (MQ=255)
                                                                                                              gCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCg                                      <  7:49438/148‑1 (MQ=255)
                                                                                                                 gACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcgc                                  <  8:105605/149‑1 (MQ=255)
                                                                                                                       tCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGTGCGATTcg                           >  7:163851/1‑149 (MQ=255)
                                                                                                                       tCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcg                           >  1:35385/1‑149 (MQ=255)
                                                                                                                         aCGCTGCCCGCCGAACTTCCTGCCATGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAtt                                                                  >  5:97226/1‑109 (MQ=255)
                                                                                                                         aCGCTGCCCGCCGAACTTCCTGCCATGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAtt                                                                  <  6:97226/109‑1 (MQ=255)
                                                                                                                           gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                                        <  5:356191/101‑1 (MQ=255)
                                                                                                                           gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                                        >  6:356191/1‑101 (MQ=255)
                                                                                                                           gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                       >  2:315238/1‑149 (MQ=255)
                                                                                                                             tGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                       <  4:131811/147‑1 (MQ=255)
                                                                                                                             tGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                       >  3:131811/1‑147 (MQ=255)
                                                                                                                                 cctcgACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                <  6:345476/132‑1 (MQ=25)
                                                                                                                                 cctcgACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                >  5:345476/4‑135 (MQ=25)
                                                                                                                                  ctcgACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACAAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc                 >  7:252213/3‑148 (MQ=37)
                                                                                                                                     gAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAgg             <  2:35385/149‑1 (MQ=255)
                                                                                                                                      aaCTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGc            <  1:315238/149‑1 (MQ=255)
                                                                                                                                         ttCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGt                     >  4:311584/1‑137 (MQ=255)
                                                                                                                                         ttCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGt                     <  3:311584/137‑1 (MQ=255)
                                                                                                                                              tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg     >  2:185798/1‑149 (MQ=255)
                                                                                                                                              tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg     >  3:247390/1‑149 (MQ=255)
                                                                                                                                              tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCGTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAGGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg     <  2:338598/149‑1 (MQ=255)
                                                                                                                                                cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGACTTCGGCGGGTGCTCAAGGCGGTTCgg     <  1:185798/148‑1 (MQ=255)
                                                                                                                                                cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGt   >  6:43673/1‑149 (MQ=255)
                                                                                                                                                cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGATTCgg     <  5:337837/147‑1 (MQ=255)
                                                                                                                                                  ctggctggCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTc  >  3:265477/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTC  >  NC_002947/3436049‑3436341

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: