BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,828,334	A→T	5.5%	I111N (ATT→AAT)	gltS ←	glutamate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,828,334	0	A	T	5.5%	495.4 / 12.9	142	I111N (ATT→AAT)	gltS	glutamate transporter
Reads supporting (aligned to +/- strand): ref base A (64/70); new base T (3/5); total (67/75)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 161 total aligned reads displayed.

ACTTCCGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAA  >  NC_000913/3828187‑3828481

aCTTCCGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATAccaatgccaa                                                                                                                                                     >  3:289086/1‑148 (MQ=255)

    ccGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGc                                                                                                                                                 >  4:336423/1‑148 (MQ=255)

    ccGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCa                                                                                                                                                <  2:373596/149‑1 (MQ=255)

      gTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCAtt                                                                                                                                              <  2:323000/149‑1 (MQ=255)

        cGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCAtttt                                                                                                                                            >  1:12436/1‑149 (MQ=255)

          cATTGGTGAAGCCATAACGTACAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCTGCGGATCAAGCCCTAACAAGCAAGCCATACCAATGCCAATGGCATTTTg                                                                                                                                           <  4:1506/148‑1 (MQ=255)

           aTTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCATTGGCAttggta                                                                                                                                          <  4:295761/148‑5 (MQ=255)

           aTTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCATTGGCAttggta                                                                                                                                          <  4:16005/148‑5 (MQ=255)

            ttGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCGGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCa                                                                                                                                         <  4:56145/148‑1 (MQ=255)

                tGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCAc                                                                                                                                     <  3:199020/148‑1 (MQ=255)

                       aTAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGAc                                                                                                                             <  1:339824/149‑1 (MQ=255)

                         aaCGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGAcc                                                                                                                            >  3:376381/1‑148 (MQ=255)

                              tCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCaacaac                                                                                                                      <  1:396285/149‑1 (MQ=255)

                                 aTGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACca                                                                                                                    >  4:113958/1‑148 (MQ=255)

                                  tGAGCAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACcaca                                                                                                                  >  2:10404/1‑149 (MQ=255)

                                  tGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACcaca                                                                                                                  >  1:366415/1‑149 (MQ=255)

                                   gAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACaa                                                                                                                 <  1:180307/149‑1 (MQ=255)

                                       aaTTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCAcc                                                                                                                                    >  2:61600/1‑126 (MQ=255)

                                       aaTTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCaa                                                                                                             >  1:118733/1‑149 (MQ=255)

                                       aaTTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCaa                                                                                                             <  1:205786/149‑1 (MQ=255)

                                        aTTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCaa                                                                                                             >  3:40072/1‑148 (MQ=255)

                                          ttACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCaaga                                                                                                           >  3:201059/1‑148 (MQ=255)

                                              tCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATg                                                                                                      <  1:367733/149‑1 (MQ=255)

                                                 aCGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGcc                                                                                                    >  3:322819/1‑148 (MQ=255)

                                                 aCGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCa                                                                                                   <  2:12436/149‑1 (MQ=255)

                                                  cGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCa                                                                                                   >  3:109916/1‑148 (MQ=255)

                                                    cAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAAc                                                                                                 >  3:387631/1‑148 (MQ=255)

                                                       cgcCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACcac                                                                                              <  3:391931/148‑1 (MQ=255)

                                                         cccGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCaaga                                                                                                           >  2:119936/1‑133 (MQ=255)

                                                             tACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGccc                                                                                        >  4:395491/1‑148 (MQ=255)

                                                              aCCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGAc                                                                                                                             >  1:132117/1‑110 (MQ=255)

                                                              aCCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGAc                                                                                                                             <  2:132117/110‑1 (MQ=255)

                                                                cGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCg                                                                                    >  1:152233/1‑149 (MQ=255)

                                                                 gTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCgg                                                                                   >  2:101424/1‑149 (MQ=255)

                                                                    aCCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCa                                                                                 <  4:309028/148‑1 (MQ=255)

                                                                     ccgccgGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCATTGGCAttgg                                                                                                                                            <  4:357279/88‑3 (MQ=255)

                                                                     ccgccgGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCATTGGCAttgg                                                                                                                                            >  3:357279/1‑86 (MQ=255)

                                                                     ccgccgGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACg                                                                               >  2:228279/1‑149 (MQ=255)

                                                                          ggAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAAc                                                                          >  1:296256/1‑149 (MQ=255)

                                                                             aaGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTgg                                                                       <  2:305525/149‑1 (MQ=255)

                                                                             aaGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTg                                                                        <  3:113958/148‑1 (MQ=255)

                                                                                 gTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAAt                                                                   <  2:66328/149‑1 (MQ=255)

                                                                                  tAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAAc                                                                                                                                 <  3:119669/86‑1 (MQ=255)

                                                                                    aTAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGtt                                                                >  2:63097/1‑149 (MQ=255)

                                                                                    aTAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGtt                                                                <  1:101424/149‑1 (MQ=255)

                                                                                     tAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGtt                                                                >  4:21788/1‑148 (MQ=255)

                                                                                     tAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTg                                                               <  1:177919/149‑1 (MQ=255)

                                                                                     gaGAACCGGCCAACATCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGtt                                                                <  3:16205/147‑1 (MQ=255)

                                                                                      aGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGTCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTg                                                               >  4:224013/1‑148 (MQ=255)

                                                                                       gAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTgg                                                              <  4:327517/148‑1 (MQ=255)

                                                                                        aaCCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCg                                                            >  2:282693/1‑149 (MQ=255)

                                                                                            ggCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACATGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTc                                                         <  4:383081/148‑1 (MQ=255)

                                                                                            ggCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCa                                                        <  1:63097/149‑1 (MQ=255)

                                                                                             gCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCa                                                        >  3:86393/1‑148 (MQ=255)

                                                                                               cAACAGCCCCGTCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAgg                                                      <  4:275444/148‑1 (MQ=255)

                                                                                               cAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGTCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCaa                                                                    <  2:20414/134‑1 (MQ=255)

                                                                                               cAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGTCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCaa                                                                    >  1:20414/1‑134 (MQ=255)

                                                                                                       ccATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGTGTTCAGGCCAAtggt                                              <  3:336423/148‑1 (MQ=255)

                                                                                                        cATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAAtggtgg                                            >  2:39146/1‑149 (MQ=255)

                                                                                                        cATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAAtggtgg                                            >  1:73457/1‑149 (MQ=255)

                                                                                                          tCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCaagaa                                                                                                          <  2:46466/85‑1 (MQ=255)

                                                                                                          tCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGc                                           >  4:200226/1‑148 (MQ=255)

                                                                                                             gCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCgaag                                       >  1:76358/1‑149 (MQ=255)

                                                                                                                gATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCgaagaa                                     <  4:376381/148‑1 (MQ=255)

                                                                                                                 aTCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGaaa                                    <  4:105242/148‑1 (MQ=255)

                                                                                                                 aTCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAg                                   <  2:73457/149‑1 (MQ=255)

                                                                                                                  tCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAg                                   >  4:395773/1‑148 (MQ=255)

                                                                                                                  tCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAg                                   <  4:19149/148‑1 (MQ=255)

                                                                                                                    aaGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCa                                <  2:366415/149‑1 (MQ=255)

                                                                                                                    aaGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCa                                <  2:15158/149‑1 (MQ=255)

                                                                                                                      gCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGc                              <  1:202134/149‑1 (MQ=255)

                                                                                                                       cccTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGt                                                                 >  1:227588/1‑113 (MQ=255)

                                                                                                                       cccTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCa                             >  2:144068/1‑149 (MQ=255)

                                                                                                                         cTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACa                                                                                                                                <  2:354473/48‑1 (MQ=255)

                                                                                                                         cTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACa                                                                                                                                >  1:354473/1‑48 (MQ=255)

                                                                                                                         cTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCAtt                           >  1:54840/1‑149 (MQ=255)

                                                                                                                         cTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCAt                            >  3:107727/1‑148 (MQ=255)

                                                                                                                           aaCAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTa                          <  3:395773/148‑1 (MQ=255)

                                                                                                                            aCAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGcc                                 >  3:79115/1‑140 (MQ=255)

                                                                                                                            aCAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGcc                                 <  4:79115/140‑1 (MQ=255)

                                                                                                                            aCAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAAc                        <  2:169221/149‑1 (MQ=255)

                                                                                                                               aGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGTCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACgg                      <  3:224013/148‑1 (MQ=255)

                                                                                                                               aGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACgg                      <  4:109916/148‑1 (MQ=255)

                                                                                                                                     ccATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGGCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCgcgc               >  1:395942/1‑149 (MQ=255)

                                                                                                                                      cATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCgcgc               <  4:201059/148‑1 (MQ=255)

                                                                                                                                        tACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAgg            >  1:218709/1‑149 (MQ=255)

                                                                                                                                         aCCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAggg           >  1:296296/1‑149 (MQ=255)

                                                                                                                                          ccaatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCaaa                                                                           <  1:350009/84‑1 (MQ=255)

                                                                                                                                          ccaatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGAGGCCAACCACACGCCCACCGGCACGCaaa                                                                           >  2:350009/1‑84 (MQ=255)

                                                                                                                                           caatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGc                                  >  2:205102/1‑124 (MQ=255)

                                                                                                                                           caatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGAc         >  1:48657/1‑149 (MQ=255)

                                                                                                                                           caatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGAc         >  1:374411/1‑149 (MQ=255)

                                                                                                                                           caatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGAc         >  1:258829/1‑149 (MQ=255)

                                                                                                                                            aatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGAc         >  4:167616/1‑148 (MQ=255)

                                                                                                                                            aatgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACa        >  2:287660/1‑149 (MQ=255)

                                                                                                                                             atgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACat       >  2:88458/1‑149 (MQ=255)

                                                                                                                                             atgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACat       <  1:222813/149‑1 (MQ=255)

                                                                                                                                              tgccaatgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACata      <  1:141178/149‑1 (MQ=255)

                                                                                                                                                   atgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCaa  >  4:77601/1‑148 (MQ=255)

                                                                                                                                                   atgGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCaa  <  4:86393/148‑1 (MQ=255)

ACTTCCGTCGCATTGGTGAAGCCATAACGTTCAATGAACAATTTACTCCACGCAGCGCCCGTACCGTGACCGCCGGAAAGAGTAATAGAACCGGCCAACAGCCCCATCAGCGGATCAAGCCCTAACAAGCTAGCCATACCAATGCCAATGGCATTTTGCATCACCAACAGACCAACAACCACAATCAAGAAGATGCCAACCACACGCCCACCGGCACGCAAACTGGCAATGTTGGCGTTCAGGCCAATGGTGGCGAAGAAAGCCAGCATTAACGGATCGCGCAGGGACATATCAA  >  NC_000913/3828187‑3828481

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