Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,073,144 A→T 9.1% intergenic (‑97/+24) plsC ← / ← PP_1845 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,073,1440AT 9.1% 988.2 / 64.2 324intergenic (‑97/+24)plsC/PP_18451‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (144/150);  new base T (12/18);  total (156/168)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 336 total aligned reads displayed.

CCAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCGG  >  NC_002947/2072995‑2073292
                                                                                                                                                      |                                                                                                                                                    
ccAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCaa                                                                                                                                                      <  2:4526499/150‑1 (MQ=255)
ccAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCaa                                                                                                                                                      <  2:718539/150‑1 (MQ=255)
ccAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCaa                                                                                                                                                      <  1:6571241/150‑1 (MQ=255)
 cAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAg                                                                                                                                                     >  1:3695952/1‑150 (MQ=255)
 cAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAg                                                                                                                                                     >  1:848691/1‑150 (MQ=255)
    cGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGccc                                                                                                                                                  <  2:2246898/150‑1 (MQ=255)
          gaaacggCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCTTGCCCGCGaaa                                                                                                                                            <  2:4956919/144‑1 (MQ=255)
          gAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCTTGCCCGCGaaa                                                                                                                                            <  2:1084856/150‑1 (MQ=255)
          gAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaa                                                                                                                                            <  1:6744201/150‑1 (MQ=255)
           aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCTTGCCCGCGaaag                                                                                                                                           <  1:1862632/150‑2 (MQ=255)
           aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                           <  1:4234867/150‑1 (MQ=255)
             aaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAgg                                                                                                                                         <  1:7156754/150‑1 (MQ=255)
                    ccccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTg                                                                                                                                  >  2:6417705/1‑150 (MQ=255)
                    ccccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTg                                                                                                                                  >  2:2768128/1‑150 (MQ=255)
                     cccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGc                                                                                                                                 <  2:5052884/150‑1 (MQ=255)
                     cccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGc                                                                                                                                 <  1:572512/150‑1 (MQ=255)
                      ccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCg                                                                                                                                <  1:2557749/150‑1 (MQ=255)
                      ccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCg                                                                                                                                >  1:4322891/1‑150 (MQ=255)
                         gcagGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCt                                                                                                                             >  2:6173809/1‑150 (MQ=255)
                            ggAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTTGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGc                                                                                                                          <  1:3128101/150‑1 (MQ=255)
                            ggAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGc                                                                                                                          >  2:3987959/1‑150 (MQ=255)
                            ggAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGc                                                                                                                          >  2:5594935/1‑150 (MQ=255)
                               aCATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCaa                                                                                                                       >  2:3348016/1‑150 (MQ=255)
                               aCATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCaa                                                                                                                       <  2:675095/150‑1 (MQ=255)
                                 aTACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAata                                                                                                                     <  2:1673596/150‑1 (MQ=255)
                                   aCGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtg                                                                                                                   >  1:7261866/1‑150 (MQ=255)
                                    cGCAATGAATAGAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgc                                                                                                                  <  2:5231201/150‑1 (MQ=255)
                                      cAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                                <  2:5736997/150‑1 (MQ=255)
                                       aaTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGaaa                                                                                                                                            <  1:456409/122‑1 (MQ=255)
                                       aaTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctgc                                                                                                               <  2:7261866/150‑1 (MQ=255)
                                       aaTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctgc                                                                                                               >  2:3630767/1‑150 (MQ=255)
                                          gAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGa                                                                                                            <  2:5100290/150‑1 (MQ=255)
                                          gAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGa                                                                                                            <  1:2715263/150‑1 (MQ=255)
                                            aTAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATc                                                                                                          <  2:3156636/150‑1 (MQ=255)
                                             tAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATcc                                                                                                         >  1:2240716/1‑150 (MQ=255)
                                              aaaGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCa                                                                                                        <  1:3861834/150‑1 (MQ=255)
                                              aaaGCATGGGACGACTCACACCAGACGGTCGCGGGGGAGTGTACGGGTGTGCACTGGCTGGGTCAATCCTTGGTGCTGGCCTGTTACGGCCCTTTCGCGGGCTAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCa                                                                                                        <  2:93330/150‑1 (MQ=255)
                                               aaGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAg                                                                                                       >  2:4114932/1‑150 (MQ=255)
                                               aaGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAg                                                                                                       <  2:6084932/150‑1 (MQ=255)
                                                aGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGATGATCCAGc                                                                                                      >  1:6863246/1‑150 (MQ=255)
                                                     ggTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGag                                                                                                 >  1:364190/1‑149 (MQ=255)
                                                      gTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                                <  1:7192370/150‑1 (MQ=255)
                                                      gTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                                <  1:7050788/150‑1 (MQ=255)
                                                      gTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAAGCTGCTGATCCAGCTCGATc                                                                                                >  1:3780821/1‑150 (MQ=255)
                                                       tACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCg                                                                                               >  2:6957870/1‑150 (MQ=255)
                                                           aCTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCtt                                                                                           <  2:5892094/150‑1 (MQ=255)
                                                             tCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGt                                                                                         >  1:2029950/1‑150 (MQ=255)
                                                             tCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGt                                                                                         >  2:2334462/1‑150 (MQ=255)
                                                               acacCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTcc                                                                                       <  1:1242851/150‑1 (MQ=255)
                                                                  ccAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAgc                                                                                    <  2:905775/150‑1 (MQ=255)
                                                                          tCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCcag                                                                            <  2:475950/150‑1 (MQ=255)
                                                                          tCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCcag                                                                            >  1:6704601/1‑150 (MQ=255)
                                                                                 gCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCt                                                                     >  2:7421014/1‑150 (MQ=255)
                                                                                   aGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCttt                                                                   <  1:208220/150‑1 (MQ=255)
                                                                                   aGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCttt                                                                   <  1:6175775/150‑1 (MQ=255)
                                                                                       tACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgcgc                                                               >  1:193618/1‑150 (MQ=255)
                                                                                             tgtgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTc                                                         <  2:1999360/150‑1 (MQ=255)
                                                                                              gtgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaa                                                                                                                                            <  2:6024779/66‑1 (MQ=255)
                                                                                              gtgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaa                                                                                                                                            >  1:6024779/1‑66 (MQ=255)
                                                                                              gtgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCa                                                        >  2:2477894/1‑150 (MQ=255)
                                                                                                gCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGc                                                      <  2:5522332/150‑1 (MQ=255)
                                                                                                 cACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCt                                                     >  1:882872/1‑150 (MQ=255)
                                                                                                   ctggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                                   <  2:6335428/150‑1 (MQ=255)
                                                                                                   ctggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                                   <  1:1198444/150‑1 (MQ=255)
                                                                                                         ggggCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTc                                             <  2:2870038/150‑1 (MQ=255)
                                                                                                             cAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTtgt                                         <  1:5623805/150‑1 (MQ=255)
                                                                                                             cAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTtgt                                         <  2:5146403/150‑1 (MQ=255)
                                                                                                             cAAACCCTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTtgt                                         <  2:3695952/150‑1 (MQ=255)
                                                                                                               aaCCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTgc                                       <  2:1235004/150‑1 (MQ=255)
                                                                                                                  cTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGcc                                    >  1:2966912/1‑150 (MQ=255)
                                                                                                                   ttGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCa                                   <  1:4839400/150‑1 (MQ=255)
                                                                                                                     ggTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGc                                 <  1:1375899/150‑1 (MQ=255)
                                                                                                                      gTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCa                                <  1:76309/150‑1 (MQ=255)
                                                                                                                          aGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGtt                            <  2:4433150/150‑1 (MQ=255)
                                                                                                                          aGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGtt                            >  2:2960880/1‑150 (MQ=255)
                                                                                                                          aGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGtt                            <  2:6516184/150‑1 (MQ=255)
                                                                                                                           ggCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTc                           <  2:1589961/150‑1 (MQ=255)
                                                                                                                           ggCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTc                           <  2:5921151/150‑1 (MQ=255)
                                                                                                                             ccTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAg                         >  2:7373005/1‑150 (MQ=255)
                                                                                                                              cTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGc                        <  1:6417705/150‑1 (MQ=255)
                                                                                                                               tGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGct                       >  2:4312120/1‑150 (MQ=255)
                                                                                                                                gTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGctt                      >  1:1690181/1‑150 (MQ=255)
                                                                                                                                 ttACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGcttc                     <  2:7106963/150‑1 (MQ=255)
                                                                                                                                 ttACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCACCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCAAGTTCAGctta                     >  1:5209257/1‑149 (MQ=255)
                                                                                                                                    cGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCcagca                                                                          >  2:5955713/1‑94 (MQ=255)
                                                                                                                                    cGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCttcttc                  >  2:3001874/1‑150 (MQ=255)
                                                                                                                                     ggCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa                 <  2:848691/150‑1 (MQ=255)
                                                                                                                                     ggCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa                 <  2:4847612/150‑1 (MQ=255)
                                                                                                                                        ccttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCt              >  2:1834832/4‑150 (MQ=255)
                                                                                                                                        ccttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCt              >  1:3532082/4‑150 (MQ=255)
                                                                                                                                         cttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc             >  1:587696/3‑150 (MQ=255)
                                                                                                                                         cTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc             >  1:4635516/1‑150 (MQ=255)
                                                                                                                                         cTTTCGCGGGCAAGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTcc            >  1:3730919/1‑150 (MQ=255)
                                                                                                                                          ttttCGCGGGCTTGCCCGCGAAAGGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc             >  1:2869102/2‑150 (MQ=255)
                                                                                                                                          ttttCGCGGGCTTGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTc             >  2:1053816/2‑150 (MQ=255)
                                                                                                                                          ttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAg                                                                                     <  2:5416989/76‑1 (MQ=37)
                                                                                                                                           tttCGCGGGCAAGGCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTcc            >  2:3787457/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCGTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTccc           >  1:501485/1‑150 (MQ=255)
                                                                                                                                                    cAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCg   <  2:2638340/150‑1 (MQ=255)
                                                                                                                                                     aaGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCgg  >  1:1678500/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCGG  >  NC_002947/2072995‑2073292

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