Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,224,539 G→A 6.5% intergenic (‑607/+214) PP_3703 ← / ← PP_3704 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,224,5390GA 6.5% 178.0 / 2.4 61intergenic (‑607/+214)PP_3703/PP_3704hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base A (2/2);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

ACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAG  >  NC_002947/4224390‑4224675
                                                                                                                                                     |                                                                                                                                        
aCAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGCACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCcagcag                                                                                                                                          >  2:664269/1‑150 (MQ=255)
 cAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGt                                                                                                                                         <  2:1278154/150‑1 (MQ=255)
 cAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGt                                                                                                                                         <  1:301975/150‑1 (MQ=255)
     tcgtcgCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCAc                                                                                                                                     >  1:255862/1‑150 (MQ=255)
           cccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCtgt                                                                                                                                 <  1:1175279/148‑1 (MQ=255)
           cccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCtgt                                                                                                                                 >  2:1175279/1‑148 (MQ=255)
            ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCa                                                                                                                              >  2:263856/1‑150 (MQ=255)
            ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCa                                                                                                                              >  1:1250927/1‑150 (MQ=255)
             cTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCAGCTGTGCAt                                                                                                                             >  1:601598/1‑150 (MQ=255)
              tATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCAtt                                                                                                                            <  2:844527/150‑1 (MQ=255)
               aTCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTg                                                                                                                           >  2:686427/1‑150 (MQ=255)
                 cAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCt                                                                                                                         <  1:686427/150‑1 (MQ=255)
                          cGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCAc                                                                                                                >  1:231917/1‑150 (MQ=255)
                          cGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCAc                                                                                                                >  1:232010/1‑150 (MQ=255)
                                 gTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                         <  1:1113177/150‑1 (MQ=255)
                                 gTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                         <  1:614356/150‑1 (MQ=255)
                                    ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                      >  2:808587/1‑150 (MQ=255)
                                        aGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                  <  2:231917/150‑1 (MQ=255)
                                        aGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                  <  2:232010/150‑1 (MQ=255)
                                          ccTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAGCCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGa                                                                                                <  2:1339748/150‑1 (MQ=255)
                                            tGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGc                                                                                              <  1:263856/150‑1 (MQ=255)
                                             gggagcggGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCt                                                                                             >  2:452363/5‑150 (MQ=255)
                                                   ggggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGAc                                                                                       >  1:23880/1‑150 (MQ=255)
                                                    gggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACa                                                                                      <  1:224014/150‑1 (MQ=255)
                                                    gggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACa                                                                                      <  1:911780/150‑1 (MQ=255)
                                                    gggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACa                                                                                      <  2:673286/150‑1 (MQ=255)
                                                           acCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                               <  2:23880/150‑1 (MQ=255)
                                                            ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTggg                                                                              >  2:757811/1‑150 (MQ=255)
                                                              tCGCTCCCGTAAGAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGGGTACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTa                                                                            <  1:664269/150‑1 (MQ=255)
                                                                gCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCAcc                                                                                                                                    >  1:687168/1‑92 (MQ=255)
                                                                gCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCAcc                                                                                                                                    <  2:687168/92‑1 (MQ=255)
                                                                gCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATc                                                                          >  2:1005928/1‑150 (MQ=255)
                                                                  tCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATcgc                                                                        >  1:210252/1‑150 (MQ=255)
                                                                     cGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGcc                                                                     <  2:439472/150‑1 (MQ=255)
                                                                      gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGt                                                                             >  2:716071/1‑141 (MQ=255)
                                                                      gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGt                                                                             <  1:716071/141‑1 (MQ=255)
                                                                      gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGccc                                                                    >  1:868247/1‑150 (MQ=255)
                                                                            aaGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGt                                                              >  1:1311242/1‑150 (MQ=255)
                                                                                cGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCg                                                          <  2:255862/150‑1 (MQ=255)
                                                                                  gCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGga                                                        <  2:868247/150‑1 (MQ=255)
                                                                                  gCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGga                                                        <  1:1005928/150‑1 (MQ=255)
                                                                                         aCGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGCTGATa                                                 <  2:1091509/150‑1 (MQ=255)
                                                                                          cGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATAc                                                <  2:601598/150‑1 (MQ=255)
                                                                                          cGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATAc                                                <  2:210252/150‑1 (MQ=255)
                                                                                           gTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACa                                               >  2:1049983/1‑150 (MQ=255)
                                                                                               tAAGCGGCGCGACGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAAc                                           <  2:969154/150‑1 (MQ=255)
                                                                                                     gAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                      <  2:1198862/85‑1 (MQ=255)
                                                                                                     gAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                      >  1:1198862/1‑85 (MQ=255)
                                                                                                             cGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCa                             <  1:757811/150‑1 (MQ=255)
                                                                                                              gggAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCAt                            >  2:237992/1‑150 (MQ=255)
                                                                                                               ggAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATg                           <  1:452363/150‑1 (MQ=255)
                                                                                                                       ccGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTg                   >  1:1131041/1‑150 (MQ=255)
                                                                                                                              aaaTAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgc            <  2:1311242/150‑1 (MQ=255)
                                                                                                                              aaaTAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgc            <  1:237992/150‑1 (MQ=255)
                                                                                                                                aTAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgcgc          >  2:506605/1‑150 (MQ=255)
                                                                                                                                 tAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGcc         >  2:616123/1‑150 (MQ=255)
                                                                                                                                        aGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCaaggaag  <  1:506605/150‑1 (MQ=255)
                                                                                                                                        aGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCaaggaag  <  1:616123/150‑1 (MQ=255)
                                                                                                                                             aaccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCAt                            >  1:578595/3‑119 (MQ=255)
                                                                                                                                             aaccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCAt                            <  2:578595/117‑1 (MQ=255)
                                                                                                                                             aaccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTg                   >  1:842570/3‑128 (MQ=255)
                                                                                                                                             aaccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTg                   <  2:842570/126‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                        
ACAGATCGTCGCCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAG  >  NC_002947/4224390‑4224675

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