Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,928,771 G→A 7.6% R112C (CGC→TGC)  fbpA ← iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,928,7710GA 7.6% 185.5 / 5.0 66R112C (CGC→TGC) fbpAiron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (40/21);  new base A (2/3);  total (42/24)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CTTGGCCGTGCGCAGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATGCGCTGCAT  >  NC_002947/5928622‑5928918
                                                                                                                                                     |                                                                                                                                                   
cttGGCCGTGCGCAGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGcgcg                                                                                                                                                     <  1:412215/150‑1 (MQ=255)
 ttGGCCGTGCGCAGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCgg                                                                                                                                                    <  2:686206/150‑1 (MQ=255)
   ggCCGTGCGCAGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTGCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTa                                                                                                                                                  >  2:497427/1‑150 (MQ=255)
             aGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCggg                                                                                                                                        >  2:1306347/1‑150 (MQ=255)
             aGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCggg                                                                                                                                        >  1:670902/1‑150 (MQ=255)
               ccacaGCCGCCCTTCCCACCGTTGGCCACCCAGGCCTTCGTAGGTGCTCAACTCCTCGGGTTTGCCCCGCTCGGTGGAGTAGATGATGGTGCGGGCCCGCAGGCTCAGGCCGGTGCAGTCATGGGACGAGGCGCAGTACCGCGcctcgtc                                                                                                                                      <  1:665290/149‑6 (MQ=255)
                cacaGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATg                                                                                                                                     <  2:1366916/150‑1 (MQ=255)
                cacaGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCAGTACCGCGcctcgtcc                                                                                                                                     <  2:1077055/150‑7 (MQ=255)
                   aGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCGAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTc                                                                                                                                  >  1:1143757/1‑150 (MQ=255)
                              ccACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGAt                                                                                                                                      >  2:538286/1‑135 (MQ=255)
                              ccACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGAt                                                                                                                                      <  1:538286/135‑1 (MQ=255)
                               cACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACg                                                                                                                      >  2:490456/1‑150 (MQ=255)
                                aCTTTTTGCTAGCCAGGGCTTCGTAGGTGCTCAACTCCTAGGGTTTGACCCGTCAGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGt                                                                                                                     <  1:162983/150‑1 (MQ=255)
                                aCTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGATCGATGACGt                                                                                                                     >  2:790303/1‑150 (MQ=255)
                                   gTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCgg                                                                                                                  <  2:1143757/150‑1 (MQ=255)
                                          cccgtctgCTTCGTAGGGGCTCATCTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGTAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGAt                                                                                                           <  1:871187/143‑1 (MQ=255)
                                            ccAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCg                                                                                                         >  2:1163305/1‑150 (MQ=255)
                                            ccAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGAAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTCCGCCGGGATCTTCCGGTCGATGACGTCGGACTCGATCg                                                                                                         >  2:109847/1‑150 (MQ=255)
                                              aGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCAGTAccg                                                                                                                                               <  1:622976/110‑3 (MQ=255)
                                              aGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCAGTAccg                                                                                                                                               >  2:622976/1‑108 (MQ=255)
                                                  cTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCTTGGTGTTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATcagatcgg                                                                                                   >  2:871187/1‑143 (MQ=255)
                                                   ttCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCAGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAg                                                                                                  >  2:910925/1‑150 (MQ=255)
                                                          gTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGccc                                                                                           >  2:1298651/1‑150 (MQ=255)
                                                                aaCTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGc                                                                                     >  2:1002126/1‑150 (MQ=255)
                                                                     ctcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGAt                                                                                                                          >  2:579080/1‑108 (MQ=255)
                                                                     ctcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGAt                                                                                                                          <  1:579080/108‑1 (MQ=255)
                                                                     ctcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGc                                                                                >  2:261584/1‑150 (MQ=255)
                                                                     ctcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGc                                                                                >  2:7753/1‑150 (MQ=255)
                                                                     ctcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCAAGATCGGTCTGGGCGATGACGTCGGACTCGATCGGCGGCAGAATGCCCATCTGCTCGGc                                                                                >  1:728344/1‑150 (MQ=255)
                                                                     ctcGGGTTTGACCAGCTAGGCGGAGAAGACCATGGTGCGGGCGCGCTTGGGCGGGCCGGTTGAGGCGTGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGc                                                                                >  2:375166/1‑150 (MQ=255)
                                                                      tcGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGcc                                                                               >  1:279404/1‑150 (MQ=255)
                                                                       cGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCt                                                                              <  1:1163305/150‑1 (MQ=255)
                                                                        gggTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTg                                                                             <  1:1298651/150‑1 (MQ=255)
                                                                         ggTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGc                                                                            <  2:1045733/150‑1 (MQ=255)
                                                                          gTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGCTGTTCTGCTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGcc                                                                           >  2:226723/1‑150 (MQ=255)
                                                                                   gCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGtt                                                                                                                                   <  2:393272/85‑1 (MQ=255)
                                                                                   gCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGtt                                                                                                                                   >  1:393272/1‑85 (MQ=255)
                                                                                   gCTCGGTGGAGTAGATGATCGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGc                                                                  <  1:1141756/150‑1 (MQ=255)
                                                                                        gTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGGTGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCa                                                             >  2:1174673/1‑150 (MQ=255)
                                                                                        gTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCa                                                             <  1:7753/150‑1 (MQ=255)
                                                                                          ggAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATNGGCTGCNGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATc                                                           <  1:777105/150‑1 (MQ=255)
                                                                                             gTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGAc                                                        <  1:1002126/150‑1 (MQ=255)
                                                                                             gTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGAc                                                        <  2:878327/150‑1 (MQ=255)
                                                                                             gTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGAc                                                        <  1:198284/150‑1 (MQ=255)
                                                                                             gTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGAc                                                        >  1:364066/1‑150 (MQ=255)
                                                                                                 atgatgGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTc                                                    <  2:364066/150‑1 (MQ=255)
                                                                                                   gatgGTGACGGCGCGCGAGCTGAGGCCGGCCCAGGCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCCACGGTcag                                                  >  2:524879/1‑150 (MQ=255)
                                                                                                       gTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCt                                                                                                                                 <  2:982865/67‑1 (MQ=255)
                                                                                                       gTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCt                                                                                                                                 >  1:982865/1‑67 (MQ=255)
                                                                                                       gTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCt                                                                                                                                 >  1:272680/1‑67 (MQ=255)
                                                                                                       gTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCt                                                                                                                                 <  2:272680/67‑1 (MQ=255)
                                                                                                        tGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTcagcagca                                             >  2:1212600/1‑150 (MQ=255)
                                                                                                        tGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGAGCGGGACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTcagcagca                                             >  1:684655/1‑150 (MQ=255)
                                                                                                                     gCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGcc                                >  2:628367/1‑150 (MQ=255)
                                                                                                                     gCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGACCAGCTAAACGATGTGCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGcc                                >  2:984590/1‑150 (MQ=255)
                                                                                                                      cTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCg                               >  2:42556/1‑150 (MQ=255)
                                                                                                                       tCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGt                              >  1:348880/1‑150 (MQ=255)
                                                                                                                           gCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATc                          >  2:1001095/1‑150 (MQ=255)
                                                                                                                                tCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCt                     >  1:819207/1‑150 (MQ=255)
                                                                                                                                   aGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCgg                  >  2:123213/1‑150 (MQ=255)
                                                                                                                                   aGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCgg                  >  1:968402/1‑150 (MQ=255)
                                                                                                                                       aTGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCtt              >  2:998729/1‑150 (MQ=255)
                                                                                                                                        tGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTg             >  1:559245/1‑150 (MQ=255)
                                                                                                                                          ggACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGAt           >  1:612847/1‑150 (MQ=255)
                                                                                                                                          ggACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGAt           <  2:968402/150‑1 (MQ=255)
                                                                                                                                           gACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGACTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGt                                                     <  1:771558/107‑1 (MQ=255)
                                                                                                                                           gACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGACTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGt                                                     >  2:771558/1‑107 (MQ=255)
                                                                                                                                            aCGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATgc         >  1:788488/1‑150 (MQ=255)
                                                                                                                                                   gcgGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATGCGCTGCAt  >  1:1181063/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTTGGCCGTGCGCAGGCACAGACGGCCTTCCCACTGTTTGTCAGCCAGGGCTTCGTAGGTGCTCAACTCCTCGGGTTTGACCCGCTCGGTGGAGTAGATGATGGTGCGGGCGCGCAGGCTCAGGCCGGTCCAGTCATGGGACGAGGCGCGGTACTGCGCCGGGATGTTCTGGTCGATGACGTCGGACTCGATCGGCTGCAGAATGCCCATCTGCTCGGCCTGCCACAGGTTGCCTGCATCGACGGTCAGCAGCAGGTCGGCCACGCCGTTATCGCCCTCGGCCTTGATGCGCTGCAT  >  NC_002947/5928622‑5928918

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: