Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,544,689 A→T 5.4% 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 5.4% 262.8 / 2.5 96R52R (CGT→CGAbioHpimeloyl‑ACP methyl ester carboxylesterase
Reads supporting (aligned to +/- strand):  ref base A (55/34);  new base T (3/3);  total (59/37)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CGCCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGAT  >  NC_000913/3544542‑3544835
                                                                                                                                                   |                                                                                                                                                  
cgcCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTAc                                                                                                                                                   <  2:1217764/149‑1 (MQ=255)
  ccTGAACCCGCTCGGGATGGGTAAACGCAATCTGGCGTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTGCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGc                                                                                                                                                 <  2:82428/149‑1 (MQ=255)
  ccTGAACACGCTCGGGATGGGTTAACGGAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCGGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGc                                                                                                                                                 <  2:263258/149‑1 (MQ=255)
  ccTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGc                                                                                                                                                 >  2:661218/1‑149 (MQ=255)
  ccTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGc                                                                                                                                                 <  2:607346/149‑1 (MQ=255)
  ccTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGc                                                                                                                                                 <  2:1430819/149‑1 (MQ=255)
   cTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGcc                                                                                                                                                >  2:51533/1‑149 (MQ=255)
      aaCACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCgc                                                                                                                                             >  1:76025/1‑148 (MQ=255)
           gCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAAACAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGcag                                                                                                                                        >  2:735222/1‑147 (MQ=255)
            cTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgg                                                                                                                                       <  2:1531441/149‑1 (MQ=255)
             tagggATGGGTTAAGCCAATCTGGCGTGCAACCAGCCCGCCCAGCCTCCTTCCTAACCAATTGTCTTTATCAGGTCCCTGTAGCAGCACGTCTTCGGCCATATCAGCAAGTGACAGAGCCCCAAATCCCCGGCTTCGGGCGTAGCCggg                                                                                                                                      <  2:843286/147‑1 (MQ=255)
               gggATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCa                                                                                                                                    >  2:1354696/1‑149 (MQ=255)
                 gATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCTAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAgg                                                                                                                                  >  1:405377/1‑149 (MQ=255)
                 gATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAgg                                                                                                                                  >  1:45886/1‑149 (MQ=255)
                  aTGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCCGGTGCCTGGTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGt                                                                                                                                 <  1:1262428/149‑1 (MQ=255)
                    gggTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTca                                                                                                                               <  1:5212/149‑1 (MQ=255)
                       ttAACGCAATCTGGCTTGCCACCAGCCCGCGCAGACTCCAGCCTAAGCAAATGGCTTAATCAGGTGCCTGTTGCAGCACGTCTTCGGACATATCAGCAAGAGACAGAGCACAAAATCCCAGTCTACCCCCGAAGCCGGGAAGGacaaca                                                                                                                            >  2:857961/1‑149 (MQ=255)
                       ttAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGCAGCCGGGCAGGTcaaca                                                                                                                            >  2:463985/1‑149 (MQ=255)
                       ttAACGCAAGCTGGCGTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaaca                                                                                                                            <  2:1196107/149‑1 (MQ=255)
                         aaCGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAg                                                                                                                          >  2:441792/1‑149 (MQ=255)
                         aaCGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCGGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAg                                                                                                                          <  1:243326/149‑1 (MQ=255)
                           cGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTACCCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGCCCGAAGCCGGGCAGGTCAACAAGGt                                                                                                                        >  1:1383627/1‑149 (MQ=255)
                           cGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGt                                                                                                                        >  2:24718/1‑149 (MQ=255)
                           cGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGt                                                                                                                        <  2:979121/149‑1 (MQ=255)
                            gCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTAACAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTg                                                                                                                       >  1:558004/1‑149 (MQ=255)
                               aTCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGGCGCCGGGCCGGTCAACCAGGTGCAg                                                                                                                    >  2:1509814/1‑149 (MQ=255)
                                tCTGGCTTGCCACCAGCCCGCCCAGACTCCCGCCTACCCAATTGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCTGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGc                                                                                                                   <  1:1032026/149‑1 (MQ=255)
                                   ggCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTa                                                                                                                <  2:788160/149‑1 (MQ=255)
                                      ttGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGCAGCCGGGCAGGTCAACAAGGTGCAGCGTaaaa                                                                                                             >  2:405707/1‑149 (MQ=255)
                                                    cccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTc                                                                                               >  2:1392077/1‑149 (MQ=255)
                                                     ccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTcc                                                                                              >  1:248179/1‑149 (MQ=255)
                                                       aGAATCCATCATAACAAATTGGCGTTATCAGGTGCCTGTTGCAGCACGGATTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCGCGGCTACGCCAGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTc                                                                                            <  2:1075367/149‑1 (MQ=255)
                                                          cTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTc                                                                                         <  2:1358207/149‑1 (MQ=255)
                                                          cTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTc                                                                                         <  2:662318/149‑1 (MQ=255)
                                                            ccAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGAGTAAAATGCGAGGTAAGTTCATCGTCaa                                                                                       >  2:377372/1‑149 (MQ=255)
                                                               gCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGc                                                                                    >  2:1309471/1‑149 (MQ=255)
                                                               gCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGc                                                                                    >  2:1191528/1‑149 (MQ=255)
                                                                ccTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAAATCGAGCTAAGTTCCTCGTCAATGCa                                                                                   >  2:668893/1‑149 (MQ=255)
                                                                   aaCCAAATGGCTTTATCAGGTGCCTGTTGCAGCACTGCTTCGGACATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCAGAAGCAGGGCAGCTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTACTCGACAATGCAAcc                                                                                >  2:1481326/1‑148 (MQ=255)
                                                                    aCCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGc                                                                               >  1:227729/1‑149 (MQ=255)
                                                                    aCCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGc                                                                               >  1:351692/1‑149 (MQ=255)
                                                                      cAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGccc                                                                             >  2:663156/1‑148 (MQ=255)
                                                                      cAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCACCGGCTCCGCCCGAAGCCGGGCAGGTCAACAAGGTCCAGCGTAAAATGCGAGCGACGTTCCTCGTCAATGCAACgacc                                                                             >  1:295752/1‑146 (MQ=255)
                                                                            gCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTc                                                                       >  2:1495624/1‑149 (MQ=255)
                                                                             gtttATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCg                                                                      <  1:67358/148‑1 (MQ=255)
                                                                                tATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCa                                                                   >  1:957258/1‑149 (MQ=255)
                                                                                tATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCa                                                                   >  1:213459/1‑149 (MQ=255)
                                                                                 aTCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACCCTTCGGCAt                                                                  >  2:1358650/1‑149 (MQ=255)
                                                                                 aTCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCAt                                                                  <  1:800172/149‑1 (MQ=255)
                                                                                  tCAGGTGCCTGTTGCAGCACGGCTTCGGCAATATCAGCAAGTGACAGCGCACCAAATCCGCGGCTACGCCCGAAGCCGGGCAGGGCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCAtt                                                                 <  2:330154/149‑1 (MQ=255)
                                                                                    aGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACGTCGGCATTcc                                                               >  2:1072469/1‑148 (MQ=255)
                                                                                       tGCCTGGTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTc                                                            <  2:966222/149‑1 (MQ=255)
                                                                                         ccTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTccc                                                          >  1:622439/1‑149 (MQ=255)
                                                                                           tGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCa                                                        >  2:515873/1‑149 (MQ=255)
                                                                                             ttGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATc                                                      >  2:1128451/1‑149 (MQ=255)
                                                                                             ttGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAAGGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATc                                                      >  1:97342/1‑149 (MQ=255)
                                                                                              tGCAGCACGGCTGCGGCCATATCAGCAATTGACTGCTCACCTAATCCCGGGCTTCGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGGAACGCCTCACTTCGGCATTCTGTCCCCATcc                                                     <  2:1007031/149‑1 (MQ=255)
                                                                                                cagcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGt                                                   <  2:838158/149‑1 (MQ=255)
                                                                                                cagcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGt                                                   <  1:1392077/149‑1 (MQ=255)
                                                                                                 agcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTg                                                  >  1:1193491/1‑149 (MQ=255)
                                                                                                  gcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgc                                                 <  1:771450/149‑1 (MQ=255)
                                                                                                  gcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgc                                                 <  1:174499/149‑1 (MQ=255)
                                                                                                     cGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagc                                              >  1:712520/1‑149 (MQ=255)
                                                                                                       gCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTAAGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGAGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGTCATACAGTCCCCATCCGTgcagcag                                            >  2:450883/1‑149 (MQ=255)
                                                                                                        cTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagcagc                                           <  1:1244273/149‑1 (MQ=255)
                                                                                                            ggCCATATCAGCAAGTGACAGCGAAGCAAATCCCCGTCTACGCCCGAAGCCGGGCAGGTATACAAGTTGCAGCGTAAAATGAGAGCTAAGTTCCTCGTCAATGCAACGCCACTCTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACaa                                       <  2:728567/149‑1 (MQ=255)
                                                                                                             gCCATATCAGCAAGTGACAGCGCTCCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTTAGTTCCTCGTCAATGCACCGCCACACTTCGGCATTCAGTCCCCCTCCGTGCAGCAGCACAAg                                      <  1:857961/149‑1 (MQ=255)
                                                                                                             gCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTGCTCGTCAATGCAACGCCACCCTTCGGCAGTCAGTCCCCATCCGTGCAGCAGCACAAg                                      >  2:674434/1‑149 (MQ=255)
                                                                                                                atatCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATg                                   >  1:1192319/1‑149 (MQ=255)
                                                                                                                atatCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGATCCTCGTCAATGCAACGCCACACTGCGGCATTCAGTACCCATCCGTGCAGCAGCACAAGATg                                   >  1:907867/1‑149 (MQ=255)
                                                                                                                  atCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCGGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGaa                                 >  2:946512/1‑149 (MQ=255)
                                                                                                                  atCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGaa                                 <  1:728605/149‑1 (MQ=255)
                                                                                                                  atCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGAAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGaa                                 >  2:876151/1‑149 (MQ=255)
                                                                                                                   tCAGCAAGTGACAGCGCACGAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAAc                                <  1:1174559/149‑1 (MQ=255)
                                                                                                                   tCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAAc                                <  1:661218/149‑1 (MQ=255)
                                                                                                                     aGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCCCACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACAt                              >  2:977500/1‑149 (MQ=255)
                                                                                                                         aGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCCGCAGCACAAGATGAACATTccc                          >  2:1268868/1‑149 (MQ=255)
                                                                                                                          gTGACAGCGCCCCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTcccc                         >  1:434062/1‑149 (MQ=255)
                                                                                                                             aCAGCGCACCAAATCCCCGGCTACGCCCGACGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTgc                      >  2:1421921/1‑148 (MQ=255)
                                                                                                                             aCAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACTCCACCCTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACATGATGAACATTCCCCtta                      >  2:312002/1‑147 (MQ=255)
                                                                                                                              cAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGAc                     >  2:574491/1‑149 (MQ=255)
                                                                                                                              cAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGAc                     <  1:1191528/149‑1 (MQ=255)
                                                                                                                               aGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACCCTTCGGCATTCAGTCCCCATCCGTGCAGCAGCCCCAGATGAACCTTCCCCTGccc                    >  1:1140494/1‑149 (MQ=255)
                                                                                                                                   cACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTg                <  2:1201372/149‑1 (MQ=255)
                                                                                                                                    aCCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTgg               >  1:1195054/1‑149 (MQ=255)
                                                                                                                                    aCCAAATCCCCGGCTACGCCCGAAGACGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCGGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTgg               <  1:188652/149‑1 (MQ=255)
                                                                                                                                      cAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTc             >  1:1555690/1‑149 (MQ=255)
                                                                                                                                      cAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCGCGTCAAGGCAACGCCACCCTTCGGCATTCAGTCCCCGTCCGTGCAGCAGCACGAGATGATCATTCCCCTGCCGTTTgggg             >  1:673914/1‑147 (MQ=255)
                                                                                                                                       acctgcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCTACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCt            >  1:1097362/6‑149 (MQ=255)
                                                                                                                                         aTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCCGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCGCACTTAGGCAGTCAGTCCCCCTCCGTGCAGCAGCCCAAGATGCACCTTCCCCTGACCTTTAGTCAGc          >  2:842144/1‑149 (MQ=255)
                                                                                                                                         aTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACCCTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGc          >  2:1264934/1‑149 (MQ=255)
                                                                                                                                          tCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcc         <  2:622439/149‑1 (MQ=255)
                                                                                                                                           gcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcca        >  2:497310/2‑149 (MQ=255)
                                                                                                                                           ccccGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACCCTTCGGCATTCAGTCCCCCTCCGTGCAGCAGCACAAGATGAACATTCCCCTGCCCTTTGGTCTGccc        >  1:1026323/1‑148 (MQ=255)
                                                                                                                                             ccGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGccacc      >  1:586269/1‑149 (MQ=255)
                                                                                                                                                 cTACGCCCGAAGCCGGGCAGGTCATCAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACTCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTTGTCTGCCACCAGAt  <  1:377372/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGAT  >  NC_000913/3544542‑3544835

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: