Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,023,215 G→A 20.6% intergenic (+74/+70) PP_5273 → / ← PP_5274 oxidoreductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,023,2150GA20.6% 131.6 / 30.5 68intergenic (+74/+70)PP_5273/PP_5274oxidoreductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/32);  new base A (7/7);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 5.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CCTGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAACTCGTCATCG  >  NC_002947/6023067‑6023363
                                                                                                                                                    |                                                                                                                                                    
ccTGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGt                                                                                                                                                     <  1:1614904/150‑1 (MQ=255)
ccTGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGt                                                                                                                                                     <  2:38683/150‑1 (MQ=255)
ccTGGCGGCAAACACCTGCGCTGGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGt                                                                                                                                                     <  2:93187/150‑1 (MQ=255)
  tGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGc                                                                                                                                                   >  1:1130823/1‑150 (MQ=255)
  tGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGc                                                                                                                                                   <  2:1528469/150‑1 (MQ=255)
         aaGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTcc                                                                                                                                            <  1:1565019/150‑1 (MQ=255)
         aaGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGACGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTcc                                                                                                                                            <  1:1319723/150‑1 (MQ=255)
           gCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTCCca                                                                                                                                          <  1:668036/150‑1 (MQ=255)
              ccTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCGGCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAg                                                                                                                                       <  1:442980/150‑1 (MQ=255)
               cTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAgg                                                                                                                                      <  2:1130823/150‑1 (MQ=255)
               cTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAgg                                                                                                                                      >  1:168893/1‑150 (MQ=255)
                       tccACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGc                                                                                                                              >  2:683973/2‑150 (MQ=255)
                       gCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGc                                                                                                                              >  2:1311078/1‑150 (MQ=255)
                       gCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTCCcaca                                                                                                                                        <  2:513420/140‑1 (MQ=255)
                       gCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTCCcaca                                                                                                                                        >  1:513420/1‑140 (MQ=255)
                         cACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCaa                                                                                                                            <  2:168893/150‑1 (MQ=255)
                         cACTGCTGGCGCTCGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGGGGGCGTGCACGCTCCCACAGGGTCGGTGCaa                                                                                                                            >  2:320860/1‑150 (MQ=255)
                         cACGGCTGGCTCGGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCaa                                                                                                                            <  2:1583168/150‑1 (MQ=255)
                           cggctgGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGt                                                                                                                          >  2:646048/3‑150 (MQ=255)
                               tgGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCcaca                                                                                                                                        >  2:1326097/1‑132 (MQ=255)
                               tgGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCcaca                                                                                                                                        <  1:1326097/132‑1 (MQ=255)
                                      gggAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTgg                                                                                                               >  1:273730/1‑150 (MQ=255)
                                       ggAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTggg                                                                                                              <  1:409501/150‑1 (MQ=255)
                                              tgtgGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaacaac                                                                                                       <  2:927934/150‑1 (MQ=255)
                                                tgGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAg                                                                                                     <  2:134912/150‑1 (MQ=255)
                                                 ggCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAgg                                                                                                    <  2:273730/150‑1 (MQ=255)
                                                            aaaGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaacaac                                                                                                       >  2:741061/1‑136 (MQ=255)
                                                            aaaGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaacaac                                                                                                       <  1:741061/136‑1 (MQ=255)
                                                            aaaGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaacaa                                                                                                        <  1:1478338/135‑1 (MQ=255)
                                                            aaaGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaacaa                                                                                                        >  2:1478338/1‑135 (MQ=255)
                                                            aaaCAGCTGACTTGAAAGTCCCCATCCTTTTGGCAGCAACTCTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGTAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGa                                                                                         <  1:391771/150‑1 (MQ=255)
                                                               tagCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCa                                                                                      <  1:683973/149‑1 (MQ=255)
                                                                    gACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTcc                                                                                                                                            >  2:59925/1‑91 (MQ=255)
                                                                    gACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTcc                                                                                                                                            <  1:59925/91‑1 (MQ=255)
                                                                     aCTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaa                                                                                                           <  2:234628/123‑1 (MQ=255)
                                                                     aCTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaa                                                                                                           >  1:234628/1‑123 (MQ=255)
                                                                     aCTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCgtgt                                                                                <  1:168873/150‑1 (MQ=255)
                                                                     aCTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCgtgt                                                                                <  2:1189114/150‑1 (MQ=255)
                                                                      cTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGa                                                                                         <  2:1330216/140‑1 (MQ=255)
                                                                      cTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGATTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGa                                                                                         >  1:1330216/1‑140 (MQ=255)
                                                                       ttGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAgg                                                                                                                                      >  1:1267758/1‑94 (MQ=255)
                                                                       ttGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAgg                                                                                                                                      <  2:1267758/94‑1 (MQ=255)
                                                                       ttGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCa                                                                              <  1:1220715/150‑1 (MQ=255)
                                                                        tGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTTAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTc                                                                                                                                             >  2:1231056/1‑86 (MQ=255)
                                                                        tGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTTAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTc                                                                                                                                             <  1:1231056/86‑1 (MQ=255)
                                                                         gAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGt                                                                            <  1:320860/150‑1 (MQ=255)
                                                                            aGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCt                                                                         >  1:75777/1‑150 (MQ=255)
                                                                               ccccATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTATGAGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTcc                                                                      >  1:349154/1‑150 (MQ=255)
                                                                                    tCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAg                                                                 <  1:1311078/150‑1 (MQ=255)
                                                                                                            ttAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGcc                                         <  1:336858/150‑1 (MQ=255)
                                                                                                            ttAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGACGGTGCCAGTTTCAAGTCTGGGGAACAACAGGCCCTGCAGCGACCATTGTGTCAGTCCATCCAGAAGGGTTTCCGCCCTTCGATAAGGGcc                                         >  2:129929/1‑150 (MQ=255)
                                                                                                                   cTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCACGTTTCAAGTCTGGGGCACCACATGCCCACCGGCGGCCAGCAGGGCAGGCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTccc                                  >  2:1419515/1‑150 (MQ=255)
                                                                                                                     aaaaGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAg                                >  1:1338311/1‑150 (MQ=255)
                                                                                                                                   ccccTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCtt                                                    >  2:1213784/1‑116 (MQ=255)
                                                                                                                                   ccccTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCtt                                                    <  1:1213784/116‑1 (MQ=255)
                                                                                                                                    cccTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaac                                                                                                          >  2:943574/1‑61 (MQ=255)
                                                                                                                                    cccTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaac                                                                                                          <  1:943574/61‑1 (MQ=255)
                                                                                                                                    cccTGTGGGAGCGGGCACGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAg                                                                 >  2:912973/1‑102 (MQ=255)
                                                                                                                                    cccTGTGGGAGCGGGCACGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAg                                                                 <  1:912973/102‑1 (MQ=255)
                                                                                                                                     ccTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATcgc                <  2:75777/150‑1 (MQ=255)
                                                                                                                                       tgtgGGAGCGGGCACGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCgcgc              >  2:520188/1‑150 (MQ=255)
                                                                                                                                        gtgGGAGCGGGCACGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGaa                                                                  >  1:696817/1‑97 (MQ=255)
                                                                                                                                        gtgGGAGCGGGCACGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGaa                                                                  <  2:696817/97‑1 (MQ=255)
                                                                                                                                         tgGGAGCGGGCACGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCaa            >  1:1489820/1‑150 (MQ=255)
                                                                                                                                          gggAGCGGGCACGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCATTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAAc           >  1:1173897/1‑150 (MQ=255)
                                                                                                                                            gAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAACTc         >  2:1491423/1‑150 (MQ=255)
                                                                                                                                               cGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTTTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAACTCGTc      <  2:952115/150‑1 (MQ=255)
                                                                                                                                                   cGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAACTCGTCATCg  <  1:646048/150‑1 (MQ=255)
                                                                                                                                                    gTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATc                                                                                    >  2:1632355/1‑67 (MQ=255)
                                                                                                                                                    gTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATc                                                                                    <  1:1632355/67‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAACTCGTCATCG  >  NC_002947/6023067‑6023363

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