Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,868,876 G→C 8.5% intergenic (+415/‑516) PP_2520 → / → PP_2522 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,868,8760GC 8.5% 176.6 / 6.9 59intergenic (+415/‑516)PP_2520/PP_2522hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base C (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

AAAGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTCGCGCGCG  >  NC_002947/2868728‑2869021
                                                                                                                                                    |                                                                                                                                                  
aaaGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCg                                                                                                                                                    >  2:164449/1‑149 (MQ=255)
 aaGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCgg                                                                                                                                                   <  3:186124/149‑1 (MQ=255)
  aGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCggg                                                                                                                                                  <  2:30694/149‑1 (MQ=255)
       ttCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCAtt                                                                                                                                             >  2:192082/1‑149 (MQ=255)
         cAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGATCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCCCGGCATTga                                                                                                                                           <  6:132038/149‑2 (MQ=255)
           gCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGt                                                                                                                                          >  4:3810/1‑148 (MQ=255)
            ccTAGTTCAGGCCATTGGCCATTCGGGTCAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGa                                                                                                                                        >  1:28930/1‑149 (MQ=255)
                gTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTg                                                                                                                                    <  3:17334/149‑1 (MQ=255)
                 ttCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGc                                                                                                                                   >  3:403440/1‑149 (MQ=255)
                   cAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCg                                                                                                                                  <  8:223386/148‑1 (MQ=255)
                    ctgccatttggcaatGGGGTGTAGCCATTTTGGAGCGGCGGGGGCATGTCGCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGa                                                                                                                                 <  8:281163/146‑1 (MQ=255)
                    aGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGTGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGa                                                                                                                                 <  1:225859/148‑1 (MQ=255)
                     ggccattggccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGata                                                                                                                               >  1:118751/1‑149 (MQ=255)
                        cattggccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATatc                                                                                                                             >  1:196439/1‑148 (MQ=255)
                          ttggccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATatcaat                                                                                                                          <  7:62827/149‑1 (MQ=255)
                          ttggccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATatcaat                                                                                                                          <  3:3810/149‑1 (MQ=255)
                            ggccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATatcaatc                                                                                                                         <  1:164449/148‑1 (MQ=255)
                                attCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTcg                                                                                                                    >  1:96893/1‑149 (MQ=255)
                                         aaGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGa                                                                                                                                 <  1:153202/127‑1 (MQ=255)
                                         aaGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGa                                                                                                                                 >  2:153202/1‑127 (MQ=255)
                                            ccATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGcc                                                                                                        >  2:189531/1‑149 (MQ=255)
                                              aTTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGACATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGcccc                                                                                                      >  4:198430/1‑149 (MQ=255)
                                                  gggAGCGGCTGGGGCATGAGCCGTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACg                                                                                                  >  7:256412/1‑149 (MQ=255)
                                                   ggAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGt                                                                                                 >  2:94376/1‑149 (MQ=255)
                                                       gggCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCATTCAGTCGCGCTTCGCGGCCCCCACGTTGt                                                                                              <  8:121073/147‑1 (MQ=255)
                                                       cGGCTGGGTCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGGTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGt                                                                                              <  7:147628/148‑1 (MQ=255)
                                                         gCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCAcc                                                                                                  <  6:407005/142‑2 (MQ=255)
                                                         gCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCAcc                                                                                                  >  5:407005/1‑141 (MQ=255)
                                                         gCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCACCAATGCCCCGGCATtgatgaacg                                                                                                                                     <  6:46672/107‑4 (MQ=39)
                                                         gCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCACCAATGCCCCGGCATtgatgaacg                                                                                                                                     >  5:46672/1‑104 (MQ=39)
                                                              ggCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATgcgc                                                                                      <  2:196439/149‑1 (MQ=255)
                                                               gCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATgcgc                                                                                      >  8:118391/1‑148 (MQ=255)
                                                                       cTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTc                                                                            <  2:96893/149‑1 (MQ=255)
                                                                                aGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCtt                                                                   >  6:45573/1‑149 (MQ=255)
                                                                                   cGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCa                                                                <  4:403440/149‑1 (MQ=255)
                                                                                     ttCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAgg                                                              >  4:56829/1‑149 (MQ=255)
                                                                                     ttCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAgg                                                              >  6:100928/1‑149 (MQ=255)
                                                                                       cAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGc                                                             >  8:26872/1‑148 (MQ=255)
                                                                                        aaCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCaa                                                           >  2:50537/1‑149 (MQ=255)
                                                                                         aCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGTGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAAc                                                          <  5:100928/149‑1 (MQ=255)
                                                                                         aCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCTCGATTTTGTCCGACGGCTTTCAGGCAAc                                                          <  7:26872/149‑1 (MQ=255)
                                                                                         aCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAAc                                                          <  1:50537/149‑1 (MQ=255)
                                                                                         aCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGAATTTGTCCGACGGCTTTCAGGCAAc                                                          >  6:417072/1‑149 (MQ=255)
                                                                                                ggTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTcac                                                    <  5:45573/148‑1 (MQ=255)
                                                                                                   gCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACAtct                                                >  6:50648/1‑149 (MQ=255)
                                                                                                    cAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACAtct                                                <  8:256412/148‑1 (MQ=255)
                                                                                                    cAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGACATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACAtct                                                <  3:198430/148‑1 (MQ=255)
                                                                                                      gCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGa                                             >  8:163343/1‑149 (MQ=255)
                                                                                                       cTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGAt                                            >  7:247405/1‑149 (MQ=255)
                                                                                                                 aGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCgg                                  >  5:168277/1‑149 (MQ=255)
                                                                                                                 aGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCgg                                  <  7:163343/149‑1 (MQ=255)
                                                                                                                  gcgcGGGCGGCCGCGTAACCCGTTCATCAATGCCGTGGCATTGGTGATCTGCGATATCAATCAGTCGCTCTTCGCGGCCCACACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATc                                           <  3:32595/139‑1 (MQ=255)
                                                                                                                  gcgcGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCCATTTTGTCCGACGGCTTTCAGGCAACCATCACAtctctaaa                                           >  4:32595/1‑135 (MQ=255)
                                                                                                                      gggcggCCGCGTAACCCGTTCATCCATGCCGGGGCATTGGTGATCTGCGATATCCATCAGTCGCGCGTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGc                             <  6:378219/149‑1 (MQ=255)
                                                                                                                        gcggcCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCt                            <  1:94376/149‑1 (MQ=255)
                                                                                                                            ccGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTc                       >  7:51040/1‑149 (MQ=255)
                                                                                                                            ccGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTc                       <  3:56829/149‑1 (MQ=255)
                                                                                                                            ccGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTc                       <  2:118751/149‑1 (MQ=255)
                                                                                                                                     ccGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAatat                >  6:321240/1‑147 (MQ=255)
                                                                                                                                            ccaATGCCCCGGCATTGATGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTc         <  2:24689/146‑1 (MQ=255)
                                                                                                                                            ccaATGCCCCGGCATTGATGATCTGCGAAATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTc         >  1:24689/2‑147 (MQ=255)
                                                                                                                                                  ccGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTcgcgcgcg  >  2:226495/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAAGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTCGCGCGCG  >  NC_002947/2868728‑2869021

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: