Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,550,718 C→T 6.8% G304G (GGC→GGT)  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,7180CT 6.8% 208.2 / 4.7 73G304G (GGC→GGT) PP_4881iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base T (2/3);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

TCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAGC  >  NC_002947/5550570‑5550867
                                                                                                                                                    |                                                                                                                                                     
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:1191107/150‑1 (MQ=255)
                                    cGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtct                                                                                                                  <  1:422458/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                                                                                                                                            <  2:360729/106‑4 (MQ=38)
                                                      tGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCAGGCCGCGGTCAGAAAccatg                                                                                                                                            >  1:360729/1‑103 (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                                                                                      >  1:191926/1‑133 (MQ=255)
                                                                                 gCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGa                                                                                      <  2:191926/133‑1 (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                                                                                    >  1:1337630/1‑119 (MQ=255)
                                                                                                 ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAg                                                                                    <  2:1337630/119‑1 (MQ=255)
                                                                                                 ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGct                                                     >  2:1050398/1‑150 (MQ=255)
                                                                                                 ccAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGct                                                     >  2:711887/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                                                                                                                 <  1:657256/70‑1 (MQ=255)
                                                                                                                     cACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGAtctc                                                                                                                 >  2:657256/1‑70 (MQ=255)
                                                                                                                          cacaAGGGCATGGTTTCTGACCGCGGCCTGAAACGGTTCGACGAGCTGCGCTCGCCGAAGATCTTGCCAGCTGACGTGAGCGATGCCGAGGAAACGAACGAGCTTGAGCGTGAGGAAGACCTGCTCTGATGGCTGCTGCACTGTCTGAAc                            >  2:203064/1‑150 (MQ=255)
                                                                                                                                 gCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATg                                                                                             >  2:758905/1‑78 (MQ=255)
                                                                                                                                 gCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATg                                                                                             <  1:758905/78‑1 (MQ=255)
                                                                                                                                  cATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCg                    <  1:277619/150‑1 (MQ=255)
                                                                                                                                  aacggTTTCAGGCCGCGGTCAGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGt                        >  2:546065/4‑146 (MQ=255)
                                                                                                                                  aacggTTTCAGGCCGCGGTCAGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGt                        <  1:546065/143‑1 (MQ=255)
                                                                                                                                   aTGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCgg                   >  1:601585/1‑150 (MQ=255)
                                                                                                                                      gTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTAc                <  1:711887/150‑1 (MQ=255)
                                                                                                                                       tttCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                                        <  1:1245908/77‑1 (MQ=255)
                                                                                                                                       tttCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAg                                                                                        >  2:1245908/1‑77 (MQ=255)
                                                                                                                                           tggccGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTc           >  1:1422821/4‑150 (MQ=255)
                                                                                                                                           tGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGa                                                                 >  1:48183/1‑96 (MQ=255)
                                                                                                                                           tGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGa                                                                 <  2:48183/96‑1 (MQ=255)
                                                                                                                                              ccGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTa        >  1:72012/1‑150 (MQ=255)
                                                                                                                                               cgcgGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTAc       >  1:1297492/1‑150 (MQ=255)
                                                                                                                                                  ggCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGa    <  2:353545/150‑1 (MQ=255)
                                                                                                                                                  ggCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGa    <  1:1050398/150‑1 (MQ=255)
                                                                                                                                                   gCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAg   >  2:268020/1‑150 (MQ=255)
                                                                                                                                                    ccTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGGCCGAACCGGTGCCGGTACGCTTCGTACCGAGc  >  1:205123/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAGCTGCGCCCGCCGAAGATCTCGCCAGCTGACCTGGGCAATGCCGAGGAAGCCATCGAACTTGAGCGTGAGGTCGGCCTGCTCTGATGACTGCCGCCCTGTCCGAACCGGTGCCGGTACGCTTCGTACCGAGC  >  NC_002947/5550570‑5550867

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