Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,550,711 A→G 8.0% D302G (GAC→GGC)  PP_4881 → iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,550,7110AG 8.0% 191.1 / 5.1 72D302G (GAC→GGC) PP_4881iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (38/28);  new base G (2/4);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

TCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTAC  >  NC_002947/5550564‑5550853
                                                                                                                                                   |                                                                                                                                              
tCACTATCTCCGGCGCAGGCGTAGTCAAGGCCAGAACACAAAAGCAACCACCGCAACACCTCCTTCACGGTGTTGCCGGCGACAGATACGAGCATGTGATCACGCAGACCACCGCCGAGTACCCACTGCACACGGGCATGGTTTCTGAcc                                                                                                                                              >  2:161905/1‑150 (MQ=255)
  aCTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACcgc                                                                                                                                            >  1:1401865/1‑150 (MQ=255)
    tATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGGGAAGAAGGCCCGCGCATAAAAGGCGATTTTTACACCGAGGACCCACTGCACAAGGGCATGGTTTCTGACCGCgg                                                                                                                                          >  1:1028427/1‑150 (MQ=255)
    tATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCgg                                                                                                                                          >  2:1191107/1‑150 (MQ=255)
      tctcCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGcc                                                                                                                                        <  1:1388769/150‑1 (MQ=255)
               cAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccat                                                                                                                               <  1:714104/150‑3 (MQ=255)
                          aaGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTg                                                                                                                    >  1:1370020/1‑150 (MQ=255)
                             gCCAGCAAACACCCGAAAGAAGCACAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATAACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTgcgc                                                                                                                 >  1:207519/1‑150 (MQ=255)
                               cAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGcgccc                                                                                                               >  2:422458/1‑150 (MQ=255)
                                 gCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTc                                                                                                                             >  1:321692/1‑134 (MQ=255)
                                 gCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTc                                                                                                                             <  2:321692/134‑1 (MQ=255)
                                    aaCACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGa                                                                                                          >  1:263165/1‑150 (MQ=255)
                                     acacCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGaa                                                                                                         >  1:1276401/1‑150 (MQ=255)
                                         ccGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtc                                                                                                     >  2:277619/1‑150 (MQ=255)
                                          cGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtct                                                                                                    <  1:422458/150‑1 (MQ=255)
                                          cGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtct                                                                                                    <  1:1191107/150‑1 (MQ=255)
                                           gAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                   <  2:1184540/150‑1 (MQ=255)
                                           gAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                   >  2:387656/1‑150 (MQ=255)
                                           gAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                   >  1:285978/1‑150 (MQ=255)
                                               gaagCCCAGGCCCTGCTCCGCTGGATGGCCACCGAAGAAGGCCAGCGCGGGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCa                                                                                               <  2:662810/150‑1 (MQ=255)
                                               gaagCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCa                                                                                               >  1:973990/1‑150 (MQ=255)
                                                    ccAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCTCCCGCCGAAGATCTCGCCAGCTGa                                                                                          >  1:453948/1‑150 (MQ=255)
                                                     cAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGAc                                                                                         >  2:885340/1‑150 (MQ=255)
                                                       ggCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAt                                                                                                      >  1:1258773/1‑135 (MQ=255)
                                                       ggCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAt                                                                                                      <  2:1258773/135‑1 (MQ=255)
                                                         cccTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTgg                                                                                     <  1:451572/150‑1 (MQ=255)
                                                            tGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatg                                                                                                                              >  1:360729/1‑103 (MQ=38)
                                                            tGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatg                                                                                                                              <  2:360729/106‑4 (MQ=38)
                                                                         ggCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGc                                                                     >  1:353545/1‑150 (MQ=255)
                                                                            cAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCAt                                                                  >  1:272691/1‑150 (MQ=255)
                                                                                   gaagGCCAGCGCGTGATCACCCAGACCACCGCCGAGTCCCCACTCCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGAGCAATGCCGAGGAAGCCATCGAACtt                                                           <  2:1028427/150‑1 (MQ=255)
                                                                                     agGCCAGCGCGTGATCGCCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGa                                                         <  2:1370020/150‑1 (MQ=255)
                                                                                       gCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGa                                                                        <  2:191926/133‑1 (MQ=255)
                                                                                       gCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGa                                                                        >  1:191926/1‑133 (MQ=255)
                                                                                         cAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGt                                                     <  2:608865/150‑1 (MQ=255)
                                                                                          aGCGCGTGATCACCCAGACCACCGCCGAGTAACCACTGCAAAAGGACATGGTTTATGACCGCGACCTGAAAGCGGTCGACAAGCTTCGCCCTCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTg                                                    >  2:1230675/1‑150 (MQ=255)
                                                                                            cgcgTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAg                                                  <  2:1401865/150‑1 (MQ=255)
                                                                                              cgTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGGCAGCGGCTTTAATCCGTTCCACGAGCTCCTTCTGTGGAGGATCTTTCCGGGTTTGCTGGGTAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGt                                                <  2:1230167/150‑1 (MQ=255)
                                                                                               gTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTc                                               >  1:1256625/1‑150 (MQ=255)
                                                                                                  ataaCCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGc                                            <  2:207519/147‑1 (MQ=255)
                                                                                                  aTCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGATCTCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTAACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGc                                            <  1:968569/150‑1 (MQ=255)
                                                                                                   acACCCAGACCACCGCCGAGTACACACTGGACAACGCCATGGTTTACAAACGCGGCCTGACACCGGTCGAGGTCATGAACCCGCGGAATGGTGCGCCGGCTGACCTTGGCAATGCCGAGGAAGCGATCGACCTTGAGCGTGAGTTCGTcc                                           >  2:519730/2‑150 (MQ=255)
                                                                                                     aCCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTg                                         >  1:143479/1‑150 (MQ=255)
                                                                                                      cccAAACAAACGCCGAGGACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGc                                        >  2:1348384/1‑150 (MQ=255)
                                                                                                       ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAg                                                                      <  2:1337630/119‑1 (MQ=255)
                                                                                                       ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAg                                                                      >  1:1337630/1‑119 (MQ=255)
                                                                                                       ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGct                                       >  2:711887/1‑150 (MQ=255)
                                                                                                       ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGct                                       >  2:1050398/1‑150 (MQ=255)
                                                                                                         aGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTCCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGctct                                     >  2:174438/1‑150 (MQ=255)
                                                                                                           accaccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtga                                   <  2:285978/150‑1 (MQ=255)
                                                                                                            ccaccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCATTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                          >  1:258452/1‑110 (MQ=255)
                                                                                                            ccaccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCATTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                          <  2:258452/110‑1 (MQ=255)
                                                                                                             caccGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCtgatg                                 >  1:92903/1‑150 (MQ=255)
                                                                                                               ccgccgAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGAc                               <  2:272691/150‑1 (MQ=255)
                                                                                                                   cgAGTACCCACTGCACAAGGGCATGGTTTGTGACCGCGGCCTGACTCCGTTCGACGACGCGCTCGCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTgcc                           <  1:526499/150‑1 (MQ=255)
                                                                                                                   cgAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTgcc                           <  1:1348384/150‑1 (MQ=255)
                                                                                                                      gTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTgccgcc                        >  1:551913/1‑150 (MQ=255)
                                                                                                                          ccACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGt                    <  2:1256625/150‑1 (MQ=255)
                                                                                                                           cACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                   >  2:657256/1‑70 (MQ=255)
                                                                                                                           cACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                   <  1:657256/70‑1 (MQ=255)
                                                                                                                                cacaAGGGCATGGTTTCTGACCGCGGCCTGAAACGGTTCGACGAGCTGCGCTCGCCGAAGATCTTGCCAGCTGACGTGAGCGATGCCGAGGAAACGAACGAGCTTGAGCGTGAGGAAGACCTGCTCTGATGGCTGCTGCACTGTCTGAAc              >  2:203064/1‑150 (MQ=255)
                                                                                                                                       gCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATg                                                                               <  1:758905/78‑1 (MQ=255)
                                                                                                                                       gCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATg                                                                               >  2:758905/1‑78 (MQ=255)
                                                                                                                                        cATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCg      <  1:277619/150‑1 (MQ=255)
                                                                                                                                        aacggTTTCAGGCCGCGGTCAGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGt          <  1:546065/143‑1 (MQ=255)
                                                                                                                                        aacggTTTCAGGCCGCGGTCAGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGt          >  2:546065/4‑146 (MQ=255)
                                                                                                                                         aTGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCgg     >  1:601585/1‑150 (MQ=255)
                                                                                                                                            gTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTAc  <  1:711887/150‑1 (MQ=255)
                                                                                                                                             tttCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                          >  2:1245908/1‑77 (MQ=255)
                                                                                                                                             tttCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                          <  1:1245908/77‑1 (MQ=255)
                                                                                                                                                 tGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGa                                                   <  2:48183/96‑1 (MQ=255)
                                                                                                                                                 tGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGa                                                   >  1:48183/1‑96 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTAC  >  NC_002947/5550564‑5550853

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: