Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,544,693 C→G 6.6% G51A (GGG→GCG) ‡ bioH ← pimeloyl‑ACP methyl ester carboxylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,544,6930CG 6.6% 351.8 / 12.5 118G51A (GGG→GCG) ‡bioHpimeloyl‑ACP methyl ester carboxylesterase
Reads supporting (aligned to +/- strand):  ref base C (66/44);  new base G (3/5);  total (69/49)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

Only 100 of 121 total aligned reads displayed.

CCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATT  >  NC_000913/3544544‑3544841
                                                                                                                                                     |                                                                                                                                                    
    aaCACGCGCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGaa                                                                                                                                                  <  2:362583/150‑1 (MQ=255)
      cacGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGc                                                                                                                                                >  2:559669/1‑150 (MQ=255)
       acGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGcc                                                                                                                                               >  1:1273451/1‑150 (MQ=255)
        cGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCg                                                                                                                                              >  2:657005/1‑150 (MQ=255)
        cGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCg                                                                                                                                              >  2:637289/1‑150 (MQ=255)
         gCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCGGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCgg                                                                                                                                             >  1:1168173/1‑150 (MQ=255)
          cTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGGCTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCggg                                                                                                                                            >  1:557497/1‑150 (MQ=255)
          cTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCggg                                                                                                                                            >  1:32731/1‑150 (MQ=255)
          cTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCggg                                                                                                                                            >  1:1669106/1‑150 (MQ=255)
           tCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgggg                                                                                                                                           <  2:579201/150‑2 (MQ=255)
           tCGGGATGGGTTAACGCAATCTGGCTTGCCACCACCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAAACCCCGGCGTCGGGCGTAGCCgggg                                                                                                                                           <  2:1461837/150‑2 (MQ=255)
            ctggATGGTTTAACGCAAGCTGGCTTGCCACCAGCCCGCCCAGACTCCAGACTAACCAAATTGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgggga                                                                                                                                          <  2:1052628/148‑3 (MQ=255)
            cGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgggga                                                                                                                                          <  2:1685231/150‑3 (MQ=255)
             gggATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAAAGGCTTTATCAGGTGCCTGTTGCAGCCCGGCTACGGCCATATCAGCAAGTGACAGCGCACCAAAACCCCGGCAACGCCCGACGCCGGGccc                                                                                                                                         >  1:603399/1‑148 (MQ=255)
              ggATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGCCTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGCCAGCGCACCACATCCCCGGCTACGCCCGAAGCCGGGCAgg                                                                                                                                        >  1:1224809/1‑150 (MQ=255)
              ggATGGGTTAACGCAATCTGGCTTGCCACAAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAgg                                                                                                                                        >  1:153440/1‑150 (MQ=255)
                  gggTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaa                                                                                                                                    <  1:510675/150‑1 (MQ=255)
                      tAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAg                                                                                                                                >  1:458128/1‑150 (MQ=255)
                       aaCGCAATCTGGCTTGCCACCAGCCAGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAgg                                                                                                                               <  1:999839/150‑1 (MQ=255)
                               cTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGt                                                                                                                       >  2:1710847/1‑150 (MQ=255)
                               cTGGCTTGCCACAAGCCCGCCCAGACTCCCGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCACATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCgg                                                                                                                       >  1:1111663/1‑149 (MQ=255)
                                       ccaccaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGCCAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAGGGTGCAGCGTAAAATGCg                                                                                                               >  1:167365/1‑150 (MQ=255)
                                       ccaccaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGCTCAACAAGGTGCAGCGTAAAATGCg                                                                                                               >  2:154888/1‑150 (MQ=255)
                                       ccaccaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGAGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCg                                                                                                               >  1:565491/1‑150 (MQ=255)
                                        caccaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCTGGTCAACAAGGTGCAGCGTAAAATGCGa                                                                                                              >  2:66722/1‑150 (MQ=255)
                                        caccaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGa                                                                                                              >  1:414229/1‑150 (MQ=255)
                                         accaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAg                                                                                                             <  1:26400/150‑1 (MQ=255)
                                           caGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCt                                                                                                           >  2:1395977/1‑150 (MQ=255)
                                           caGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCt                                                                                                           <  2:818347/150‑1 (MQ=255)
                                            aGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAACTGCGAGCTa                                                                                                          >  1:717314/1‑150 (MQ=255)
                                            aGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTa                                                                                                          >  2:316607/1‑150 (MQ=255)
                                             gcccgcccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGACATATCAGAAAGTGCCAGCGCCCCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGTGCTaa                                                                                                         >  2:388037/1‑150 (MQ=255)
                                              cccgcccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAg                                                                                                        >  1:657924/1‑150 (MQ=255)
                                                  cccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGGTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTcc                                                                                                    <  1:1630075/150‑1 (MQ=255)
                                                   ccAGCCACCAGCCTAACCAAATGGCTTTATCTGGGGCCTGTTGCAGCACGGCTTCGGCCATATCAGCATGCGACAGCGCACCAAATCCGCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGGGCAGCGTAAAATGCGAGCTAAGTTCCt                                                                                                   <  2:430709/150‑1 (MQ=255)
                                                     aGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCg                                                                                                 >  1:1448386/1‑150 (MQ=255)
                                                     aGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCg                                                                                                 <  2:1465861/150‑1 (MQ=255)
                                                     aGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCg                                                                                                 >  2:586191/1‑150 (MQ=255)
                                                       aCTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTc                                                                                               <  1:1550924/150‑1 (MQ=255)
                                                        cTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCa                                                                                              >  2:1274560/1‑150 (MQ=255)
                                                        cTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCa                                                                                              >  1:389843/1‑150 (MQ=255)
                                                           cAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATg                                                                                           <  1:975358/150‑1 (MQ=255)
                                                            aGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGc                                                                                          <  1:1466819/150‑1 (MQ=255)
                                                               cTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAAc                                                                                       <  1:188713/150‑1 (MQ=255)
                                                               cTAACCAAATGGCTTTATCAGCTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAAc                                                                                       <  2:1722782/150‑1 (MQ=255)
                                                                tacacAAATGGCTTTATCAGGAGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTAAATGCAACg                                                                                      <  1:1308626/146‑1 (MQ=255)
                                                                     aaaTGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGTTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCcaca                                                                                 <  1:177389/150‑1 (MQ=255)
                                                                           cTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCgg                                                                           <  2:414229/150‑1 (MQ=255)
                                                                           cTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCgg                                                                           <  1:214042/150‑1 (MQ=255)
                                                                           cTTTAGCAGGTGCCGGGTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGGCAATGCAACGCCACACTTCgg                                                                           <  2:277845/150‑1 (MQ=255)
                                                                            tttATCATGTGCCTGTTGCACGACGGCTTCGGCCATATCATCAACTGTCAGCCCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGc                                                                          <  2:225156/150‑1 (MQ=255)
                                                                            tttATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCCGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGc                                                                          >  1:1311474/1‑150 (MQ=255)
                                                                                tCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGAGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTc                                                                      <  2:565491/150‑1 (MQ=255)
                                                                                tCAGGTGCCTGGTGCATCACGGCTTCGGCCATATCAGCATGTGAGAGCGGACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGTTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCTCACTTCGGCATTc                                                                      <  1:1420842/150‑1 (MQ=255)
                                                                                 cAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCa                                                                     <  1:1024250/150‑1 (MQ=255)
                                                                                     tGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTcc                                                                 >  1:1221503/1‑150 (MQ=255)
                                                                                         tGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACACCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCAt                                                             >  1:1109136/1‑150 (MQ=255)
                                                                                          gTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATc                                                            >  2:1184170/1‑150 (MQ=255)
                                                                                            tGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCg                                                          >  2:1547948/1‑150 (MQ=255)
                                                                                              cagcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCACCAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTg                                                        >  1:1651001/1‑150 (MQ=255)
                                                                                                gcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgca                                                      <  1:316607/150‑1 (MQ=255)
                                                                                                gcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgca                                                      >  1:206073/1‑150 (MQ=255)
                                                                                                  aCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCACTGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagc                                                    >  1:1554038/1‑150 (MQ=255)
                                                                                                  aCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagc                                                    <  1:1184170/150‑1 (MQ=255)
                                                                                                    ggCGTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGGCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAAAGCCACACTTCGGCTTTCAGTCCCCATCCGTgcagcag                                                  <  2:1232102/150‑1 (MQ=255)
                                                                                                     ccTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGAGCAGCGTAAAATGCGAGCAAAGTTACACGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagcagc                                                 >  2:1125590/2‑150 (MQ=255)
                                                                                                      cTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagcagca                                                <  1:1464078/150‑1 (MQ=255)
                                                                                                       ttCGGCCATATCAGCAAGTGACTGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGcac                                               <  2:1091137/150‑1 (MQ=255)
                                                                                                          ggCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCCCACTTCGGCATTCAGTCCCCATCCGTGCAGCAACACAAg                                            >  1:1142164/1‑150 (MQ=255)
                                                                                                           gCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGa                                           >  2:535084/1‑150 (MQ=255)
                                                                                                             cATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATg                                         >  1:1362561/1‑150 (MQ=255)
                                                                                                              atatCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGa                                        <  2:1587212/150‑1 (MQ=255)
                                                                                                                 tCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACa                                     <  1:1499893/150‑1 (MQ=255)
                                                                                                                    gCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTc                                  <  1:1678269/150‑1 (MQ=255)
                                                                                                                        gTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCt                              >  1:892180/1‑150 (MQ=255)
                                                                                                                         tGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCATCGCCACACTTAGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTACACTg                             >  2:300797/1‑150 (MQ=255)
                                                                                                                            cAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGAcc                          >  2:1406280/1‑150 (MQ=255)
                                                                                                                              gcgcACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCtt                        >  1:424362/1‑150 (MQ=255)
                                                                                                                               cgcACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCttt                       >  1:766652/1‑150 (MQ=255)
                                                                                                                                gcACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTg                      >  2:1137446/1‑150 (MQ=255)
                                                                                                                                gcACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTg                      <  1:154888/150‑1 (MQ=255)
                                                                                                                                gcACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCATGACCTTTg                      >  1:927240/1‑150 (MQ=255)
                                                                                                                                   ccAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTc                   >  2:1604904/1‑150 (MQ=255)
                                                                                                                                     aaaTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTg                 <  2:1580078/150‑1 (MQ=255)
                                                                                                                                      cctgcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGc                >  2:322507/5‑150 (MQ=255)
                                                                                                                                       ctgcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcc               >  1:374440/4‑150 (MQ=255)
                                                                                                                                       aTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcc               <  2:657309/150‑1 (MQ=255)
                                                                                                                                       aTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcc               <  1:649875/150‑1 (MQ=255)
                                                                                                                                       aTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcc               <  1:525016/150‑1 (MQ=255)
                                                                                                                                        tgcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcgc              >  1:926035/3‑148 (MQ=255)
                                                                                                                                          cccGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGccacc            >  1:1099260/1‑150 (MQ=255)
                                                                                                                                          cccGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGccacc            >  1:682546/1‑150 (MQ=255)
                                                                                                                                          cccGGCTACGCCCGAAGCAGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGccacc            <  2:1074775/150‑1 (MQ=255)
                                                                                                                                           ccGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGccacca           >  2:592126/1‑150 (MQ=255)
                                                                                                                                             ggCTACGCGCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCACCATCCGTTCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGa         <  2:1082405/150‑1 (MQ=255)
                                                                                                                                                tACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGGCTGCCACCAGATGt      <  1:724784/150‑1 (MQ=255)
                                                                                                                                                 aCGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGtt     <  2:1433285/150‑1 (MQ=255)
                                                                                                                                                   gCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTAt   <  1:622874/150‑1 (MQ=255)
                                                                                                                                                    cccGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTAtt  >  1:1509494/1‑150 (MQ=255)
                                                                                                                                                    cccGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTAtt  <  1:1309624/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATT  >  NC_000913/3544544‑3544841

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