Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,868,877 G→C 8.5% intergenic (+416/‑515) 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,8770GC 8.5% 168.6 / 6.8 59intergenic (+416/‑515)PP_2520/PP_2522hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base C (2/3);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

AAGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTCGCGCGCG  >  NC_002947/2868729‑2869021
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                                                 >  2:153202/1‑127 (MQ=255)
                                        aaGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGa                                                                                                                                 <  1:153202/127‑1 (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)
                                                                      cTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑Tc                                                                            <  2:96893/149‑1 (MQ=255)
                                                                               aGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCtt                                                                   >  6:45573/1‑149 (MQ=255)
                                                                                  cGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCa                                                                <  4:403440/149‑1 (MQ=255)
                                                                                    ttCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAgg                                                              >  4:56829/1‑149 (MQ=255)
                                                                                    ttCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAgg                                                              >  6:100928/1‑149 (MQ=255)
                                                                                      cAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGc                                                             >  8:26872/1‑148 (MQ=255)
                                                                                       aaCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCaa                                                           >  2:50537/1‑149 (MQ=255)
                                                                                        aCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGTGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAAc                                                          <  5:100928/149‑1 (MQ=255)
                                                                                        aCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCTCGATTTTG‑TCCGACGGCTTTCAGGCAAc                                                          <  7:26872/149‑1 (MQ=255)
                                                                                        aCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAAc                                                          <  1:50537/149‑1 (MQ=255)
                                                                                        aCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGAATTTG‑TCCGACGGCTTTCAGGCAAc                                                          >  6:417072/1‑149 (MQ=255)
                                                                                               ggTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTcac                                                    <  5:45573/148‑1 (MQ=255)
                                                                                                  gCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACAtct                                                >  6:50648/1‑149 (MQ=255)
                                                                                                   cAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACAtct                                                <  8:256412/148‑1 (MQ=255)
                                                                                                   cAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGACATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACAtct                                                <  3:198430/148‑1 (MQ=255)
                                                                                                     gCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGa                                             >  8:163343/1‑149 (MQ=255)
                                                                                                      cTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGAt                                            >  7:247405/1‑149 (MQ=255)
                                                                                                                aGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCgg                                  <  7:163343/149‑1 (MQ=255)
                                                                                                                aGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCgg                                  >  5:168277/1‑149 (MQ=255)
                                                                                                                 gcgcGGGCGGCCGCGTAACCCGTTCATCAATGCCGTGGCATTGGTGATCTGCGATATCAATCAGTCGCTCTTCGCGGCCCACACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATc                                           <  3:32595/139‑1 (MQ=255)
                                                                                                                 gcgcGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCCATTTTG‑TCCGACGGCTTTCAGGCAACCATCACAtctctaaa                                           >  4:32595/1‑135 (MQ=255)
                                                                                                                     gggcggCCGCGTAACCCGTTCATCCATGCCGGGGCATTGGTGATCTGCGATATCCATCAGTCGCGCGTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGc                             <  6:378219/149‑1 (MQ=255)
                                                                                                                       gcggcCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCt                            <  1:94376/149‑1 (MQ=255)
                                                                                                                           ccGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTc                       <  2:118751/149‑1 (MQ=255)
                                                                                                                           ccGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTc                       >  7:51040/1‑149 (MQ=255)
                                                                                                                           ccGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTc                       <  3:56829/149‑1 (MQ=255)
                                                                                                                                    ccGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAatat                >  6:321240/1‑147 (MQ=255)
                                                                                                                                           ccaATGCCCCGGCATTGATGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTc         <  2:24689/146‑1 (MQ=255)
                                                                                                                                           ccaATGCCCCGGCATTGATGATCTGCGAAATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTc         >  1:24689/2‑147 (MQ=255)
                                                                                                                                                 ccGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTcgcgcgcg  >  2:226495/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AAGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTG‑TCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTCGCGCGCG  >  NC_002947/2868729‑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: ‑