BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,107,710	C→T	8.9%	A452A (GCC→GCT)	clsC →	cardiolipin synthase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,107,710	0	C	T	8.9%	294.7 / 9.7	100	A452A (GCC→GCT)	clsC	cardiolipin synthase C
Reads supporting (aligned to +/- strand): ref base C (47/43); new base T (2/7); total (49/51)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GAGAGCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACGG  >  NC_000913/1107561‑1107855

gagagCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGcc                                                                                                                                                   <  2:1198224/150‑1 (MQ=255)

gagagCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGcc                                                                                                                                                   <  2:143413/150‑1 (MQ=255)

  gagCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGccac                                                                                                                                                 >  2:447682/1‑150 (MQ=255)

  gagCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGccac                                                                                                                                                 >  2:442560/1‑150 (MQ=255)

   agCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGccacc                                                                                                                                                <  2:103962/150‑1 (MQ=255)

   agCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGccacc                                                                                                                                                <  1:746184/150‑1 (MQ=255)

      gAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGt                                                                                                                                             >  1:25003/1‑150 (MQ=255)

      gAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGt                                                                                                                                             >  1:414243/1‑150 (MQ=255)

      gAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATCATGCGGCCTCGCAGCTCCGTCTGGACAGGTGGGGACTGATCAACTCGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAgg                                                                                                                                             >  1:813375/1‑149 (MQ=255)

       aaaCGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGtt                                                                                                                                            <  2:1275899/150‑1 (MQ=255)

         aggCTGGCACAGTTTATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTc                                                                                                                                          <  2:859858/148‑1 (MQ=255)

          cGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCt                                                                                                                                         <  2:796564/150‑1 (MQ=255)

              ggCACAGTTAATTGATAAACGCTTTTTTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGaa                                                                                                                                     >  2:247319/1‑150 (MQ=255)

                   aGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCggg                                                                                                                                <  1:1028934/150‑1 (MQ=255)

                    gTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGt                                                                                                                               >  1:425814/1‑150 (MQ=255)

                      tAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTa                                                                                                                             <  2:836789/150‑1 (MQ=255)

                       aaTTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttct                                                                                                                            <  2:839956/150‑3 (MQ=255)

                       aaTTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACACGCTTCCAGAAACTGGTGGCGGGttct                                                                                                                            <  2:1074675/150‑3 (MQ=255)

                        aTTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttctt                                                                                                                           <  1:555425/150‑4 (MQ=255)

                           gATACACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGTCAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                        <  2:463792/150‑7 (MQ=255)

                           gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGTCAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                        <  2:468916/150‑7 (MQ=255)

                           gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                        <  1:586821/150‑7 (MQ=255)

                           gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCTTCCAGAAACTGGTGGCGGGttcttttt                                                                                                                        <  2:436751/150‑7 (MQ=255)

                           gATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTc                                                                                                                        >  2:717751/1‑150 (MQ=255)

                                        ttCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGa                                                                                                           <  1:1151870/150‑1 (MQ=255)

                                         tCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGat                                                                                                          >  1:1100745/1‑150 (MQ=255)

                                          cAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGata                                                                                                         >  2:406043/1‑150 (MQ=255)

                                          cAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGata                                                                                                         >  2:337650/1‑150 (MQ=255)

                                                cAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGccc                                                                                                   >  1:64696/1‑150 (MQ=255)

                                                 aGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCggg                                                                                                                                >  1:15404/1‑120 (MQ=255)

                                                 aGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCggg                                                                                                                                <  2:15404/120‑1 (MQ=255)

                                                 aGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGATGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCg                                                                                                  >  1:1193979/1‑150 (MQ=255)

                                                    atgatgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTgg                                                                                               <  1:1526520/150‑1 (MQ=255)

                                                     tgatgCGGCCTGGCAGCTCCGTCTGGACAGGGGGGGACGGATCAACTGGGTTGATCGGCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGa                                                                                              >  2:926600/1‑150 (MQ=255)

                                                      gatgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGaa                                                                                             >  2:1347181/1‑150 (MQ=255)

                                                      gatgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGaa                                                                                             >  2:1389521/1‑150 (MQ=255)

                                                        tgCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATg                                                                                           <  2:1524066/150‑1 (MQ=255)

                                                          cGGCCTGGCAGCTCGGTGTGGACAGGTGGGGACGTATCAAGTGTTTTGATCGTCATGCTATGAAAGTGATTATTCTCTAAAAAGAACCCGCGACAAGTTTCTGGAAGAGGGTTATTGTCAGACTGGCGTCGATATTGCCGGTGGAATTGt                                                                                         <  2:421598/150‑1 (MQ=255)

                                                          cGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGt                                                                                         <  2:181662/150‑1 (MQ=255)

                                                           ggCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGtt                                                                                        >  1:1168541/1‑150 (MQ=255)

                                                             ccTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTAt                                                                                      >  2:24928/1‑150 (MQ=255)

                                                             ccTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTAt                                                                                      <  1:510239/150‑1 (MQ=255)

                                                              cTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTAtt                                                                                     <  1:59705/150‑1 (MQ=255)

                                                                 gCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTa                                                                                  <  2:351816/150‑1 (MQ=255)

                                                                  cAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTaa                                                                                 >  1:674201/1‑150 (MQ=255)

                                                                   aGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTaaa                                                                                <  2:425814/150‑1 (MQ=255)

                                                                    gCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAAt                                                                               >  1:1161829/1‑150 (MQ=255)

                                                                       ccGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAAtgtg                                                                            >  2:931305/1‑150 (MQ=255)

                                                                        cGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTgg                                                                           <  1:406043/150‑1 (MQ=255)

                                                                                  ggTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGc                                                                 >  1:175826/1‑150 (MQ=255)

                                                                                  ggTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGc                                                                 >  1:451556/1‑150 (MQ=255)

                                                                                     ggggACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTgg                                                              >  2:1530244/1‑150 (MQ=255)

                                                                                          cGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCtt                                                         >  2:312167/1‑150 (MQ=255)

                                                                                           ggATCAACTGGGTTGATCGTGATGCAAAGTATTAGATTATTCTGTAAATAGTACCGGCCACCAGTTTCTGGAAGCGGGTTATTGTCAGACTGGCGTCGATATTGCGCGTGGATTGGTTATTGTAAATGTGTCTTAACGGGCTGGTGGttt                                                        <  2:918169/150‑1 (MQ=255)

                                                                                            gATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTa                                                       >  2:510868/1‑150 (MQ=255)

                                                                                            gATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTa                                                       <  1:24928/150‑1 (MQ=255)

                                                                                               cAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCg                                                    <  2:969659/150‑1 (MQ=255)

                                                                                                aaCTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGt                                                   >  1:254669/1‑150 (MQ=255)

                                                                                                aaCTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGt                                                   >  1:768311/1‑150 (MQ=255)

                                                                                                     ggTTGATCGTCATGCAGAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCa                                              <  2:1091747/150‑1 (MQ=255)

                                                                                                      gTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAATAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGGCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCAc                                             <  2:255058/150‑1 (MQ=255)

                                                                                                      gTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCAc                                             >  1:1136980/1‑150 (MQ=255)

                                                                                                          aTCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCtt                                         >  1:1525443/1‑150 (MQ=255)

                                                                                                             gTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaac                                      <  1:447682/150‑1 (MQ=255)

                                                                                                             gTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaac                                      <  1:442560/150‑1 (MQ=255)

                                                                                                                aTGCAAAGAAAGTGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaacaac                                   <  1:1390920/150‑1 (MQ=255)

                                                                                                                aTGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGTAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTTGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTaacaac                                   <  1:1396557/150‑1 (MQ=255)

                                                                                                                 tGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAAcc                                  <  1:1389521/150‑1 (MQ=255)

                                                                                                                 tGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAAcc                                  >  2:707308/1‑150 (MQ=255)

                                                                                                                  gCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCg                                 <  1:717751/150‑1 (MQ=255)

                                                                                                                   cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCg                                                    <  2:889709/130‑1 (MQ=255)

                                                                                                                   cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCg                                                    >  1:889709/1‑130 (MQ=255)

                                                                                                                   cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCgg                                >  1:224893/1‑150 (MQ=255)

                                                                                                                   cAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCgg                                <  1:295534/150‑1 (MQ=255)

                                                                                                                    aaagaaagAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGt                               >  2:1184448/1‑150 (MQ=255)

                                                                                                                    aaagaaagAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGt                               <  2:1168541/150‑1 (MQ=255)

                                                                                                                      agaaagaGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGttt                             <  2:674201/150‑1 (MQ=255)

                                                                                                                         aagaGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGGTTAACAACCGGTTTAcc                          >  2:1281358/1‑150 (MQ=255)

                                                                                                                          agatATTATTCTCAACAAAGAACCCGCCTCCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTTGCGTACCGGGCTGGTGCTTTATCGTTCTGCCGCTTAACAACCGGTTTACCa                         <  2:813375/150‑1 (MQ=255)

                                                                                                                          agagATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCa                         >  2:1077285/1‑150 (MQ=255)

                                                                                                                           gagaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                        >  1:783326/1‑150 (MQ=255)

                                                                                                                           gagaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                        <  1:1132095/150‑1 (MQ=255)

                                                                                                                             tattattCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                      <  2:974262/149‑1 (MQ=255)

                                                                                                                             gaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                      <  1:121147/150‑1 (MQ=255)

                                                                                                                             gaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                      <  1:1532119/150‑1 (MQ=255)

                                                                                                                             gaTTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCagag                      <  2:25003/150‑1 (MQ=255)

                                                                                                                                  ttCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACaa                 >  1:1426347/1‑150 (MQ=255)

                                                                                                                                    ctcAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTACCGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACaaaa               >  1:791520/1‑150 (MQ=255)

                                                                                                                                    ctcAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACaaaa               >  1:1269805/1‑150 (MQ=255)

                                                                                                                                      cAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAAc             >  2:1245960/1‑150 (MQ=255)

                                                                                                                                      cAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAAc             >  1:548383/1‑150 (MQ=255)

                                                                                                                                       aaaaaaGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTACCGGGATGGTGCTTTATCGTTCTCACGATTAACAACCGGTTTACCAGAGAACAAAAACt            <  1:707308/150‑1 (MQ=255)

                                                                                                                                         aaaaGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTc          <  1:247319/150‑1 (MQ=255)

                                                                                                                                          ccataaCCCGCTTCCAGAAACTGGTGGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCa         >  1:1374988/5‑150 (MQ=255)

                                                                                                                                          ccataaCCCGCTTCCAGAAACTGGTGGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCa         >  1:1214353/5‑150 (MQ=255)

                                                                                                                                             gAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTa      >  1:1445145/1‑150 (MQ=255)

                                                                                                                                               aCCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                        <  2:462739/130‑1 (MQ=255)

                                                                                                                                               aCCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCag                        >  1:462739/1‑130 (MQ=255)

                                                                                                                                                cccGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACg   >  1:205512/1‑150 (MQ=255)

                                                                                                                                                 ccGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACgg  >  2:524583/1‑150 (MQ=255)

GAGAGCGAAACGCTGGCACAGTTAATTGATAAACGCTTTATTCAGAGCCAGTATGATGCGGCCTGGCAGCTCCGTCTGGACAGGTGGGGACGGATCAACTGGGTTGATCGTCATGCAAAGAAAGAGATTATTCTCAAAAAAGAACCCGCCACCAGTTTCTGGAAGCGGGTTATGGTCAGACTGGCGTCGATATTGCCCGTGGAATGGTTATTGTAAATGTGGCTTAACGGGCTGGTGCTTTATCGTTCTCACGCTTAACAACCGGTTTACCAGAGAACAAAAACTTCAGTAACGG  >  NC_000913/1107561‑1107855

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