Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,544,689 A→T 7.0% R52R (CGT→CGA 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,6890AT 7.0% 242.4 / 6.4 71R52R (CGT→CGAbioHpimeloyl‑ACP methyl ester carboxylesterase
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base T (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GCCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGT  >  NC_000913/3544543‑3544837
                                                                                                                                                  |                                                                                                                                                    
gcCTGAACACGCTCGGGATGGGTTAACGCAAACTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACg                                                                                                                                                    >  4:678279/1‑149 (MQ=255)
   tGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGccc                                                                                                                                                 <  1:679095/149‑1 (MQ=255)
      acacGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGaa                                                                                                                                              <  4:667263/149‑1 (MQ=255)
        acGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGc                                                                                                                                            <  1:669298/149‑1 (MQ=255)
        acGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGc                                                                                                                                            <  3:400275/149‑1 (MQ=255)
         cGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGcc                                                                                                                                           <  2:131516/149‑1 (MQ=255)
            tCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCggg                                                                                                                                        <  3:456438/149‑1 (MQ=255)
             cGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgggg                                                                                                                                       <  2:143873/149‑2 (MQ=255)
             cGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGc                                                                                                                                       <  1:306363/149‑1 (MQ=255)
                gATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAgg                                                                                                                                    >  3:46918/1‑149 (MQ=255)
                     gTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaac                                                                                                                               >  1:154520/1‑149 (MQ=255)
                      ttAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaaca                                                                                                                              >  1:20363/1‑149 (MQ=255)
                      ttAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaaca                                                                                                                              >  1:200474/1‑149 (MQ=255)
                       tAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaacaa                                                                                                                             <  3:7241/149‑1 (MQ=255)
                        aaCGCAATCTGGCTTGCTACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAg                                                                                                                            <  2:518274/149‑1 (MQ=255)
                         aCGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAgg                                                                                                                           <  4:46918/149‑1 (MQ=255)
                          cGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGt                                                                                                                          >  4:345057/1‑149 (MQ=255)
                          cGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGt                                                                                                                          >  2:251316/1‑149 (MQ=255)
                           gCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTg                                                                                                                         >  2:523509/1‑149 (MQ=255)
                             aaTCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCa                                                                                                                       >  1:462668/1‑149 (MQ=255)
                               tCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATACCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGc                                                                                                                     >  2:410803/1‑149 (MQ=255)
                                          accaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGa                                                                                                          >  4:248326/1‑149 (MQ=255)
                                          accaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGa                                                                                                          >  3:123924/1‑149 (MQ=255)
                                              gcccgcccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTa                                                                                                      >  4:108483/1‑149 (MQ=255)
                                               cccgcccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTaa                                                                                                     <  4:352051/149‑1 (MQ=255)
                                                ccgcccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAg                                                                                                    >  3:435548/1‑149 (MQ=255)
                                                   cccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCggg                                                                                                                                        <  4:526429/110‑1 (MQ=255)
                                                   cccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCggg                                                                                                                                        >  3:526429/1‑110 (MQ=255)
                                                    ccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTcc                                                                                                >  4:445146/1‑149 (MQ=255)
                                                    ccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTcc                                                                                                >  1:456585/1‑149 (MQ=255)
                                                      aGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTc                                                                                              <  2:424603/149‑1 (MQ=255)
                                                          tCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCa                                                                                          >  1:694865/1‑149 (MQ=255)
                                                           ccAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCaa                                                                                         <  4:123924/149‑1 (MQ=255)
                                                                cTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCaa                                                                                    <  3:504303/149‑1 (MQ=255)
                                                                 tAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAAc                                                                                   >  2:212720/1‑149 (MQ=255)
                                                                 tAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAAc                                                                                   <  2:456585/149‑1 (MQ=255)
                                                                   aCCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTGAAATGCGAGCTAAGTTCCTCGTCAATGCAACGc                                                                                 >  2:484171/1‑149 (MQ=255)
                                                                      aaaTGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCcac                                                                              >  4:745981/1‑149 (MQ=255)
                                                                           gCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTc                                                                         <  4:302652/149‑1 (MQ=255)
                                                                                 tCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTGAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCAtt                                                                   <  1:484171/149‑1 (MQ=255)
                                                                                      tGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTc                                                              <  3:731315/149‑1 (MQ=255)
                                                                                         cTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTcccc                                                           <  3:345057/149‑1 (MQ=255)
                                                                                          tGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCa                                                          >  3:592612/1‑149 (MQ=255)
                                                                                             tGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATcc                                                       <  3:678279/149‑1 (MQ=255)
                                                                                                 gcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgc                                                   >  2:601484/1‑149 (MQ=255)
                                                                                                 gcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgc                                                   >  3:712460/1‑149 (MQ=255)
                                                                                                    cGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcag                                                 <  2:284658/148‑1 (MQ=255)
                                                                                                    cGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcag                                                 >  1:284658/1‑148 (MQ=255)
                                                                                                       cTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagcagc                                             >  2:52414/1‑149 (MQ=255)
                                                                                                       cTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagcagc                                             <  1:204980/149‑1 (MQ=255)
                                                                                                       cTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagcagc                                             >  4:719687/1‑149 (MQ=255)
                                                                                                           ggCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACCCTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACaa                                         >  2:109669/1‑149 (MQ=255)
                                                                                                            gCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAg                                        >  1:50770/1‑149 (MQ=255)
                                                                                                             ccATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTa                                                                                                                  <  1:720476/74‑1 (MQ=255)
                                                                                                             ccATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTa                                                                                                                  >  2:720476/1‑74 (MQ=255)
                                                                                                             ccATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGa                                       >  1:386619/1‑149 (MQ=255)
                                                                                                              cATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGAt                                      <  3:334847/149‑1 (MQ=255)
                                                                                                               atatCCGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCCAGGCACCGCCACACTCCGGCATTCAGTCCCCCTCCGTGCAGCAGCACAAGATg                                     >  2:167404/1‑149 (MQ=255)
                                                                                                               atatCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATg                                     >  3:9032/1‑149 (MQ=255)
                                                                                                                    aGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACAt                                >  2:279859/1‑149 (MQ=255)
                                                                                                                     gCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACAtt                               >  1:770245/1‑149 (MQ=255)
                                                                                                                      cAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTc                              <  1:52414/149‑1 (MQ=255)
                                                                                                                              aGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGAcc                      >  3:689260/1‑149 (MQ=255)
                                                                                                                                   aCCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTgg                 >  1:723500/1‑149 (MQ=255)
                                                                                                                                    ccAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGt                <  1:251316/149‑1 (MQ=255)
                                                                                                                                      aaaTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCt              <  1:246548/149‑1 (MQ=255)
                                                                                                                                        aTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGc            <  3:72244/149‑1 (MQ=255)
                                                                                                                                         tgcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCATTGGTCTGcc           >  2:525996/3‑149 (MQ=255)
                                                                                                                                         tCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcc           >  2:38752/1‑149 (MQ=255)
                                                                                                                                               gCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGa     <  2:279231/149‑1 (MQ=255)
                                                                                                                                               gCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGa     <  3:445146/149‑1 (MQ=255)
                                                                                                                                                  aCGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGt  <  2:196790/149‑1 (MQ=255)
                                                                                                                                                  aCGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGt  <  4:592612/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGT  >  NC_000913/3544543‑3544837

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