Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,111,153 C→G 5.9% A116P (GCG→CCG)  pfpI ← protease PfpI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,111,1530CG 5.9% 217.0 / 4.3 68A116P (GCG→CCG) pfpIprotease PfpI
Reads supporting (aligned to +/- strand):  ref base C (36/28);  new base G (2/2);  total (38/30)
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.71e-01

GTGATGAGGTTGCCATCGACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTAGAGGTCAAGCTGCGCACTGTC  >  NC_002947/3111020‑3111300
                                                                                                                                     |                                                                                                                                                    
gTGATGAGGTTGCCATCGACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAAc                                                                                                                                       >  3:210111/1‑149 (MQ=255)
 tgatgaGGTTGCCATCGACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAAcc                                                                                                                                      <  8:345258/149‑1 (MQ=255)
  gatgaGGTTGCCATCGACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCg                                                                                                                                     >  4:240512/1‑149 (MQ=255)
   atgaGGTTGCCATCGACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTTGATATCAGCAGCCACGCGCCGTGGCAGATAACCGc                                                                                                                                    >  3:26752/1‑149 (MQ=255)
        ggtGCCATCGACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCg                                                                                                                               >  5:220407/3‑149 (MQ=255)
        gTTGCCATCGACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCg                                                                                                                               >  3:327247/1‑149 (MQ=255)
                 gACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTgccgg                                                                                                                      <  6:282308/149‑1 (MQ=255)
                 gACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTgccgg                                                                                                                      >  8:182324/1‑149 (MQ=255)
                  accaccACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTgccggc                                                                                                                     >  6:54507/1‑149 (MQ=255)
                   ccaccaCTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGGCGCGTGGctg                                                                                                                                            <  4:89155/125‑3 (MQ=255)
                   ccaccaCTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGGCGCGTGGctg                                                                                                                                            >  3:105173/1‑123 (MQ=255)
                   ccaccaCTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGGCGCGTGGctg                                                                                                                                            <  4:105173/125‑3 (MQ=255)
                   ccaccaCTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGGCGCGTGGctg                                                                                                                                            >  3:89155/1‑123 (MQ=255)
                   ccaccaCTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTgccggcc                                                                                                                    >  7:194231/1‑149 (MQ=255)
                   ccaccaCTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTgccggcc                                                                                                                    <  4:210111/149‑1 (MQ=255)
                    caccacTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTgccggccg                                                                                                                   >  1:237029/1‑149 (MQ=255)
                       cacTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGca‑t                                                                                                               >  7:103783/1‑149 (MQ=255)
                            gTTAATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATACCCGCCAGCGGCCTGCCGGCCGCA‑TCATgg                                                                                                          >  4:61759/1‑149 (MQ=255)
                                aTCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCttt                                                                                                      <  7:182324/149‑1 (MQ=255)
                                 tCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCtttt                                                                                                     >  8:169840/1‑149 (MQ=255)
                                   cacccacGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGGCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTa                                                                                                   <  3:61759/149‑1 (MQ=255)
                                     cccacGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTAcc                                                                                                 <  5:178487/149‑1 (MQ=255)
                                      ccacGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTAc                                                                                                  >  2:243851/1‑147 (MQ=255)
                                               gcccgcGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTgcg                                                                                       <  3:240512/149‑1 (MQ=255)
                                                cccgcGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTg                                                                                                                                                <  8:140794/92‑1 (MQ=255)
                                                cccgcGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTg                                                                                                                                                >  7:140794/1‑92 (MQ=255)
                                                  cgcgTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCAc                                                                                    >  3:267517/1‑149 (MQ=255)
                                                    cgTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCg                                                                                  <  4:327247/149‑1 (MQ=255)
                                                         cGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGAt                                                                             >  5:325498/1‑149 (MQ=255)
                                                           aaTGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCgg                                                                                 >  7:113024/1‑143 (MQ=255)
                                                           aaTGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCgg                                                                                 <  8:113024/143‑1 (MQ=255)
                                                             tGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAg                                                                         <  5:323758/149‑1 (MQ=255)
                                                              gTCATCTTGCAAGGTCCTGTAGCTGGTCACTCGCTTGCCTTTGCCCAACCCGCTCAATATCCCCAGCCACGCGCCGGGGCAGATAACCGCCAGCGGCCCGCCGGCCGCA‑TCATGGCTCTTTACCAGTTTCTGCGCACCGGGGATCAAGc                                                                        <  6:139482/149‑1 (MQ=255)
                                                                 aTCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCAGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTAACAGTTTCTGCGCACCGGGGATCAAGCGa                                                                      >  2:177676/1‑148 (MQ=255)
                                                                       cAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGGCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCgtg                                                                <  2:237029/148‑1 (MQ=255)
                                                                         aGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCACTCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTc                                                              <  1:284390/149‑1 (MQ=255)
                                                                              ttGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGt                                                         >  6:123822/1‑148 (MQ=255)
                                                                                   gCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGAt                                                                             <  8:215593/123‑1 (MQ=255)
                                                                                   gCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGAt                                                                             >  7:215593/1‑123 (MQ=255)
                                                                                        tCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTc                                                                                           >  5:285180/1‑104 (MQ=255)
                                                                                        tCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTc                                                                                           <  6:285180/104‑1 (MQ=255)
                                                                                        tCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACAccg                                              >  2:91890/1‑149 (MQ=255)
                                                                                          aTTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACAccgcc                                            >  7:266205/1‑149 (MQ=255)
                                                                                              gTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCggg                                         <  4:267517/148‑1 (MQ=255)
                                                                                              gTTTGCCTTTGGCCAACCAGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACCCAGCCggg                                         >  4:320846/1‑148 (MQ=255)
                                                                                                    cTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGtt                                                        >  7:70458/1‑127 (MQ=255)
                                                                                                    cTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGtt                                                        <  8:70458/127‑1 (MQ=255)
                                                                                                    cTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCt                                                      <  6:129501/129‑1 (MQ=255)
                                                                                                    cTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCt                                                      >  5:129501/1‑129 (MQ=255)
                                                                                                    cTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCAc                                   >  3:214905/1‑148 (MQ=255)
                                                                                                     tttGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCAcc                                  >  2:360430/1‑148 (MQ=255)
                                                                                                      ttGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCAcc                                  >  2:33889/1‑147 (MQ=255)
                                                                                                        ggCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAg                               <  5:123822/148‑1 (MQ=255)
                                                                                                         gCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGc                             <  1:177676/149‑1 (MQ=255)
                                                                                                           cAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTgc                                                                                        <  4:106341/88‑1 (MQ=255)
                                                                                                           cAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTgc                                                                                        >  3:106341/1‑88 (MQ=255)
                                                                                                           cAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGcgt                           >  4:169514/1‑149 (MQ=255)
                                                                                                           cAACCCGCTCGATATAAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCTCA‑TCATGGCTTTTTACAAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGcgt                           >  4:163805/1‑149 (MQ=255)
                                                                                                              cccGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTCCACACCGCCGGGCAGCACCAAGGcgtcg                         >  4:285580/1‑148 (MQ=255)
                                                                                                                cGCTTGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTag                      <  4:26752/149‑1 (MQ=255)
                                                                                                                            gcagcCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGcggtgtgc                      >  5:181194/1‑131 (MQ=255)
                                                                                                                            gcagcCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGcggtgtgc                      <  6:181194/137‑7 (MQ=255)
                                                                                                                            gcagcCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTag                      >  7:236118/1‑137 (MQ=255)
                                                                                                                            gcagcCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTag                      <  8:236118/137‑1 (MQ=255)
                                                                                                                            gcagcCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTAGAGGTCAAGCTgc          <  4:214905/149‑1 (MQ=255)
                                                                                                                            gcagcCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTAGAGGTCAAGCTgc          <  8:194231/149‑1 (MQ=255)
                                                                                                                            gcagcCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTAGAGGTCAAGCTgc          <  1:360430/149‑1 (MQ=255)
                                                                                                                                  aCGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTAGAGGTCAAGCTGCGCACTg    >  5:374957/1‑149 (MQ=255)
                                                                                                                                    gcgcCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTAGAGGTCAAGCTGCGCACTGTc  <  7:169840/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
GTGATGAGGTTGCCATCGACCACCACTTGTTCATCCACCCACGTACCGCCCGCGTTGCGAATGTCATCTTGCAAGGTCTTGTAGCTGGTCATTCGTTTGCCTTTGGCCAACCCGCTCGATATCAGCAGCCACGCGCCGTGGCAGATAACCGCCAGCGGCCTGCCGGCCGCA‑TCATGGCTTTTTACCAGTTTCTGCGCACCGGGGATCAAGCGAATCGTGTCGGAGTTCTGCACACCGCCGGGCAGCACCAAGGCGTCGTAGAGGTCAAGCTGCGCACTGTC  >  NC_002947/3111020‑3111300

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