Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,198,945 A→T 8.7% intergenic (‑246/+31) ethA ← / ← PP_2806 FAD‑containing monooxygenase EthA/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,198,9450AT 8.7% 170.9 / 6.1 59intergenic (‑246/+31)ethA/PP_2806FAD‑containing monooxygenase EthA/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (23/30);  new base T (3/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

ACCACGTCACTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGATCA  >  NC_002947/3198797‑3199094
                                                                                                                                                    |                                                                                                                                                     
accacGTCACTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAg                                                                                                                                                      >  2:290064/1‑150 (MQ=255)
 ccacGTCACTGCTATTTCAGGATGATTTGCACCTGTTACCTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCCGCTTTCAGAGCCGGCGGGGGAGGGGATCGCCGTGCAGc                                                                                                                                                     <  1:1204682/150‑1 (MQ=255)
        aCTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACgcag                                                                                                                                              <  1:150453/150‑1 (MQ=255)
          tGCTATTTCAGGATGATTTGCACCTGTTACCTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACgcagca                                                                                                                                            >  2:861257/1‑150 (MQ=255)
                     gatgatTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGATAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATGGCCGGGCACCGCCGCAGCACGGCGAATccc                                                                                                                                 >  2:595316/1‑149 (MQ=255)
                      atgatTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCtg                                                                                                                                <  2:1337447/150‑1 (MQ=255)
                        gatTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCtgtg                                                                                                                              <  1:290064/150‑1 (MQ=255)
                        gatTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCtgtg                                                                                                                              <  2:1338670/150‑1 (MQ=255)
                          tttGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTggg                                                                                                                            >  1:492856/1‑150 (MQ=255)
                          tttGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGAACACCATGCGCTGACGCAGCACGGCGAATCCTGTggg                                                                                                                            >  2:889912/1‑150 (MQ=255)
                                       ctcttccGATC‑ACGATCAGGCTGGGTAAGTCTGCTGGGTATAGTACTCAAAGCGATGTAAGTGGGGCGCGGCGTCTTTCAGAGCCGGAGCGGGAGGGGATCGCCGTGCAGCGCCGCAGCACGGCGAATCCTGTTGGGCCTCAGTAACTgg                                                                                                              <  2:667004/145‑1 (MQ=255)
                                       cTTCGCCGATAAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTg                                                                                                               <  1:1273001/150‑1 (MQ=255)
                                           gCCGATCAACGATCAGACTGGGTAAGTCTTGTAGCGAGCGTCATCGGAACGTTCGAAGTGGGGCGCGGCGGCTTTCAGAGCCGGAGGGGGAGGGGAACGGCTTGCAGAGACGCAGCACGGGGAATCCTGGGGGGCCTGAGTAACTGGcca                                                                                                           >  1:691189/1‑150 (MQ=255)
                                            ccGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGccac                                                                                                          >  2:762051/1‑150 (MQ=255)
                                               aTCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAg                                                                                                       >  2:436136/1‑150 (MQ=255)
                                                       tCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATACTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGa                                                                                               >  2:499650/1‑150 (MQ=255)
                                                           aCTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCGGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGc                                                                                           <  2:253443/150‑1 (MQ=255)
                                                                 tCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCa                                                                                     >  1:712374/1‑150 (MQ=255)
                                                                     tCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcg                                                                                 >  2:1238414/1‑150 (MQ=255)
                                                                         gTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtct                                                                             >  2:1001266/1‑150 (MQ=255)
                                                                          tAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtctc                                                                            <  1:436136/150‑1 (MQ=255)
                                                                           aGATATCGTCATCGGAACGTGCGAAGGGGGGCGCTGCGGCTTTGAGAGCCGGAGGGGGAGGGGTTGGGCGTGCAGCGACGCAGCAGGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtctcg                                                                           <  2:132415/150‑1 (MQ=255)
                                                                            gATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGATTCGCCGTGCTGCGTCGCTGCACgg                                                                                                                                         <  1:1212427/87‑1 (MQ=37)
                                                                            gATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGATTCGCCGTGCTGCGTCGCTGCACgg                                                                                                                                         >  2:1212427/1‑87 (MQ=37)
                                                                              tatCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGATTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtctcggtc                                                                        <  2:107164/150‑1 (MQ=255)
                                                                               atCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtctcggtct                                                                       >  1:947535/1‑150 (MQ=255)
                                                                                      tCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTc                                                                >  1:592311/1‑150 (MQ=255)
                                                                                         gAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAg                                                             >  2:909384/1‑150 (MQ=255)
                                                                                               tCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGc                                                       <  1:762051/150‑1 (MQ=255)
                                                                                                  taGGGGGGCTCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACg                                                    >  2:1125192/2‑150 (MQ=255)
                                                                                                     gggggTGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGTACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTTTCGGCCTCCAGGCGTGCACGGTc                                                 <  2:400246/150‑1 (MQ=255)
                                                                                                       gggCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCtgtg                                                                                                                              <  1:692272/71‑1 (MQ=255)
                                                                                                       gggCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCtgtg                                                                                                                              >  2:692272/1‑71 (MQ=255)
                                                                                                            cTGCGGCTTTCAGCGCGGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGa                                          <  2:853847/150‑1 (MQ=255)
                                                                                                            cTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGa                                          <  2:1257915/150‑1 (MQ=255)
                                                                                                               cGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCg                                       <  1:1001266/150‑1 (MQ=255)
                                                                                                                   tttCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTg                                                        >  1:830441/1‑129 (MQ=255)
                                                                                                                   tttCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTg                                                        <  2:830441/129‑1 (MQ=255)
                                                                                                                    ttCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCa                                  <  2:492856/150‑1 (MQ=255)
                                                                                                                       agagCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGc                                                                                                                          >  1:345205/1‑59 (MQ=255)
                                                                                                                       agagCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGc                                                                                                                          <  2:345205/59‑1 (MQ=255)
                                                                                                                       agagCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGc                               <  1:1125192/150‑1 (MQ=255)
                                                                                                                       agagCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGc                               >  1:400618/1‑150 (MQ=255)
                                                                                                                         agCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtctcg                                                                           <  2:940631/104‑1 (MQ=255)
                                                                                                                         agCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGcggtctcg                                                                           >  1:940631/1‑104 (MQ=255)
                                                                                                                                gggggaggggATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAg                      <  2:712374/150‑1 (MQ=255)
                                                                                                                                    gaggggaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAg                      >  2:33038/1‑146 (MQ=255)
                                                                                                                                    gaggggaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAg                      <  1:33038/146‑1 (MQ=255)
                                                                                                                                     aggggaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGt                 <  1:861257/150‑1 (MQ=255)
                                                                                                                                      ggggaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTACCTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTc                >  1:564612/1‑150 (MQ=255)
                                                                                                                                         gaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGaa             <  2:947535/150‑1 (MQ=255)
                                                                                                                                         gaTCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACGGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGaa             <  2:625278/150‑1 (MQ=255)
                                                                                                                                           tCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACa           <  1:1238414/150‑1 (MQ=255)
                                                                                                                                             tccGTGCAGCGACGCAGCACGGCGAATCCTCTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGa         <  2:564612/149‑1 (MQ=255)
                                                                                                                                              ccGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGAt        >  1:321429/1‑150 (MQ=255)
                                                                                                                                               cGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTc                                                 <  1:1203919/108‑1 (MQ=255)
                                                                                                                                               cGTGCTGCGTCGCTGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTc                                                 >  2:1203919/1‑108 (MQ=255)
                                                                                                                                               cGTGCTGCGTCGCTGCACGGCG‑ATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTg      >  1:827494/1‑150 (MQ=255)
                                                                                                                                                   cAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGATc   <  1:909384/150‑1 (MQ=255)
                                                                                                                                                    aGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGATCa  >  2:135820/1‑150 (MQ=255)
                                                                                                                                                    aGCGACGCAGCACGGCGAATACTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGATCa  <  1:499650/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCACGTCACTGCTATTTCAGGATGATTTGCACCTGTTACTTCGCCGATCAACGATCAGACTGGGTCAGTCTGGTAGATATCGTCATCGGAACGTTCGAAGGGGGGCGCTGCGGCTTTCAGAGCCGGAGGGGGAGGGGATCGCCGTGCAGCGACGCAGCACGGCGAATCCTGTGGGGCCTCAGTAACTGGCCACCAGCAACTTGACTGCCTGGCAGGCGGTCTCGGTCTCGGCCTCCAGGCGTGCACGGTCGCCGTGAGCGGACCAGGCGAGCATCAGGCCGTCGAACAGATTGATCA  >  NC_002947/3198797‑3199094

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