Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,642,343 G→A 13.1% intergenic (‑129/+23) PP_1442 ← / ← lon‑I hypothetical protein/DNA‑binding ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,642,3430GA13.1% 205.6 / 17.7 82intergenic (‑129/+23)PP_1442/lon‑Ihypothetical protein/DNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base G (36/35);  new base A (5/6);  total (41/41)
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.91e-01

AGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCAGTTCGAAAATC  >  NC_002947/1642203‑1642492
                                                                                                                                             |                                                                                                                                                     
aggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCCTTCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGaaa                                                                                                                                              <  1:1530294/150‑2 (MQ=18)
 ggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAaga                                                                                                                                             <  1:1062217/150‑1 (MQ=255)
 ggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAaga                                                                                                                                             <  1:1173095/150‑1 (MQ=255)
      ggTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCACCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCg                                                                                                                                        <  2:165378/150‑1 (MQ=255)
       gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGaa                                                                                                                                               <  2:434161/142‑1 (MQ=255)
       gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGaa                                                                                                                                               >  1:434161/1‑142 (MQ=255)
          aTGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGc                                                                                                                                    <  2:686038/150‑1 (MQ=255)
                  tCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAaga                                                                                                                                             <  2:280819/133‑1 (MQ=255)
                  tCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAaga                                                                                                                                             >  1:280819/1‑133 (MQ=255)
                   ccGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGaa                                                                                                                           >  2:997933/1‑150 (MQ=255)
                     gAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                                         <  1:25801/150‑1 (MQ=255)
                        ggCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCa                                                                                                                      >  1:1712616/1‑150 (MQ=255)
                           tGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACtt                                                                                                                   >  1:230207/1‑150 (MQ=255)
                             gCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTc                                                                                                                 >  1:777629/1‑150 (MQ=255)
                               gATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGc                                                                                                               >  1:1415092/1‑150 (MQ=255)
                                aTGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGcc                                                                                                              <  2:165044/150‑1 (MQ=255)
                                aTGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGcc                                                                                                              <  2:1692956/150‑1 (MQ=255)
                                    ccTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGt                                                                                                          >  1:669076/1‑150 (MQ=255)
                                    ccTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGt                                                                                                          <  2:230207/150‑1 (MQ=255)
                                     cTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTc                                                                                                         >  1:1431739/1‑150 (MQ=255)
                                      tGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTATGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTACGTGTTGACCCTCGAAGAGGTCGTTGCAGGCTAGAACAGCACTTTCGCCACGTCg                                                                                                        >  2:990211/1‑150 (MQ=255)
                                          gATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGCTATGGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCaa                                                                                                    >  1:1453224/1‑150 (MQ=255)
                                           aTACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCaaa                                                                                                   <  1:997933/150‑1 (MQ=255)
                                               ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                                               >  2:1145843/1‑150 (MQ=255)
                                               ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                                               >  1:1033838/1‑150 (MQ=255)
                                                   tgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTgg                                                                                           <  2:792547/150‑1 (MQ=255)
                                                      tCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTTGCGa                                                                                        >  1:1767474/1‑150 (MQ=255)
                                                        gAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGaaa                                                                                      <  2:923915/150‑1 (MQ=255)
                                                        gAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGCACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCAc                                                                                                            >  2:1145309/1‑128 (MQ=255)
                                                        gAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGCACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCAc                                                                                                            <  1:1145309/128‑1 (MQ=255)
                                                         aGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTACGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAAt                                                                                     >  2:694559/1‑150 (MQ=255)
                                                              aTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACt                                                                                >  2:1491539/1‑150 (MQ=255)
                                                              aTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCACGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACt                                                                                <  1:959000/150‑1 (MQ=255)
                                                                ttGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTgg                                                                              >  1:608827/1‑149 (MQ=255)
                                                                   tgctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAg                                                                           <  2:489296/150‑1 (MQ=255)
                                                                     ctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGc                                                                         <  2:1415092/150‑1 (MQ=255)
                                                                      tgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAGTGTCAGGcc                                                                        >  1:877177/1‑150 (MQ=255)
                                                                      tgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGcc                                                                        <  1:1356384/150‑1 (MQ=255)
                                                                        ccTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCa                                                                            >  2:1014197/1‑144 (MQ=255)
                                                                        ccTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCa                                                                            <  1:1014197/144‑1 (MQ=255)
                                                                          tgtgCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTcg                                                                    >  1:1301589/1‑150 (MQ=255)
                                                                          tgtgCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTcg                                                                    >  2:1307635/1‑150 (MQ=255)
                                                                            tgCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTcgcg                                                                  >  2:1596652/1‑150 (MQ=255)
                                                                               cGGTCCTTTCGCGAGTAAGCCCGCGCCTGCATGTGTAGCGTACCCGTTCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACGGTCAGGCCTTCGCGCAg                                                               <  2:1834967/150‑1 (MQ=255)
                                                                                gccctcttCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAgg                                                              >  1:660738/7‑150 (MQ=17)
                                                                                    tcttCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGt                                                          >  1:325993/3‑150 (MQ=18)
                                                                                    cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGt                                                          >  2:293978/1‑150 (MQ=255)
                                                                                     tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc                                                         <  2:1640170/150‑1 (MQ=255)
                                                                                     tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc                                                         <  1:1307635/150‑1 (MQ=255)
                                                                                      ttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTcc                                                        <  2:473920/150‑1 (MQ=255)
                                                                                         gcgAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGc                                                     >  2:380049/1‑150 (MQ=255)
                                                                                           gAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAg                                                   >  1:1439017/1‑150 (MQ=255)
                                                                                           gAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGAAAGAGCTGCAAAGCGCTTGGCGCAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAg                                                   >  1:498077/1‑150 (MQ=255)
                                                                                             gTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGACATGCCTTCGCGCAGGTAGTCCGGCAGtt                                                 <  1:990211/150‑1 (MQ=255)
                                                                                              tAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCGAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCAGGCAGttc                                                >  2:1039123/1‑150 (MQ=255)
                                                                                                 gCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                                             <  2:1301589/150‑1 (MQ=255)
                                                                                                 gCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                                             <  2:1767474/150‑1 (MQ=255)
                                                                                                 gCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                                             >  2:1305898/1‑150 (MQ=25)
                                                                                                  cccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCa                                                                            >  1:634854/1‑118 (MQ=255)
                                                                                                  cccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCa                                                                            <  2:634854/118‑1 (MQ=255)
                                                                                                   ccGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGa                                           >  1:154075/1‑150 (MQ=25)
                                                                                                    cgcgCCTGCAGGGGTAGCGTACCCGTGCAGGTACGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                                          <  1:694559/150‑1 (MQ=255)
                                                                                                         cTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTcgc                                     >  2:1173514/1‑150 (MQ=255)
                                                                                                          tGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCa                                                                <  2:890133/122‑1 (MQ=255)
                                                                                                          tGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGGCAGGCCTTCGCGCa                                                                >  1:890133/1‑122 (MQ=255)
                                                                                                            cAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTcgccgc                                  >  2:1146795/1‑150 (MQ=255)
                                                                                                                 gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg                             <  1:1491539/150‑1 (MQ=255)
                                                                                                                 gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg                             <  1:1145843/150‑1 (MQ=255)
                                                                                                                    gCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGc                          >  2:28669/1‑150 (MQ=255)
                                                                                                                           ccGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGc                    >  1:1342504/1‑150 (MQ=255)
                                                                                                                            cctGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCa                   <  2:154075/148‑1 (MQ=255)
                                                                                                                             gTGCAGGTTCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGt                               <  2:964104/137‑1 (MQ=255)
                                                                                                                             gTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGt                               >  1:964104/1‑137 (MQ=255)
                                                                                                                             ctGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAg                  <  1:1305898/149‑1 (MQ=255)
                                                                                                                               gCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGa                <  1:828175/150‑1 (MQ=255)
                                                                                                                                cAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGAt               <  2:325993/150‑1 (MQ=255)
                                                                                                                                 atgTGCGGGCTTGCCCGAGGGGAGGGAGGTGAAGGCTAGCTCGGCATTGTCTCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATc              <  2:1453224/148‑1 (MQ=255)
                                                                                                                                    tGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCAGt           <  2:1033838/150‑1 (MQ=255)
                                                                                                                                     gCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCAGtt          <  1:1183851/150‑1 (MQ=255)
                                                                                                                                     gCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCAGtt          <  1:1184102/150‑1 (MQ=255)
                                                                                                                                      cGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCAGTTc         >  1:1336478/1‑150 (MQ=255)
                                                                                                                                           ttGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCAGTTCGaaaa    >  1:1079698/1‑150 (MQ=255)
                                                                                                                                             gCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCAGTTCGAAAAtc  >  2:87988/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
AGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCAGTTCGAAAATC  >  NC_002947/1642203‑1642492

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: