Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,238,144 A→T 8.6% M147K (ATG→AAG) ‡ PP_1975 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,238,1440AT 8.6% 163.2 / 5.5 58M147K (ATG→AAG) ‡PP_1975EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base T (3/2);  total (29/29)
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.22e-01

CTTCTCAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGCACAGCAGCG  >  NC_002947/2237996‑2238291
                                                                                                                                                    |                                                                                                                                                     
cttctCAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACa                                                                                                                                                       >  8:78325/1‑149 (MQ=255)
   ctcAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTTTTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCg                                                                                                                                                    >  6:211356/1‑149 (MQ=255)
      aaGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCa                                                                                                                                                 >  6:171751/1‑149 (MQ=255)
      aaGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCa                                                                                                                                                 <  1:52964/149‑1 (MQ=255)
            gTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCt‑‑t                                                                                                                                         >  2:163934/1‑149 (MQ=255)
                aaCAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGgc                                                                                                                                     >  4:347834/1‑149 (MQ=255)
                   aaGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGcgcg                                                                                                                                  >  2:294069/1‑149 (MQ=255)
                   aaGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGACGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGagcg                                                                                                                                  >  4:302423/1‑149 (MQ=255)
                     gCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGTCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCggt                                                                                                                                <  2:57558/149‑1 (MQ=255)
                      ccGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCggtg                                                                                                                               <  7:78325/149‑1 (MQ=255)
                            aTCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGGGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTc                                                                                                                         <  2:238540/149‑1 (MQ=255)
                             tCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTTTTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCg                                                                                                                        <  5:211356/149‑1 (MQ=255)
                               aTTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGAGCGGTGGTGCTCggg                                                                                                                      >  7:324369/1‑149 (MQ=255)
                                ttAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCggg                                                                                                                      >  1:250187/1‑148 (MQ=255)
                                   gCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGc                                                                                                                  <  6:177402/149‑1 (MQ=255)
                                   gCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGc                                                                                                                  >  4:272215/1‑149 (MQ=255)
                                    cAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGATGTTCGCCATgacc                                                                                                                                            >  7:16208/1‑121 (MQ=255)
                                    cAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGATGTTCGCCATgacc                                                                                                                                            <  8:16208/124‑4 (MQ=255)
                                            gTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGc                                                                                                         <  1:163934/149‑1 (MQ=255)
                                                   atTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAggg                                                                                                  >  7:61142/1‑149 (MQ=255)
                                                      gATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTg                                                                                                >  6:312061/1‑148 (MQ=255)
                                                         gTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGGGTAGGTTGGGCGGGGGTGGTCGGGCCGCTTGTTCAGCGGCTGGGGGat                                                                                              >  1:124212/1‑148 (MQ=255)
                                                             aagaTGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAg                                                                                        >  6:287582/1‑149 (MQ=255)
                                                               gaTGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCg                                                                                      >  5:153849/1‑149 (MQ=255)
                                                                aTGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGACGCTGGAGCCGCGGCGGGGGGAAGCCAGCcgg                                                                                    >  3:233074/1‑149 (MQ=255)
                                                                        ccAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGTATCAgcg                                                                             >  7:35110/1‑149 (MQ=255)
                                                                        ccAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAgcg                                                                             >  3:358905/1‑149 (MQ=255)
                                                                         cAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAgcgc                                                                            <  5:287582/149‑1 (MQ=255)
                                                                         cAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAgcgc                                                                            <  3:272215/149‑1 (MQ=255)
                                                                          aGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAgcgcg                                                                           >  7:29136/1‑149 (MQ=255)
                                                                             cccAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGc                                                                                                         <  7:262827/116‑1 (MQ=255)
                                                                             cccAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGc                                                                                                         >  8:262827/1‑116 (MQ=255)
                                                                                       ggTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGaaa                                                               <  3:347834/148‑1 (MQ=255)
                                                                                         tCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAAccc                                                            >  1:141467/1‑149 (MQ=255)
                                                                                           agccagccGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTg                                                          <  2:250187/149‑1 (MQ=255)
                                                                                             ccagccGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGc                                                         >  8:194816/1‑148 (MQ=255)
                                                                                                gccgccCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAgg                                                     >  6:283742/1‑149 (MQ=255)
                                                                                                    cccAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACg                                                 <  5:312061/149‑1 (MQ=255)
                                                                                                         ggtgggtCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACa                                            <  8:29136/149‑1 (MQ=255)
                                                                                                           tgggtCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAgg                                          <  1:294069/149‑1 (MQ=255)
                                                                                                             ggtCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCAc                                                                       <  2:337089/118‑1 (MQ=255)
                                                                                                             ggtCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCAc                                                                       >  1:337089/1‑118 (MQ=255)
                                                                                                              gtCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGaa                                                                <  2:336061/124‑1 (MQ=255)
                                                                                                              gtCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGaa                                                                >  1:336061/1‑124 (MQ=255)
                                                                                                              gtCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATc                                       <  6:153849/149‑1 (MQ=255)
                                                                                                                 cGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCAt                                     <  7:194816/148‑1 (MQ=255)
                                                                                                                    tGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGTATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGctc                                 <  8:35110/149‑1 (MQ=255)
                                                                                                                           ggCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAg                                           <  3:158210/132‑1 (MQ=255)
                                                                                                                           ggCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGt                          <  7:267003/149‑1 (MQ=255)
                                                                                                                               aaaGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCa                      <  2:124212/149‑1 (MQ=255)
                                                                                                                               aaaGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCa                      <  5:68482/149‑1 (MQ=255)
                                                                                                                                aaGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAg                     <  2:141467/149‑1 (MQ=255)
                                                                                                                                aaGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAg                     <  5:171751/149‑1 (MQ=255)
                                                                                                                                aaGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAg                     <  8:61142/149‑1 (MQ=255)
                                                                                                                                      aaggcATGGCGATGTTCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATg                                    <  6:61616/124‑1 (MQ=255)
                                                                                                                                      aaggcATGGCGATGTTCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATg                                    >  5:61616/5‑128 (MQ=255)
                                                                                                                                      aaggcATGGCGATGTTCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATc               >  3:351772/5‑149 (MQ=255)
                                                                                                                                                gAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGCAcag       >  7:354586/1‑147 (MQ=255)
                                                                                                                                                   cATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCCCCGGTGAAACCCTGCAAGGCTCGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGCACAGCAGcg  >  5:356906/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTCTCAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCT‑‑TTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGCACAGCAGCG  >  NC_002947/2237996‑2238291

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: