Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,227,419 C→T 7.1% G547D (GGC→GAC)  PP_4606 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,227,4190CT 7.1% 167.7 / 2.4 58G547D (GGC→GAC) PP_4606ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (22/31);  new base T (2/2);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTGCGCAGCTTGTCGGT  >  NC_002947/5227271‑5227567
                                                                                                                                                    |                                                                                                                                                    
ctcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGc                                                                                                                                                      <  4:244455/149‑1 (MQ=255)
ctcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGc                                                                                                                                                      <  2:139164/149‑1 (MQ=255)
 tcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGc                                                                                                                                                      >  6:268133/1‑148 (MQ=255)
   aGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGTCAc                                                                                                                                                   >  5:71892/1‑149 (MQ=255)
   aGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCAc                                                                                                                                                   <  3:235498/149‑1 (MQ=255)
   aGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTAGACTTTGTAAACGATACCGGCGTTAGGTAACCAGGCCTGGTCAc                                                                                                                                                   >  5:66554/1‑149 (MQ=255)
     gTGCCACCGAGGAGTCGAGCACGACATTGCCGTTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCAc                                                                                                                                                   <  3:76298/147‑1 (MQ=255)
       gCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGAt                                                                                                                                               >  1:256640/1‑149 (MQ=255)
        ccACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATg                                                                                                                                              >  6:283524/1‑149 (MQ=255)
          aCCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTc                                                                                                                                            >  8:172815/1‑149 (MQ=255)
           ccGAGGAGTCGAGCACTACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCg                                                                                                                                           <  6:240923/149‑1 (MQ=255)
                  gTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCTAGTTGGGCTTGACCGACTCGCTGTAGCTGCCATAGAACGAGTGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTgg                                                                                                                                    <  5:343596/149‑1 (MQ=255)
                       gCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTg                                                                                                                               <  8:8998/149‑1 (MQ=255)
                           gACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGCAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACg                                                                                                                           <  4:178560/149‑1 (MQ=255)
                             cATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACggg                                                                                                                         <  6:288929/149‑1 (MQ=255)
                                   cGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACcg                                                                                                                   >  8:137473/1‑149 (MQ=255)
                                        aGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGcc                                                                                                              <  2:256640/149‑1 (MQ=255)
                                         gcggcgCGATGGTCGAGTTGGGCTTGAATGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCa                                                                                                             <  6:149474/149‑1 (MQ=255)
                                              gagATGGTAGAGTTGTGCTTGAACGAATCGCTGTTGCTGCCATATAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTa                                                                                                        <  8:241185/147‑1 (MQ=255)
                                                     tCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCg                                                                                                 <  7:172815/149‑1 (MQ=255)
                                                        aGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCa                                                                                                             <  5:78559/134‑1 (MQ=255)
                                                        aGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCa                                                                                                             >  6:78559/1‑134 (MQ=255)
                                                              gCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTgg                                                                                        <  5:283524/149‑1 (MQ=255)
                                                                   atggACTCGCTGTAGCTGCCATAGAACGATAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAAc                                                                                    <  4:80239/145‑1 (MQ=255)
                                                                   aaCGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAAcg                                                                                   >  8:198003/1‑149 (MQ=255)
                                                                         tCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCAcc                                                                             >  4:48541/1‑149 (MQ=255)
                                                                            cTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGc                                                                          <  2:71852/149‑1 (MQ=255)
                                                                                 gCTGCCCTAGAACCACAGCTGGTCGTCGACTTTGTAAACGATCCCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACaa                                                                     <  8:233128/149‑1 (MQ=255)
                                                                                 gCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACaa                                                                     <  6:235428/149‑1 (MQ=255)
                                                                                     ccATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGAt                                                                 >  2:176330/1‑149 (MQ=255)
                                                                                      cATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTCCCCAGGCCTGGACACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATc                                                                <  7:198003/149‑1 (MQ=255)
                                                                                       aTAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATc                                                                >  3:97245/1‑148 (MQ=255)
                                                                                       aTAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATc                                                                >  7:1827/1‑148 (MQ=255)
                                                                                           aaCGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTg                                                           >  3:78008/1‑149 (MQ=255)
                                                                                                  gCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCg                                                    <  7:137473/149‑1 (MQ=255)
                                                                                                      gtcgtcGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAgg                                                 <  5:303637/148‑1 (MQ=255)
                                                                                                      gtcgtcGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAgg                                                 <  1:324446/148‑1 (MQ=255)
                                                                                                       tcgtcgACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGTCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTg                                               <  6:66554/149‑1 (MQ=255)
                                                                                                       tcgtcgACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGTCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTg                                               <  6:71892/149‑1 (MQ=255)
                                                                                                           cgACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACAGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAgc                                           <  7:295578/149‑1 (MQ=255)
                                                                                                            gACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCTCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAgcg                                          <  5:268133/149‑1 (MQ=255)
                                                                                                             aCTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAgcgc                                         >  4:86287/1‑149 (MQ=255)
                                                                                                                  gTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCt                                    >  1:188201/1‑149 (MQ=255)
                                                                                                                   tAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTc                                      <  2:228726/146‑1 (MQ=255)
                                                                                                                   tAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTc                                      >  1:228726/1‑146 (MQ=255)
                                                                                                                    aaaCGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGGACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTg                                   >  6:9100/1‑148 (MQ=255)
                                                                                                                          tACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGaa                             >  2:360416/1‑148 (MQ=255)
                                                                                                                            ccGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACa                           <  3:48541/148‑1 (MQ=255)
                                                                                                                                   gAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAgcgaccag                   >  4:314384/1‑144 (MQ=255)
                                                                                                                                    aGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCgg                  >  5:191151/1‑149 (MQ=255)
                                                                                                                                     ggTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGt                 <  8:1827/149‑1 (MQ=255)
                                                                                                                                     ggTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGt                 >  1:216529/1‑149 (MQ=255)
                                                                                                                                       tACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTgc               <  6:332484/149‑1 (MQ=255)
                                                                                                                                          ccAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGGACGCGCTGCAcctcgat                          <  2:182128/135‑7 (MQ=255)
                                                                                                                                          ccAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGGACGCGCTGCAcctcgat                          >  1:182128/1‑129 (MQ=255)
                                                                                                                                           cAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTGCGCAg           <  1:176330/149‑1 (MQ=255)
                                                                                                                                            aGGCCTGGCCACTGATGTCGGTATTGGCCTTGGACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTGCGCAGc          <  5:9100/149‑1 (MQ=255)
                                                                                                                                            aGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTGCGCAGc          <  2:364298/149‑1 (MQ=255)
                                                                                                                                              gCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTGCGCAGCtt        <  1:360416/149‑1 (MQ=255)
                                                                                                                                                cTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTGCGCAGCTTGt      >  7:68296/1‑149 (MQ=255)
                                                                                                                                                cTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTGCGCAGCTTGt      >  2:18953/1‑149 (MQ=255)
                                                                                                                                                    ccACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTGCGCAGCTTGTCGGt  <  1:148892/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTCGGTATTGGCCTTGAACGGGCGACCGCGGCCAGCGTACTGGTCGTACTGCTGGAAACGCGCACCTGCTACAAGGATCCAGTGCTCATCGAGGTGCAGCGCGTCCTGCACGAACAGCGCATCGGTGCGCAGCTTGTCGGT  >  NC_002947/5227271‑5227567

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