BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	242,344	C→G	5.2%	P320P (CCG→CCC)	fadE ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	242,344	0	C	G	5.2%	329.9 / 6.5	96	P320P (CCG→CCC)	fadE	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand): ref base C (44/47); new base G (3/2); total (47/49)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

GGCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCA  >  NC_000913/242196‑242486

ggCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAACACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCg                                                                                                                                               >  2:1253949/1‑150 (MQ=255)

  cACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGc                                                                                                                                             >  1:998017/1‑150 (MQ=255)

  cACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGc                                                                                                                                             >  2:255373/1‑150 (MQ=255)

    cTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgt                                                                                                                                           >  2:974553/1‑150 (MQ=255)

      ccaccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCCCGCCCGGCgtgg                                                                                                                                         >  2:1502601/1‑150 (MQ=255)

      ccaccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtgg                                                                                                                                         >  2:449182/1‑150 (MQ=255)

         ccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGtggtgg                                                                                                                                      >  2:1191042/1‑150 (MQ=255)

         ccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGtggtgg                                                                                                                                      >  1:1494473/1‑150 (MQ=255)

           aGCATCCGCCAGCCTTGCCCGTCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaatt                                                                                                                                    <  1:206178/150‑6 (MQ=255)

           aGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaatt                                                                                                                                    <  1:1314108/150‑6 (MQ=255)

                 cGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATc                                                                                                                              >  2:456035/1‑150 (MQ=255)

                   ccAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGGCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAg                                                                                                                            <  1:1368545/150‑1 (MQ=255)

                      gCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAgcgc                                                                                                                         <  2:200458/150‑1 (MQ=255)

                      gCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAgcgc                                                                                                                         >  1:633761/1‑150 (MQ=255)

                        cTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGcac                                                                                                                       <  2:951632/150‑1 (MQ=255)

                            cccGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGt                                                                                                                   >  2:1274444/1‑150 (MQ=255)

                            cccGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGt                                                                                                                   >  2:1233739/1‑150 (MQ=255)

                                 ccATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGc                                                                                                              <  2:1025975/150‑1 (MQ=255)

                                  cATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGcc                                                                                                             >  1:631061/1‑150 (MQ=255)

                                  cATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGcc                                                                                                             <  1:753979/150‑1 (MQ=255)

                                    ttttCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTa                                                                                                           >  1:1168915/1‑150 (MQ=255)

                                      ttCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGTATCAGCGCACAGGTAATGCCTaaa                                                                                                         <  1:255373/150‑1 (MQ=255)

                                           cccGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATc                                                                                                                              >  1:306647/1‑124 (MQ=255)

                                           cccGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATc                                                                                                                              <  2:306647/124‑1 (MQ=255)

                                           cccGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtcttc                                                                                                    <  1:691743/150‑1 (MQ=255)

                                             cgccgATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTg                                                                                                  <  1:1274444/150‑1 (MQ=255)

                                                 gatgatGTAATCGATCGGCACGAAGACGTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAcc                                                                                              <  2:301598/150‑1 (MQ=255)

                                                 gatgatGTAAGCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAcc                                                                                              <  2:124580/150‑1 (MQ=255)

                                                     atgTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACACCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAccgccg                                                                                          >  1:42813/1‑150 (MQ=255)

                                                        tAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGt                                                                                       <  2:304319/150‑1 (MQ=255)

                                                            cgatcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAtt                                                                                   >  2:1379104/1‑150 (MQ=255)

                                                             gatcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAttt                                                                                  >  2:1597/1‑150 (MQ=255)

                                                                  gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCg                                                                             <  1:897899/150‑1 (MQ=255)

                                                                  gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCg                                                                             <  1:892180/150‑1 (MQ=255)

                                                                  gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCg                                                                             <  1:1356641/150‑1 (MQ=255)

                                                                  gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCg                                                                             >  2:1609092/1‑150 (MQ=255)

                                                                     cGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGt                                                                          >  1:1299619/1‑150 (MQ=255)

                                                                      gAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTc                                                                         >  2:1277191/1‑150 (MQ=255)

                                                                         gACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGga                                                                      <  2:1168915/150‑1 (MQ=255)

                                                                          aCATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgag                                                                     <  2:1501265/150‑1 (MQ=255)

                                                                           cATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                    >  2:34149/1‑150 (MQ=255)

                                                                              cTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGtt                                                                 >  2:889523/1‑150 (MQ=255)

                                                                              cTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGtt                                                                 >  2:1204523/1‑150 (MQ=255)

                                                                              cTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAACTTTTCCGGGTCGGAGAGtt                                                                 >  1:882064/1‑150 (MQ=255)

                                                                              cTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAACTTTTCCGGGTCGGAGAGtt                                                                 <  2:882064/150‑1 (MQ=255)

                                                                               tttACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGttt                                                                <  1:91193/150‑1 (MQ=255)

                                                                               tttACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGttt                                                                <  1:85302/150‑1 (MQ=255)

                                                                                 tACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTaa                                                              >  2:921722/1‑150 (MQ=255)

                                                                                  aCCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTaaa                                                             >  1:963671/1‑150 (MQ=255)

                                                                                   ccGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAccg                                                                                             >  1:1403758/1‑117 (MQ=255)

                                                                                   ccGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAccg                                                                                             <  2:1403758/117‑1 (MQ=255)

                                                                                   ccGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAAc                                                            <  1:1379104/150‑1 (MQ=255)

                                                                                    cgcgcgTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACg                                                           <  1:1191042/150‑1 (MQ=255)

                                                                                          tCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGcc                                                     >  1:403063/1‑150 (MQ=255)

                                                                                            ggTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCa                                                   >  2:536095/1‑150 (MQ=255)

                                                                                            ggTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCa                                                   <  2:1494473/150‑1 (MQ=255)

                                                                                             gTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCaa                                                  >  2:780462/1‑150 (MQ=255)

                                                                                               ccGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGc                                                >  2:1275026/1‑150 (MQ=255)

                                                                                                    cTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGt                                           >  2:436013/1‑150 (MQ=255)

                                                                                                    cTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGt                                           <  2:998017/150‑1 (MQ=255)

                                                                                                    cTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGt                                           <  1:1123969/150‑1 (MQ=255)

                                                                                                    cTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCAGAGAGTTTAAACGCCAGCCCAAGCACGGt                                           >  2:204918/1‑150 (MQ=255)

                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCg                                         >  2:385419/1‑150 (MQ=255)

                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCg                                         <  1:13340/150‑1 (MQ=255)

                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCg                                         >  2:71928/1‑150 (MQ=255)

                                                                                                         aCGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCaa                                      <  2:1321138/150‑1 (MQ=255)

                                                                                                            gTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCg                                   >  2:341072/1‑150 (MQ=255)

                                                                                                              aCGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGt                                 <  1:79511/150‑1 (MQ=255)

                                                                                                               cGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTg                                <  1:151706/150‑1 (MQ=255)

                                                                                                               cGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTg                                <  2:589992/150‑1 (MQ=255)

                                                                                                                gTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGc                               <  1:974553/150‑1 (MQ=255)

                                                                                                                 ttCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGagag                                                                   <  1:694240/113‑1 (MQ=255)

                                                                                                                 ttCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGagag                                                                   >  2:694240/1‑113 (MQ=255)

                                                                                                                   cAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAg                            <  2:407544/150‑1 (MQ=255)

                                                                                                                   cAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAg                            <  2:1476784/150‑1 (MQ=255)

                                                                                                                   cAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACCGTCGCAATCGGTGCCAg                            >  2:1599832/1‑150 (MQ=255)

                                                                                                                     gCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCg                          <  1:449182/150‑1 (MQ=255)

                                                                                                                      cGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGt                         >  2:1516686/1‑150 (MQ=255)

                                                                                                                       gggAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTa                        >  2:1476976/1‑150 (MQ=255)

                                                                                                                        ggAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTaa                       <  1:1275305/150‑1 (MQ=255)

                                                                                                                          aaGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATg                     <  2:1330813/150‑1 (MQ=255)

                                                                                                                            gTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTa                   >  2:997323/1‑150 (MQ=255)

                                                                                                                                cgacgaCCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGt               <  1:5683/150‑1 (MQ=255)

                                                                                                                                    gacCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGt           <  1:1246170/150‑1 (MQ=255)

                                                                                                                                      ccaaccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc         >  1:406017/8‑150 (MQ=255)

                                                                                                                                      ccAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc         >  1:368170/1‑150 (MQ=255)

                                                                                                                                      ccAATTTCCACGCCCGGCGTGGTGGTTGGAATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc         >  2:80220/1‑150 (MQ=255)

                                                                                                                                       caaccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc        >  1:1075172/7‑150 (MQ=255)

                                                                                                                                       caaccaccACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc        <  1:1200040/144‑1 (MQ=255)

                                                                                                                                       cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc        <  1:1521388/150‑1 (MQ=255)

                                                                                                                                       cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc        <  1:171297/150‑1 (MQ=255)

                                                                                                                                        aaTTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCa       <  1:1476976/150‑1 (MQ=255)

                                                                                                                                         aTTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg      <  1:1609092/150‑1 (MQ=255)

                                                                                                                                          tttCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAgg     <  2:631061/150‑1 (MQ=255)

                                                                                                                                          ccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAgg     >  2:1219023/4‑150 (MQ=255)

                                                                                                                                             ccACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCa  <  1:1253949/150‑1 (MQ=255)

GGCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCA  >  NC_000913/242196‑242486

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: `‑`