Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,096,289 A→C 6.5% P159P (CCT→CCG)  argT ← lysine/arginine/ornithine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,096,2890AC 6.5% 176.0 / 4.4 62P159P (CCT→CCG) argTlysine/arginine/ornithine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base C (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

GGCGTCGGTGAAGGCCGGGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGCCGTCCTTCATGACCAGGCGCGCCG  >  NC_002947/5096142‑5096434
                                                                                                                                                   |                                                                                                                                                 
ggCGTCGGTGAAGGCCGTGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCTAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTaa                                                                                                                                                   <  6:210519/148‑1 (MQ=255)
   gTCGGTGAAGGCCGGGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGca                                                                                                                                               <  5:278832/149‑1 (MQ=255)
    tCGGTGAAGGCCGGGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGcac                                                                                                                                              <  2:236016/149‑1 (MQ=255)
      ggTGAAGGCCGGGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCTTACGGcaccc                                                                                                                                            <  5:138944/149‑2 (MQ=255)
       gTGAAGGCCGGGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCTTACGGcaccca                                                                                                                                           <  1:19623/149‑3 (MQ=255)
       gTGAAGGCCGGGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGcaccac                                                                                                                                           >  4:330388/1‑149 (MQ=255)
          aaGGCCGGGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGt                                                                                                                                        <  6:372971/149‑1 (MQ=255)
               cGGGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCAc                                                                                                                                   >  8:21969/1‑149 (MQ=255)
                    cTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTg                                                                                                                              >  5:254364/1‑149 (MQ=255)
                        gAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGgcgc                                                                                                                          >  8:52320/1‑149 (MQ=255)
                        gAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGgcgc                                                                                                                          >  7:86296/1‑149 (MQ=255)
                         aaCGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGcc                                                                                                                         <  6:95863/149‑1 (MQ=255)
                          aCGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCTcc                                                                                                                         >  7:353999/1‑148 (MQ=255)
                           agctatgCCCTTGCCGGTGTCTGTCTTTAGGACGCCATCTTCCAGCAGGGTGGCTTCAGCCACTGTGCGGGCGAGGCGGCCAGCTTCAAGGTCACGTTAGATTTGTTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCa                                                                                                                       <  8:325187/143‑1 (MQ=255)
                            gcgAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAg                                                                                                                      >  5:23202/1‑149 (MQ=255)
                            gcgAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCAACTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTACAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAg                                                                                                                      >  2:2605/1‑149 (MQ=255)
                              gAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCa                                                                                                                    >  5:91996/1‑149 (MQ=255)
                                       gCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTg                                                                                                            <  7:52320/148‑1 (MQ=255)
                                           gCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCg                                                                                                         >  1:25142/1‑147 (MQ=255)
                                               cAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCgg                                                                                                   >  1:52437/1‑149 (MQ=255)
                                                aGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCAGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCgg                                                                                                   <  2:127560/148‑1 (MQ=255)
                                                         gAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAg                                                                                                      >  5:61948/1‑136 (MQ=255)
                                                          aaGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGa                                                                                        <  7:21969/149‑1 (MQ=255)
                                                              cATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGaa                                                                                            <  7:324093/141‑1 (MQ=255)
                                                              cATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGaa                                                                                            >  8:324093/1‑141 (MQ=255)
                                                              cATCTTCCAGCAGGCTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTGCAGGTTGATTTCCTTCTCGGTGCCGTAAGGCCCCACGGTGGCACCCTTGGCGCCGAGCACTTCCTTGGCGAAGCGGTCGTGATAGGAGCcg                                                                                    >  8:153974/1‑149 (MQ=255)
                                                               aTCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTgg                                                                                                           >  5:296153/1‑125 (MQ=255)
                                                               aTCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTgg                                                                                                           <  6:296153/125‑1 (MQ=255)
                                                                 cTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCgcgc                                                                                 >  1:293261/1‑149 (MQ=255)
                                                                 cTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCGGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTg                                                                                                            <  7:120805/122‑1 (MQ=255)
                                                                 cTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCGGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTg                                                                                                            >  8:120805/1‑122 (MQ=255)
                                                                    ccAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCt                                                                                <  2:52437/147‑1 (MQ=255)
                                                                    ccAGCAGGGTGGCGTCACCAACAGTGCCGGCGCGGCGGCCAGCTGCACCGTCCAGGTATATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGc                                                                              <  6:19032/149‑1 (MQ=255)
                                                                       gcagGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTgg                                                                                                           >  8:93498/1‑117 (MQ=255)
                                                                            gTGGCGTCAGCCAAAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCTTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGAt                                                                      <  7:150966/149‑1 (MQ=255)
                                                                               gCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGAtctt                                                                   <  1:36351/149‑1 (MQ=255)
                                                                                cggcAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGAtcttc                                                                  <  7:68597/146‑1 (MQ=255)
                                                                                   cAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATcttctt                                                                <  6:91996/148‑1 (MQ=255)
                                                                                       cacaGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGc                                                                          <  1:277527/134‑1 (MQ=255)
                                                                                       cacaGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGc                                                                          >  2:277527/1‑134 (MQ=255)
                                                                                        acaGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCttt                                                          <  7:76096/149‑1 (MQ=255)
                                                                                          aGTGCCGTGGAGGCGTCTCCTTGCAACCGCCAGGTAGATTACGTTGTGGGTGCCGTAAGGCACCTCGGTGGCACGCTTGGCGGCCAGCTCTTCCTTGGCGTAGCGGTCGGGAATGGAGCCGCGGTGCACGCCGATCTTCTTGCCttt                                                          <  7:153974/147‑1 (MQ=255)
                                                                                              ccGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTc                                                    <  1:2605/149‑1 (MQ=255)
                                                                                                    aGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAg                                              <  1:281125/149‑1 (MQ=255)
                                                                                                     ggcggcCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAgg                                             <  3:330388/149‑1 (MQ=255)
                                                                                                      gcggcCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAgg                                             <  8:86296/148‑1 (MQ=255)
                                                                                                      gcggcCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGc                                            >  4:245361/1‑149 (MQ=255)
                                                                                                         gcCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCCCTGCACGCCGATCTTATTGCCTTTCAGTTCGTCCAGGCAGt                                         >  2:160124/1‑149 (MQ=255)
                                                                                                            aGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGc                                      >  6:239209/1‑149 (MQ=255)
                                                                                                                   aCGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATcttctt                                                                >  6:36704/1‑116 (MQ=255)
                                                                                                                   aCGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATcttctt                                                                <  5:36704/116‑1 (MQ=255)
                                                                                                                        cAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGgtggtg                           <  6:23202/148‑1 (MQ=255)
                                                                                                                         aGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGcc                         >  5:213846/1‑149 (MQ=255)
                                                                                                                           gTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGCCGt                       >  3:43417/1‑149 (MQ=255)
                                                                                                                           gTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGCCGt                       >  7:134790/1‑149 (MQ=255)
                                                                                                                            tAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGCCGTc                      <  2:139612/149‑1 (MQ=255)
                                                                                                                                  tCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCTCGCTGCACGCCTATCTTCTTGCCTTTCAGTTCTTCCAGGCTGTCGCTGACGGTGGTGCCGTCCTTCa                 <  7:184534/148‑1 (MQ=255)
                                                                                                                                     ttCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGCCGTCCTTCATGAc             >  5:72059/1‑149 (MQ=255)
                                                                                                                                     accgtggTGCCTTACGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGCCGTCCTTCATGAc             >  5:146723/6‑149 (MQ=255)
                                                                                                                                       cgtggTGCCTTAAGGCACCACGGTGGCACCCATGGCGCCCAGCACTTACTTGGCGAACCGGGCGTGAATGGAGCCGCGCTGCCCCCCCATCGTCTTGCCGTTGAGTTCGTCCAGGCTGGCGCTGAGGGTGGTGCCCGACTTCATGACCa           >  4:166182/4‑149 (MQ=255)
                                                                                                                                         tggTGCCTTACGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGCCGTCCTTCATGACCAgg         >  5:82787/2‑149 (MQ=255)
                                                                                                                                               cGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGCCTGTCGCTGACGGTGGTGCCGTCCTTCATGACCAGGCGCGcc   >  5:167057/1‑149 (MQ=255)
                                                                                                                                                gTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGCCGTCCTTCATGACCAGGCGCGCCg  <  6:72059/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCGTCGGTGAAGGCCGGGCCTACGAACGCGAAGCCCTTGCCGGCGTCAGTCTTGAGGAAGCCATCTTCCAGCAGGGTGGCGTCAGCCACAGTGCCGTCGAGGCGGCCAGCTGCAACGTCCAGGTAGATTTCGTTCTGGGTGCCGTAAGGCACCACGGTGGCACCCTTGGCGCCCAGCACTTCCTTGGCGAAGCGGTCGTGAATGGAGCCGCGCTGCACGCCGATCTTCTTGCCTTTCAGTTCGTCCAGGCTGTCGCTGACGGTGGTGCCGTCCTTCATGACCAGGCGCGCCG  >  NC_002947/5096142‑5096434

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