BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,559,349	G→T	6.9%	L60I (CTC→ATC)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,559,349	0	G	T	6.9%	455.1 / 14.7	157	L60I (CTC→ATC)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand): ref base G (79/66); new base T (5/7); total (84/73)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

GATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGG  >  NC_000913/4559200‑4559498

gaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCagag                                                                                                                                                       >  2:1086375/1‑150 (MQ=255)

 aTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAgg                                                                                                                                                      >  2:1294729/1‑150 (MQ=255)

 aTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAgg                                                                                                                                                      <  1:1077771/150‑1 (MQ=255)

   agagTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAATTGGACGTGTTGATCAATAAAACCTGGGCAGAGGAt                                                                                                                                                    >  2:410845/1‑150 (MQ=255)

    gagTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGAGCAATAAAACCTCGGCAGAGGAtg                                                                                                                                                   >  2:507436/1‑149 (MQ=255)

     agTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCt                                                                                                                                                  <  1:1735983/150‑1 (MQ=255)

      gTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTg                                                                                                                                                 <  2:843996/150‑1 (MQ=255)

      gTCGGTGCCCAGCAGACCAACCACTGACGTGACGCACGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAAAGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTg                                                                                                                                                 >  1:825720/1‑150 (MQ=255)

         ggTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGTCGAATTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCATCGCCCCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGCTCTGccc                                                                                                                                              >  1:535283/1‑150 (MQ=255)

           tGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGATCCTCTGCCCag                                                                                                                                            <  1:60951/150‑2 (MQ=255)

           ggCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATTGACGTGTTGATCAATAAAACCTGGGCAGATCCTCTGCCCag                                                                                                                                            <  1:85392/149‑2 (MQ=255)

           ggCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGATCCTCTGCCCag                                                                                                                                            <  2:466279/149‑2 (MQ=255)

            gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAACCCTGGGCAGAGGATCTGCCCACt                                                                                                                                           >  1:616400/1‑150 (MQ=255)

            gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGATCCTCTGCCCagg                                                                                                                                           <  2:677642/150‑3 (MQ=255)

            gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACt                                                                                                                                           <  1:41748/150‑1 (MQ=255)

             cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGccacctg                                                                                                                                          >  2:1594177/1‑145 (MQ=255)

              ccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTgg                                                                                                                                         >  1:301022/1‑149 (MQ=255)

              ccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTga                                                                                                                                         >  1:22754/1‑150 (MQ=255)

              ccAGCAGACCAACCACTGACGGGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTga                                                                                                                                         <  1:382829/150‑1 (MQ=255)

                 gcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGAt                                                                                                                                      <  2:1237542/150‑1 (MQ=255)

                  cagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTAAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGGTGATCAATAACACCTGGGCAGAGGATCTGCCCACTGAGATc                                                                                                                                     <  1:738305/150‑1 (MQ=255)

                   agACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGATGATCTGCCACCTGAGATCg                                                                                                                                    >  2:136902/1‑150 (MQ=255)

                         accaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAAcc                                                                                                                              >  2:270056/1‑150 (MQ=255)

                          ccaCAGCCGTGCCGCCCGCTTCCGTCAGGCGACTACGCGCCACCTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCg                                                                                                                             <  2:279150/150‑1 (MQ=255)

                           caCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGcc                                                                                                                                               <  1:1801422/131‑1 (MQ=255)

                            aCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTg                                                                                                                           >  2:382140/1‑150 (MQ=255)

                               gacgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAg                                                                                                                        >  1:999304/1‑150 (MQ=255)

                                 cgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGtt                                                                                                                      >  2:748862/1‑150 (MQ=255)

                                      cgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                                                                 >  1:136006/1‑150 (MQ=255)

                                      cgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                                                                 >  1:871290/1‑150 (MQ=255)

                                      cgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                                                                 >  1:446555/1‑150 (MQ=255)

                                      cgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGGGTGGATCAAGAAAACCTGGGCAGAGGATCTGCCCACTGAGAGCGCCAACCGTGCAGTTCGGTa                                                                                                                 >  1:1045624/1‑150 (MQ=255)

                                       gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTAc                                                                                                                <  1:1608342/150‑1 (MQ=255)

                                       gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTAc                                                                                                                <  2:734066/150‑1 (MQ=255)

                                         ccGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACaa                                                                                                              <  2:1059773/150‑1 (MQ=255)

                                          cGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAAt                                                                                                             <  2:34097/150‑1 (MQ=255)

                                            cTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGt                                                                                                           <  2:586812/150‑1 (MQ=255)

                                             ttCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTc                                                                                                          <  1:642787/150‑1 (MQ=255)

                                                 gTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGaa                                                                                                      >  1:1741235/1‑150 (MQ=255)

                                                      gCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGATGGGat                                                                                                 >  1:1284529/1‑150 (MQ=255)

                                                        gACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGatat                                                                                               >  2:637432/1‑150 (MQ=255)

                                                          cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCTATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                             >  2:923024/1‑150 (MQ=255)

                                                          cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                             >  2:1848392/1‑150 (MQ=255)

                                                          cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                             >  1:264273/1‑150 (MQ=255)

                                                          cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                             >  2:744269/1‑150 (MQ=255)

                                                          cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                             >  2:1855735/1‑150 (MQ=255)

                                                            tAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCt                                                                                           >  2:1438567/1‑150 (MQ=255)

                                                             aGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTg                                                                                          >  1:1484597/1‑150 (MQ=255)

                                                              gcgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGACAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTgg                                                                                         >  1:1383127/1‑150 (MQ=255)

                                                               cgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGc                                                                                        >  1:572406/1‑150 (MQ=255)

                                                               cgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGc                                                                                        >  1:454006/1‑150 (MQ=255)

                                                               cgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGc                                                                                        >  2:1183725/1‑150 (MQ=255)

                                                                     ttCCGGCGTTAGCGGCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGCTCTGACCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGc                                                                                  <  2:154233/150‑1 (MQ=255)

                                                                     ttCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTCCAAGGTCAGAAGGGATATTGCTGGCAACggg                                                                                  >  1:1709944/1‑149 (MQ=255)

                                                                     ttCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGc                                                                                  >  2:1500927/1‑150 (MQ=255)

                                                                         ggCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATa                                                                              >  1:972127/1‑150 (MQ=255)

                                                                               cgcgTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTg                                                                        >  2:1295644/1‑150 (MQ=255)

                                                                                  gtcgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCg                                                                     <  2:109543/150‑1 (MQ=255)

                                                                                  gtcgtGGGACCCGCTTCGCCGCCACCGCCAATCAAATGGCCGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCg                                                                     <  1:1006238/150‑1 (MQ=255)

                                                                                   tcgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGt                                                                    >  2:809788/1‑150 (MQ=255)

                                                                                      ggggACCTGCTGCGCCGCCACCGCCAATCAAATGGACGGGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAg                                                                 <  2:147958/149‑1 (MQ=255)

                                                                                         gACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcga                                                              >  2:875395/1‑150 (MQ=255)

                                                                                          aCCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgac                                                             >  1:1027908/1‑150 (MQ=255)

                                                                                           ccTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgacg                                                            >  2:1221472/1‑150 (MQ=255)

                                                                                              gCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAgg                                                         <  1:748862/150‑1 (MQ=255)

                                                                                              gCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAgg                                                         >  2:1281508/1‑150 (MQ=255)

                                                                                              gCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGGTAGCGACGAgg                                                         <  2:1637271/150‑1 (MQ=255)

                                                                                                      gccaccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCa                                                 <  2:1715472/150‑1 (MQ=255)

                                                                                                        caccgccaTTCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATATTTTTGCCGTTAGCGACGAGGACATCGCaaa                                               >  2:1544197/1‑150 (MQ=255)

                                                                                                          ccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAAtt                                             >  1:165450/1‑150 (MQ=255)

                                                                                                            gccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTcc                                           >  1:1479464/1‑150 (MQ=255)

                                                                                                             ccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTccc                                          <  1:1252579/150‑1 (MQ=255)

                                                                                                              caATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTcccc                                         <  1:1086375/150‑1 (MQ=255)

                                                                                                                 tCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGAt                                      <  1:1257802/150‑1 (MQ=255)

                                                                                                                  cAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATc                                     >  1:44125/1‑150 (MQ=255)

                                                                                                                  cAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATc                                     <  2:34800/150‑1 (MQ=255)

                                                                                                                  cAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAAGTCCCCGATc                                     >  2:715191/1‑150 (MQ=255)

                                                                                                                  aaaatttaCGTTTTGATAAATAAAACCTGGTCAGATGATCTGCCCTCTGAGATCGACAACCGTTCAGTTCTGTACAATGTCAGAAGTGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATc                                     <  1:1662631/143‑1 (MQ=255)

                                                                                                                      tGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTGCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTTGCGACGAGGACATCGCAAATTCCCCGATCTTcc                                 >  2:916613/1‑150 (MQ=255)

                                                                                                                       ggACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCg                                >  1:650528/1‑150 (MQ=255)

                                                                                                                          cGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGgcg                             <  2:1716013/150‑1 (MQ=255)

                                                                                                                                  tCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACaaa                     >  1:1362330/1‑150 (MQ=255)

                                                                                                                                    aaTAAAACCTGGGCAGATGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATg                   <  2:1632294/150‑1 (MQ=255)

                                                                                                                                    aaTAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTTCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATg                   <  1:915639/150‑1 (MQ=255)

                                                                                                                                      tAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATgcg                 <  2:891928/150‑1 (MQ=255)

                                                                                                                                         tcagtGGGCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGCAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTc              >  1:1505116/5‑150 (MQ=255)

                                                                                                                                          aCCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTcc             <  2:795954/150‑1 (MQ=255)

                                                                                                                                           ccTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTccc            >  1:787977/1‑150 (MQ=255)

                                                                                                                                           ccTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTccc            <  2:301022/150‑1 (MQ=255)

                                                                                                                                           ccTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCATCCGGCGCATACAAATGCGCTccc            <  1:1177819/150‑1 (MQ=255)

                                                                                                                                             tGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCtt          >  2:1838518/1‑149 (MQ=255)

                                                                                                                                             tGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTg          <  2:210215/150‑1 (MQ=255)

                                                                                                                                              gggCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgc         >  2:411682/1‑150 (MQ=255)

                                                                                                                                               ggCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgca        <  2:454006/150‑1 (MQ=255)

                                                                                                                                                gCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgcag       >  1:1602866/1‑150 (MQ=255)

                                                                                                                                                 cAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgcagc      <  2:1573255/150‑1 (MQ=255)

                                                                                                                                                 cAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTGCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgcagc      <  1:270056/150‑1 (MQ=255)

                                                                                                                                                  agagGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgcagca     <  1:666973/150‑1 (MQ=255)

                                                                                                                                                    agGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAgg   >  2:1800336/1‑150 (MQ=255)

                                                                                                                                                     ggATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAggg  <  1:1594177/150‑1 (MQ=255)

GATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGG  >  NC_000913/4559200‑4559498

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