Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 564,380:1 +A 13.8% intergenic (+29/‑144) rplQ → / → katA 50S ribosomal protein L17/catalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947564,3801.A13.8% 357.7 / 57.6 109intergenic (+29/‑144)rplQ/katA50S ribosomal protein L17/catalase
Reads supporting (aligned to +/- strand):  ref base . (39/55);  new base A (6/9);  total (45/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 112 total aligned reads displayed.

CCGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAG  >  NC_002947/564231‑564528
                                                                                                                                                      |                                                                                                                                                    
ccGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGc                                                                                                                                                       >  2:1635591/1‑150 (MQ=255)
ccGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGc                                                                                                                                                       >  2:1197971/1‑150 (MQ=255)
ccGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGc                                                                                                                                                       >  1:314017/1‑150 (MQ=255)
   tCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTa                                                                                                                                                   <  2:2363311/150‑1 (MQ=255)
         cggcTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGccc                                                                                                                                             <  1:1246221/150‑1 (MQ=255)
          ggcTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCg                                                                                                                                            <  2:832553/150‑1 (MQ=255)
             tACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCACTAGGCCCGGtt                                                                                                                                         <  2:725428/150‑1 (MQ=255)
             tACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGtt                                                                                                                                         <  1:2051761/150‑1 (MQ=255)
               ccTGCGCATCCTGAAGTGCGGTTTCCGGCGACGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGtttt                                                                                                                                       >  1:62894/1‑150 (MQ=255)
                 tGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGtttttt                                                                                                                                     <  1:639456/150‑1 (MQ=255)
                 tGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAGGCCCGGttttct                                                                                                                                     <  1:2278934/150‑3 (MQ=255)
                    gcATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAGGCCCGGTTTTCTGta                                                                                                                                  <  1:710038/150‑2 (MQ=255)
                         cTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATAc                                                                                                                             <  1:1466007/150‑1 (MQ=255)
                          tGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACa                                                                                                                            >  1:1958057/1‑150 (MQ=255)
                           gAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACAt                                                                                                                           <  2:1012571/150‑1 (MQ=255)
                              gTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTa                                                                                                                        >  2:2113232/1‑150 (MQ=255)
                                   gTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGtt                                                                                                                   <  1:1362708/150‑1 (MQ=255)
                                    tttCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGttt                                                                                                                  <  1:1635591/150‑1 (MQ=255)
                                         gcgcTGGCGACAACGCGCCTATGGCGTATGTCGAGTTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGACTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAtttt                                                                                                             >  1:1600402/1‑149 (MQ=255)
                                          cgcTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCt                                                                                                            <  2:1958057/150‑1 (MQ=255)
                                           gcTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTa                                                                                                           <  2:1864136/150‑1 (MQ=255)
                                            cTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAt                                                                                                          <  2:1336902/150‑1 (MQ=255)
                                                 gACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCACTAGGCCCGGtttt                                                                                                                                       <  2:1515124/116‑1 (MQ=255)
                                                 gACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCACTAGGCCCGGtttt                                                                                                                                       >  1:1515124/1‑116 (MQ=255)
                                                  aCAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGTGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGt                                                                                                    <  1:2102793/150‑1 (MQ=255)
                                                        cgcCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCa                                                                                              <  1:1532654/150‑1 (MQ=255)
                                                          ccTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATg                                                                                            <  2:663058/150‑1 (MQ=255)
                                                            tATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTa                                                                                          >  2:597537/1‑150 (MQ=255)
                                                            tATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAGGCCCGGttttc                                                                                                                                      <  1:2148334/106‑2 (MQ=38)
                                                                gCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCACTAGGCCCGGTTTTCTGtagca                                                                                                                               <  2:2271116/109‑5 (MQ=37)
                                                                 cctATGTCGAGCTCGTTGATCGTCCAGTAGGAGGTGCTGTAGAAGATGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGaa                                                                                     <  2:1600402/148‑1 (MQ=255)
                                                                 cGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGaa                                                                                     <  1:669593/150‑1 (MQ=255)
                                                                  gTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAttt                                                                                                              >  1:732936/1‑124 (MQ=255)
                                                                  gTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAttt                                                                                                              <  2:732936/124‑1 (MQ=255)
                                                                  gTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAAt                                                                                    <  2:116783/150‑1 (MQ=255)
                                                                    aTGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAAttt                                                                                  <  2:314017/150‑1 (MQ=255)
                                                                      gTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGt                                                                                <  1:49721/150‑1 (MQ=255)
                                                                             tGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTtata                                                                                                                              <  1:2409856/97‑1 (MQ=255)
                                                                             tGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTtata                                                                                                                              >  2:2409856/1‑97 (MQ=255)
                                                                               gTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAt                                                                                                          <  2:1387433/115‑1 (MQ=255)
                                                                               gTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAt                                                                                                          >  1:1387433/1‑115 (MQ=255)
                                                                                       tCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCg                                                               >  1:919772/1‑150 (MQ=255)
                                                                                           gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGtt                                                                                                                                  >  1:952255/1‑79 (MQ=255)
                                                                                           gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGtt                                                                                                                                  <  2:952255/79‑1 (MQ=255)
                                                                                           gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTCATCAATTTGTCCCTGTAATCTTGCTCGTCaa                                                           <  1:55249/150‑1 (MQ=255)
                                                                                           gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCaa                                                           >  1:1472609/1‑150 (MQ=255)
                                                                                           gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCaa                                                           >  2:5891/1‑150 (MQ=255)
                                                                                            tCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGCAAGTCATGTAATGAATTTGTCCCTGTAATCTTGc                                                                  >  2:1557114/1‑142 (MQ=255)
                                                                                            tCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGCAAGTCATGTAATGAATTTGTCCCTGTAATCTTGc                                                                  <  1:1557114/142‑1 (MQ=255)
                                                                                              ggcggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTCTCCCTGTAATCTTCCTCCTCactcc                                                        >  2:1260719/1‑146 (MQ=255)
                                                                                                cggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActc                                                      <  2:1472609/150‑1 (MQ=255)
                                                                                                  gTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctc                                                    <  1:1490984/150‑1 (MQ=255)
                                                                                                       gTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCAt                                                                                             <  2:1430648/104‑1 (MQ=255)
                                                                                                       gTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTAc                                               <  2:919772/150‑1 (MQ=255)
                                                                                                        tAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACa                                              >  1:2179115/1‑150 (MQ=255)
                                                                                                          gAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctc                                                    <  2:67429/142‑1 (MQ=255)
                                                                                                          gAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctc                                                    >  1:67429/1‑142 (MQ=255)
                                                                                                          gAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGc                                            <  1:438326/150‑1 (MQ=255)
                                                                                                          gAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGc                                            >  2:224591/1‑150 (MQ=255)
                                                                                                           aaGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGtt                                                                                                                                         <  2:398786/56‑1 (MQ=255)
                                                                                                           aaGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGtt                                                                                                                                         >  1:398786/1‑56 (MQ=255)
                                                                                                           aaGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGcc                                           <  1:2272872/150‑1 (MQ=255)
                                                                                                            aGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCg                                          >  1:880221/1‑150 (MQ=255)
                                                                                                             gctgctGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGt                                                                          >  1:2312668/1‑117 (MQ=255)
                                                                                                             gctgctGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGt                                                                          <  2:2312668/117‑1 (MQ=255)
                                                                                                               tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAt                                                                                                          >  2:1150477/1‑83 (MQ=255)
                                                                                                               tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTAt                                                                                                          <  1:1150477/83‑1 (MQ=255)
                                                                                                               tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTa                                                >  2:1064740/1‑141 (MQ=255)
                                                                                                               tgctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTa                                                <  1:1064740/141‑1 (MQ=255)
                                                                                                                gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTa                                      <  1:983239/150‑1 (MQ=255)
                                                                                                                gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTa                                      <  2:1523866/150‑1 (MQ=255)
                                                                                                                  tgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGc                                    >  1:111751/1‑150 (MQ=255)
                                                                                                                     gTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCgg                                 <  1:597537/150‑1 (MQ=255)
                                                                                                                      tAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGc                                >  2:189445/1‑150 (MQ=255)
                                                                                                                       aaGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGcc                               >  2:2025042/1‑150 (MQ=255)
                                                                                                                             ttCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGaa                         <  1:224591/150‑1 (MQ=255)
                                                                                                                                tgtctgCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGAc                      <  2:2070772/150‑1 (MQ=255)
                                                                                                                                  tctgCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCg                    >  2:585257/1‑150 (MQ=255)
                                                                                                                                    tgCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCAATTGTTTATTGCTATTGATGTAAGTCATGTAATGAATTTTTCCCTGTAATCTTGCTCGGCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGaa                  >  1:133900/1‑150 (MQ=255)
                                                                                                                                    gcagaCAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACa                                              <  1:389073/118‑1 (MQ=255)
                                                                                                                                    gcagaCAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACa                                              >  2:389073/5‑122 (MQ=255)
                                                                                                                                     taaaCAAAAAACCGGGCACTAGGCCCGGTTTTCTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAAttt                                                                                  <  2:1006451/82‑1 (MQ=21)
                                                                                                                                      cTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAAt                <  1:2025042/150‑1 (MQ=255)
                                                                                                                                       tACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg               <  1:2113232/150‑1 (MQ=255)
                                                                                                                                       aacaaaaaaCCGGGCACTAGGCCCGGTTTTCTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg               >  1:1338603/6‑150 (MQ=255)
                                                                                                                                          aaaaaaCCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActct                                                     <  2:617719/107‑1 (MQ=39)
                                                                                                                                          aaaaaaCCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActct                                                     >  1:617719/3‑109 (MQ=39)
                                                                                                                                          aGAAAACCTGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGgag            <  1:1197971/150‑1 (MQ=255)
                                                                                                                                          aGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGgag            >  2:2140808/1‑150 (MQ=255)
                                                                                                                                          aGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGgag            >  2:1548172/1‑150 (MQ=255)
                                                                                                                                           gAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActc                                                      <  1:454103/107‑1 (MQ=255)
                                                                                                                                           gAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActc                                                      >  2:454103/1‑107 (MQ=255)
                                                                                                                                             aaaCCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTCATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGcc                               <  1:53193/128‑1 (MQ=255)
                                                                                                                                             aaaCCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGcc                               >  2:53193/1‑128 (MQ=255)
                                                                                                                                              aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTTTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCATTGAAg                        <  1:902882/134‑1 (MQ=255)
                                                                                                                                              aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCATTGAAg                        >  2:902882/1‑134 (MQ=255)
                                                                                                                                               aCCGGGCGTATTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTTTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGccc       <  2:284980/150‑1 (MQ=255)
                                                                                                                                                ccGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGt                                                              <  2:614529/94‑1 (MQ=255)
                                                                                                                                                ccGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGt                                                              >  1:614529/1‑94 (MQ=255)
                                                                                                                                                    gCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTATCCGGCCCTTTAAGACCTAAATGGAGAGCCCATGAg  >  2:437837/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAG  >  NC_002947/564231‑564528

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: