Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,989,005 C→A 7.4% A508D (GCC→GAC)  oplA → 5‑oxoprolinase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,989,0050CA 7.4% 198.3 / 3.4 63A508D (GCC→GAC) oplA5‑oxoprolinase A
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base A (2/3);  total (31/32)
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.98e-01

CGGGCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGA  >  NC_002947/3988858‑3989152
                                                                                                                                                    |                                                                                                                                                   
cggGCATGCAAGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGATGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTAGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGcc                                                                                                                                                    >  8:189516/1‑149 (MQ=255)
   gCATGCCCGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACCCCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCCGCCTGTCGATGCGCCTGGACGCCggg                                                                                                                                                 >  2:49085/1‑149 (MQ=255)
   gCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCggg                                                                                                                                                 <  8:269268/149‑1 (MQ=255)
        cACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACg                                                                                                                                             <  5:178320/148‑1 (MQ=255)
            tggtgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGAc                                                                                                                                        <  4:171420/149‑1 (MQ=255)
            tggtgCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGAc                                                                                                                                        >  2:23827/1‑149 (MQ=255)
                 cgcCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTACCTACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTgg                                                                                                                                    <  1:84771/149‑1 (MQ=255)
                       gCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGGCGCat                                                                                                                               <  5:121480/147‑2 (MQ=38)
                       gCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGGCGCat                                                                                                                               <  5:222669/147‑2 (MQ=38)
                        ccGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACt                                                                                                                             >  7:143665/1‑148 (MQ=255)
                          gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCTCCGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCATGGACGCCGGGCACGTCGCCCTGGCGCACCTCaa                                                                                                                          >  2:233932/1‑149 (MQ=255)
                          gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCaa                                                                                                                          >  4:81414/1‑149 (MQ=255)
                          gaagaaATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCACAACTCaa                                                                                                                          >  5:139440/1‑149 (MQ=255)
                               aaTTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAgg                                                                                                                     <  1:319/149‑1 (MQ=255)
                                aTTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGc                                                                                                                    <  4:279567/149‑1 (MQ=255)
                                           cacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGaa                                                                                                          >  6:102450/1‑148 (MQ=255)
                                           cacCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAAc                                                                                                         <  8:143665/149‑1 (MQ=255)
                                              cGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCATTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTAc                                                                                                      <  5:7142/149‑1 (MQ=255)
                                                  cTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTc                                                                                                   <  1:23827/148‑1 (MQ=255)
                                                     aTCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGc                                                                                                <  6:143230/148‑1 (MQ=255)
                                                      tCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCaa                                                                                              <  3:81414/149‑1 (MQ=255)
                                                       ccccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGgc                                                                                                                                   <  8:286200/111‑1 (MQ=37)
                                                       ccccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCAGGTCGACGTGCCCGGCGTCCAGgc                                                                                                                                   >  7:286200/1‑111 (MQ=37)
                                                        cccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAgg                                                                                            >  8:188450/1‑149 (MQ=255)
                                                        cccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAg                                                                                             <  3:64888/148‑1 (MQ=255)
                                                         ccAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGc                                                                                           >  6:94143/1‑149 (MQ=255)
                                                                 gCCTCGGGCTTGTGCGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAAGCCTTGCCaa                                                                                   >  1:85783/1‑149 (MQ=255)
                                                                    tCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCtt                                                                                        >  4:221771/1‑141 (MQ=255)
                                                                    tCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCtt                                                                                        <  3:221771/141‑1 (MQ=255)
                                                                        gCTTGTGTGCGTTAGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCTCCGAACTACGTCAGCAAGGCCTTGCCAACCTgcgc                                                                            <  7:189516/149‑1 (MQ=255)
                                                                          ttGTGTGCGTTCGGCCAGATCCTGTCCGACGAGCGCTA‑CGACCAGCTCACCAGCCTGTCCATCCGCCTGGACGCCGGGCACTACGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACCTCAGCAAGGCCTTGCCAACCTccgcc                                                                           >  2:79556/1‑147 (MQ=255)
                                                                               tgCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTg                                                                                       >  7:79565/1‑131 (MQ=255)
                                                                               tgCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTg                                                                                       <  8:79565/131‑1 (MQ=255)
                                                                                 cGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGc                                                                   >  1:270341/1‑149 (MQ=255)
                                                                                  gTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCt                                                                  >  8:37334/1‑149 (MQ=255)
                                                                                             tCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAg                                                       <  7:37334/149‑1 (MQ=255)
                                                                                             tCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAg                                                       <  5:260114/149‑1 (MQ=255)
                                                                                              ccTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAgg                                                      >  6:109090/1‑149 (MQ=255)
                                                                                               cTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGc                                                     >  1:104974/1‑149 (MQ=255)
                                                                                                    cGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCcag                                                  >  5:153520/1‑147 (MQ=255)
                                                                                                    cGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCcag                                                  >  5:162554/1‑147 (MQ=255)
                                                                                                            gcTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCAc                                        <  2:270341/149‑1 (MQ=255)
                                                                                                              tA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTa                                      >  7:219168/1‑149 (MQ=255)
                                                                                                                    ccAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTAAGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACAcc                                  >  8:128388/1‑148 (MQ=255)
                                                                                                                      aGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTg                               <  5:94143/149‑1 (MQ=255)
                                                                                                                        cTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGga                             <  5:109090/149‑1 (MQ=255)
                                                                                                                         tCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTg                               >  5:187300/1‑146 (MQ=255)
                                                                                                                           aCCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGAt                          >  2:55412/1‑149 (MQ=255)
                                                                                                                             cAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                         <  5:253101/100‑1 (MQ=255)
                                                                                                                             cAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGaa                                                                         >  6:253101/1‑100 (MQ=255)
                                                                                                                               gCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATgcgc                      >  2:315013/1‑149 (MQ=255)
                                                                                                                                ccTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAAGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCt                     <  2:85783/149‑1 (MQ=255)
                                                                                                                                    tCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCt                 <  1:49085/149‑1 (MQ=255)
                                                                                                                                      gATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGATGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTa                                      <  1:51191/126‑1 (MQ=255)
                                                                                                                                      gATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTa                                      >  2:51191/1‑126 (MQ=255)
                                                                                                                                       aTGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTggg              <  6:162554/149‑1 (MQ=255)
                                                                                                                                       aTGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTggg              >  4:292329/1‑149 (MQ=255)
                                                                                                                                         gcgcCAGGTCGACGTGCCCGGCGTCCAGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTgggg             >  5:162412/1‑148 (MQ=255)
                                                                                                                                          cgcCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGc            <  1:233932/148‑1 (MQ=255)
                                                                                                                                           gcCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTgggg             >  2:227057/1‑146 (MQ=255)
                                                                                                                                            ccTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGa         <  2:42621/149‑1 (MQ=255)
                                                                                                                                                 aCGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCAt    <  6:153520/149‑1 (MQ=255)
                                                                                                                                                  cGCCGGGCACGTCTACCTGGCGCAACTCAACCAGGCACTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGCCCAAAAGTTGCCACTACACACTGGAGATGCGCTACCTGGGGCAGATTCAt    >  8:249860/1‑148 (MQ=255)
                                                                                                                                                   gCCGGGCACGTCGACCTGGCACAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGa  <  6:139440/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGGCATGCACGTGGTGCGCCTGGCCGAAGAAATTGGCATGCACACCGTGCTGATCCCCAAGGTGGCCTCGGGCTTGTGTGCGTTCGGCCAGATCCTGTCCGACGTGCGCTA‑CGACCAGCTCACCAGCCTGTCGATGCGCCTGGACGCCGGGCACGTCGACCTGGCGCAACTCAACCAGGCCCTCGCCGAACTACGTCAGCAAGGCCTTGCCAACCTGCGCGAAGACGGCTTTGGCGACCAGGCCAGCAGTTGCCACTACACCCTGGAGATGCGCTACCTGGGGCAGATTCATGA  >  NC_002947/3988858‑3989152

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: ‑