Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 548,804 G→T 6.4% V552F (GTC→TTC)  fusA → elongation factor G 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947548,8040GT 6.4% 198.6 / 3.8 63V552F (GTC→TTC) fusAelongation factor G 1
Reads supporting (aligned to +/- strand):  ref base G (34/25);  new base T (2/2);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

ACCAAATCCAACGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGA  >  NC_002947/548660‑548946
                                                                                                                                                |                                                                                                                                               
aCCAAATCCAACGTCGAGATCGAAGGCAAGTTAGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCAATGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGt                                                                                                                                             >  7:95097/1‑149 (MQ=255)
 ccAAATCCAACGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGtt                                                                                                                                            >  4:336486/1‑149 (MQ=255)
 ccAAATCCAACGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGtt                                                                                                                                            >  4:103692/1‑149 (MQ=255)
  cAAATCCAACGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTtg                                                                                                                                           >  3:152627/1‑149 (MQ=255)
    aaTCCAACGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGACTTCGTtgga                                                                                                                                         <  4:121211/149‑2 (MQ=255)
     aTCCAACGTCGAGATCGAAGGCACGTGCGTTCGGCAGTCGGGTGGTCGTCGTCAGTTCGGTCACTGCTGGGATCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTtggtg                                                                                                                                        <  2:367338/149‑1 (MQ=255)
        cAACGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGTCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGgtggtg                                                                                                                                     >  2:343364/1‑149 (MQ=255)
        cAACGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGgtggtg                                                                                                                                     >  2:16009/1‑149 (MQ=255)
        cAACGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGACGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGgtggtg                                                                                                                                     >  3:40906/1‑149 (MQ=255)
          aCGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCACCAACGACTTCGTtggagaacac                                                                                                                                   <  2:301107/149‑8 (MQ=255)
           cGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGa                                                                                                                                  <  1:16009/149‑1 (MQ=255)
               gagaTCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGt                                                                                                                                             >  1:70777/1‑134 (MQ=255)
               gagaTCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGt                                                                                                                                             <  2:70777/134‑1 (MQ=255)
                      aaGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGaa                                                                                                                       >  1:273534/1‑149 (MQ=255)
                           aaGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACAt                                                                                                                  <  4:152627/149‑1 (MQ=255)
                                cgttcgtCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGAt                                                                                                                                 >  5:95790/1‑129 (MQ=255)
                                cgttcgtCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGAt                                                                                                                                 <  6:95790/129‑1 (MQ=255)
                                            gggTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAgg                                                                                                 <  6:139080/149‑1 (MQ=255)
                                            gggTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAgg                                                                                                 <  3:103692/149‑1 (MQ=255)
                                                    gtggtcAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGTCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGaa                                                                                         <  1:343364/149‑1 (MQ=255)
                                                    gtggtcAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCg                                                                                                                                              <  4:128944/96‑1 (MQ=255)
                                                    gtggtcAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCg                                                                                                                                              >  3:128944/1‑96 (MQ=255)
                                                          aGTGCGTTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCa                                                                                    <  8:95097/148‑1 (MQ=255)
                                                            ttCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATGCCTAAGGAATACATCCCTGCGATCCAGAAGGGAATCGAAGAGCAGAt                                                                                 >  2:67535/1‑149 (MQ=255)
                                                                     tgctgGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACgg                                                                        >  1:24534/1‑149 (MQ=255)
                                                                      gctgGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGc                                                                       >  5:112670/1‑149 (MQ=255)
                                                                       ctgGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGaa                                                                                                   <  2:305555/120‑1 (MQ=255)
                                                                       ctgGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGaa                                                                                                   >  1:305555/1‑120 (MQ=255)
                                                                                cTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCg                                                            >  7:317018/1‑149 (MQ=255)
                                                                                      cGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATc                                                      <  3:336486/149‑1 (MQ=255)
                                                                                           cGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCg                                                            <  7:355172/138‑1 (MQ=255)
                                                                                           cGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCg                                                            >  8:355172/1‑138 (MQ=255)
                                                                                               cGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCg                                             >  4:125438/1‑149 (MQ=255)
                                                                                                   gacgacAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGcc                                          <  2:24534/148‑1 (MQ=255)
                                                                                                   gacgacAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCt                                         >  8:130646/1‑149 (MQ=255)
                                                                                                   gacgacAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCt                                         >  4:276736/1‑149 (MQ=255)
                                                                                                        cAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAgg                                    >  3:138759/1‑149 (MQ=255)
                                                                                                           gggCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGtt                                                         <  8:120672/125‑1 (MQ=255)
                                                                                                           gggCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGtt                                                         >  7:120672/1‑125 (MQ=255)
                                                                                                           gggCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATc                                                      >  1:346761/1‑128 (MQ=255)
                                                                                                           gggCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATc                                                      <  2:346761/128‑1 (MQ=255)
                                                                                                           gggCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCcg                                 >  5:289729/1‑149 (MQ=255)
                                                                                                              cAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGAATCGAAGAGCAGATGAAAAACGGCGTCTGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCcgcg                               <  1:67535/149‑1 (MQ=255)
                                                                                                                 cATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTAt                           >  3:243458/1‑149 (MQ=255)
                                                                                                                     acaGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGt                                                                                               >  4:170380/1‑78 (MQ=255)
                                                                                                                     acaGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGt                                                                                               <  3:170380/78‑1 (MQ=255)
                                                                                                                     acaGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGAAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGa                       >  8:341884/1‑149 (MQ=255)
                                                                                                                        gAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCATAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAATAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACgg                    <  4:40906/149‑1 (MQ=255)
                                                                                                                            gTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCctc                                                   <  4:118179/114‑1 (MQ=255)
                                                                                                                            gTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCctc                                                   >  3:118179/1‑114 (MQ=255)
                                                                                                                                   gtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTc                         <  1:294972/133‑1 (MQ=255)
                                                                                                                                   gtTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTc                         >  2:294972/1‑133 (MQ=255)
                                                                                                                                   caccaccAACGACTTCGTTGGAGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGACGCGGTATTCCACGGTTCGTACCac          >  7:252311/6‑148 (MQ=255)
                                                                                                                                      ctcCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCacgac       <  2:273534/148‑1 (MQ=255)
                                                                                                                                      caccAACGACTTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCacgacg      >  1:181407/3‑149 (MQ=255)
                                                                                                                                        ccAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTc    >  2:358341/1‑149 (MQ=255)
                                                                                                                                         cAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGt     <  3:125438/147‑1 (MQ=255)
                                                                                                                                         cAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTc    >  6:58677/1‑148 (MQ=255)
                                                                                                                                         cAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTc    <  6:289729/148‑1 (MQ=255)
                                                                                                                                         cAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTc    >  3:286885/1‑148 (MQ=255)
                                                                                                                                          aaCGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGa  >  8:170433/1‑149 (MQ=255)
                                                                                                                                            cGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCg                >  2:295356/1‑133 (MQ=255)
                                                                                                                                            cGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACCTCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCg                <  1:295356/133‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ACCAAATCCAACGTCGAGATCGAAGGCAAGTTCGTTCGTCAGTCGGGTGGTCGTGGTCAGTTCGGTCACTGCTGGGTTCGCTTCTCCGAGCCGGACGTGGACGACAAGGGCAACATCACAGAAGGTCTGGTGTTCTCCAACGAAGTCGTTGGTGGTGTGATTCCTAAGGAATACATCCCTGCGATCCAGAAGGGTATCGAAGAGCAGATGAAAAACGGCGT‑TGTTGCCGGTTATCCTCTGATCGGCCTGAAGGCCGCGGTATTCGACGGTTCGTACCACGACGTCGA  >  NC_002947/548660‑548946

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