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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,544,6930CG 6.5% 254.8 / 7.4 77G51A (GGG→GCG) ‡bioHpimeloyl‑acyl carrier protein methyl ester esterase
Reads supporting (aligned to +/- strand):  ref base C (35/36);  new base G (2/3);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

TGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTC  >  NC_000913/3544546‑3544842
                                                                                                                                                   |                                                                                                                                                     
tGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCg                                                                                                                                                     >  2:56292/1‑150 (MQ=255)
tGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCg                                                                                                                                                     >  1:912892/1‑150 (MQ=255)
 gAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATTTCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGa                                                                                                                                                    >  2:1161317/1‑150 (MQ=255)
   aaacGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAg                                                                                                                                                  >  1:563314/3‑150 (MQ=255)
    cacGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGc                                                                                                                                                 >  1:1166368/1‑150 (MQ=255)
    cacGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGc                                                                                                                                                 >  1:998422/1‑150 (MQ=255)
         tCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgggg                                                                                                                                            <  2:192531/150‑2 (MQ=255)
         tCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgggg                                                                                                                                            <  1:971610/150‑2 (MQ=255)
          cGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgggga                                                                                                                                           <  1:1418728/150‑3 (MQ=255)
            ggATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAgg                                                                                                                                         <  2:146488/150‑1 (MQ=255)
               tGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTca                                                                                                                                      <  1:685117/150‑1 (MQ=255)
                gggTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaa                                                                                                                                     >  2:801834/1‑150 (MQ=255)
                      aCGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGt                                                                                                                               <  1:791806/150‑1 (MQ=255)
                       cGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTg                                                                                                                              >  2:10265/1‑150 (MQ=255)
                        gCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGc                                                                                                                             <  1:1034033/150‑1 (MQ=255)
                         cAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCa                                                                                                                            >  1:803311/1‑150 (MQ=255)
                           aTCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGc                                                                                                                          >  1:93443/1‑150 (MQ=255)
                           aTCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGCTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGc                                                                                                                          <  1:833948/150‑1 (MQ=255)
                               ggCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGGAGCGTaa                                                                                                                      <  2:563314/150‑1 (MQ=255)
                               ggCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTaa                                                                                                                      <  1:52522/150‑1 (MQ=255)
                                 cTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTaaaa                                                                                                                    >  2:677105/1‑150 (MQ=255)
                                          aGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTa                                                                                                           >  2:103162/1‑150 (MQ=255)
                                           gcccgcccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTaa                                                                                                          >  2:787001/1‑150 (MQ=255)
                                                   aGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGCCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCg                                                                                                  >  2:1270735/1‑150 (MQ=255)
                                                    gACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGt                                                                                                 <  1:295090/150‑1 (MQ=255)
                                                      cTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCa                                                                                               <  2:998422/150‑1 (MQ=255)
                                                      cTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCa                                                                                               >  1:269844/1‑150 (MQ=255)
                                                       tCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCaa                                                                                              <  1:962061/150‑1 (MQ=255)
                                                         cAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCATGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATg                                                                                            <  2:1190721/150‑1 (MQ=255)
                                                         cAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATg                                                                                            >  2:721703/1‑150 (MQ=255)
                                                             cTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAAc                                                                                        >  1:1373227/1‑150 (MQ=255)
                                                              tAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGCGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACg                                                                                       >  2:699923/1‑150 (MQ=255)
                                                                aCCAAATGGCGTTATCAGGTGCCTGTTGCAGCACGGCTTCGGACATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGcc                                                                                     <  2:772048/150‑1 (MQ=255)
                                                                  cAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCcac                                                                                   >  2:1562447/1‑150 (MQ=255)
                                                                   aaaTGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCcaca                                                                                  >  1:742532/1‑150 (MQ=255)
                                                                    aaTGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCacac                                                                                 <  1:156524/150‑1 (MQ=255)
                                                                    aaTGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCacac                                                                                 >  2:293797/1‑150 (MQ=255)
                                                                       ggCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTc                                                                              <  1:677105/150‑1 (MQ=255)
                                                                              tCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTc                                                                       >  1:905771/1‑150 (MQ=255)
                                                                              tCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTc                                                                       <  1:14826/150‑1 (MQ=255)
                                                                               cAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCa                                                                      >  1:842503/1‑150 (MQ=255)
                                                                               cAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGCCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCa                                                                      >  2:1467114/1‑150 (MQ=255)
                                                                                aGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAg                                                                     <  2:196491/150‑1 (MQ=255)
                                                                                     ccTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTcccc                                                                <  1:732380/150‑1 (MQ=255)
                                                                                      cTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCa                                                               <  2:1373227/150‑1 (MQ=255)
                                                                                       tGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCAt                                                              >  2:704267/1‑150 (MQ=255)
                                                                                       tGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCAt                                                              >  1:610417/1‑150 (MQ=255)
                                                                                       tGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCAt                                                              >  1:1398643/1‑150 (MQ=255)
                                                                                             agcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgc                                                        <  1:704267/150‑1 (MQ=255)
                                                                                              gcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgca                                                       <  1:103162/150‑1 (MQ=255)
                                                                                              gcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgca                                                       >  2:1377895/1‑150 (MQ=255)
                                                                                                aCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAGTGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagc                                                     <  2:731500/150‑1 (MQ=255)
                                                                                                 cGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagca                                                    <  2:1236152/150‑1 (MQ=255)
                                                                                                   gCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagcagc                                                  <  1:1046703/150‑1 (MQ=255)
                                                                                                    cTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTgcagcagca                                                 <  2:7136/150‑1 (MQ=255)
                                                                                                     ttCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGcac                                                <  1:1407316/150‑1 (MQ=255)
                                                                                                       cGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCACGAAGCCGGGCAGGTCAACAAGCTGAAGCGTAAAATGCGACCTAAGTTCCTCGTCAATGCAGAGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACaa                                              >  2:540302/1‑150 (MQ=255)
                                                                                                        ggCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAg                                             >  1:1358576/1‑150 (MQ=255)
                                                                                                        ggCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGa                                            <  2:40358/150‑1 (MQ=255)
                                                                                                               tCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACa                                      <  2:905771/150‑1 (MQ=255)
                                                                                                               tCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACa                                      <  1:241099/150‑1 (MQ=255)
                                                                                                                    aaGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTccc                                 <  2:912892/150‑1 (MQ=255)
                                                                                                                     aGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTcccc                                <  2:1398643/150‑1 (MQ=255)
                                                                                                                     aGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTcccc                                <  1:787001/150‑1 (MQ=255)
                                                                                                                         aCAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGAc                            <  1:1578370/150‑1 (MQ=255)
                                                                                                                          cAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGAcc                           >  2:587783/1‑150 (MQ=255)
                                                                                                                           aGCGCACCAAATCCCCGGCTACGCCCGAAGCCTGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCt                          >  2:1043241/1‑150 (MQ=255)
                                                                                                                           aGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCt                          >  2:57483/1‑150 (MQ=255)
                                                                                                                               cACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTgg                      <  1:1161317/150‑1 (MQ=255)
                                                                                                                                   acctgcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAACGTGCAGCGTAAAATGCGAGCTAAGTGCCTCGTCAATGCCACGCCACACTTCCGCATTCAGTCCCCATCCGTGCAGCCGCACAAGATGAACAATCCCCTGACCTTTGGTCTg                  >  1:1409826/6‑150 (MQ=255)
                                                                                                                                   aaaTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTg                  <  2:842503/150‑1 (MQ=255)
                                                                                                                                     ctgcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcc                >  1:402992/4‑150 (MQ=255)
                                                                                                                                         ccGGCTACGCCCGAAGCCGGGCAGGTCAACAATGTGCAGCGTAAAATGCGAGCTAAGTTCCTAGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAAGATTCCCCTGACCTTTGGTCTGccacca            <  2:766999/150‑1 (MQ=255)
                                                                                                                                           ggCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACCCTTCGGCATTCAGTCCCCCTCCGTGCAGCAGCACAAGATGCACATTCCCCTGACCTTTGGTCTGCCACCAGa          >  1:1034536/1‑150 (MQ=255)
                                                                                                                                               aCGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGtt      <  2:70153/150‑1 (MQ=255)
                                                                                                                                                cGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTa     >  2:619004/1‑150 (MQ=255)
                                                                                                                                                 gCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGAAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTAt    >  1:121084/1‑150 (MQ=255)
                                                                                                                                                  cccGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTAtt   <  2:1336925/150‑1 (MQ=255)
                                                                                                                                                  cccGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTAtt   >  1:30217/1‑150 (MQ=255)
                                                                                                                                                  cccGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTAtt   >  1:631183/1‑150 (MQ=255)
                                                                                                                                                   ccGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTc  <  1:1383563/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATTC  >  NC_000913/3544546‑3544842

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