BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,559,349	G→T	5.4%	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	5.4%	445.2 / 8.3	129	L60I (CTC→ATC)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand): ref base G (51/71); new base T (2/5); total (53/76)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 131 total aligned reads displayed.

AGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGG  >  NC_000913/4559205‑4559498

 gTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGGCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTg                                                                                                                                                 <  2:1520035/150‑1 (MQ=255)

 gTCGGTGCACAGCAGACCAACCACTGATGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTg                                                                                                                                                 <  1:928893/150‑1 (MQ=255)

  tCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGc                                                                                                                                                >  1:1499327/1‑150 (MQ=255)

   cGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGcc                                                                                                                                               <  2:1211246/150‑1 (MQ=255)

   cGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGGCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGcc                                                                                                                                               <  2:56910/150‑1 (MQ=255)

    ggTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGccc                                                                                                                                              <  2:812587/150‑1 (MQ=255)

    ggTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGGGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGATCCTCTGccc                                                                                                                                              <  2:976325/150‑1 (MQ=255)

     gTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATTAAATGGACGTGTTGATCAATAAAACCTGGGCAGATCCTCTGCCCa                                                                                                                                             <  2:177447/150‑1 (MQ=255)

      tGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCAc                                                                                                                                            <  1:651124/150‑1 (MQ=255)

       gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGATCCTCTGCCCagg                                                                                                                                           <  1:85056/150‑3 (MQ=255)

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

       gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACt                                                                                                                                           >  1:20908/1‑150 (MQ=255)

        cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTg                                                                                                                                          <  2:151329/150‑1 (MQ=255)

        cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTg                                                                                                                                          <  2:148903/150‑1 (MQ=255)

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

         ccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTga                                                                                                                                         <  1:1140252/150‑1 (MQ=255)

          cagcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTgag                                                                                                                                        <  2:267936/150‑1 (MQ=255)

            gcagCCCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGAt                                                                                                                                      >  2:1523735/1‑150 (MQ=255)

            gcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCACACTGAGAt                                                                                                                                      <  2:785436/150‑1 (MQ=255)

               gACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGa                                                                                                                                   >  2:579789/1‑150 (MQ=255)

               gACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGa                                                                                                                                   >  1:588441/1‑150 (MQ=255)

                accaaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGAc                                                                                                                                  <  2:556269/150‑1 (MQ=255)

                accaaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGAc                                                                                                                                  <  2:1250151/150‑1 (MQ=255)

                  caaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACaa                                                                                                                                >  1:732864/1‑150 (MQ=255)

                  caaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACaa                                                                                                                                >  1:129071/1‑150 (MQ=255)

                     cgaCTGAGGTGTCGCCCGGTTCCGTCAGGCGACTTAGCGCCAGTTCCGGGGTGCGGGTCGTGGGACCTGCTTCGCCGCCACGGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGTTGTGCCCACTTAGTTCGACTACCg                                                                                                                             <  2:446200/148‑1 (MQ=255)

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

                        cTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGc                                                                                                                          <  1:1217738/150‑1 (MQ=255)

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

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

                               gacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCgg                                                                                                                   <  1:874450/150‑1 (MQ=255)

                               gacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCgg                                                                                                                   >  2:1365053/1‑150 (MQ=255)

                                acgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGt                                                                                                                  <  2:1197119/150‑1 (MQ=255)

                                  gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCAc                                                                                                                                            >  1:1417525/1‑122 (MQ=255)

                                  gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCAc                                                                                                                                            <  2:1417525/122‑1 (MQ=255)

                                   cccGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACa                                                                                                               <  1:421499/150‑1 (MQ=255)

                                   cccGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACa                                                                                                               >  1:60279/1‑150 (MQ=255)

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

                                          ccGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAg                                                                                                        <  2:710196/150‑1 (MQ=255)

                                          ccGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAg                                                                                                        <  2:838369/150‑1 (MQ=255)

                                          ccGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAg                                                                                                        <  2:1030870/150‑1 (MQ=255)

                                           cGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGa                                                                                                       >  1:1068502/1‑150 (MQ=255)

                                              cAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAgg                                                                                                    <  1:852715/150‑1 (MQ=255)

                                               aGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAggg                                                                                                   >  1:654136/1‑150 (MQ=255)

                                                ggCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGCTCGACAACCGTCCAGTTCGGTACAATGTCAGAAGGGa                                                                                                  >  2:1432608/1‑150 (MQ=255)

                                                ggCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGa                                                                                                  <  2:806349/150‑1 (MQ=255)

                                                   gACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACCACCGTGCAGTTCGGTACAATGTCACAAGGGatat                                                                                               >  1:816240/1‑150 (MQ=255)

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

                                                     cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATTAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                             >  1:385618/1‑150 (MQ=255)

                                                           gcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCa                                                                                       >  1:1138295/1‑150 (MQ=255)

                                                             cACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAAc                                                                                     >  2:117317/1‑150 (MQ=255)

                                                               cTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACgg                                                                                   >  2:881853/1‑150 (MQ=255)

                                                               cTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACgg                                                                                   >  1:778468/1‑150 (MQ=255)

                                                                ttCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAATGGATATTGCTGGCAACGGc                                                                                  <  2:1515531/150‑1 (MQ=255)

                                                                ttCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGc                                                                                  <  2:1394530/150‑1 (MQ=255)

                                                                     gCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATaa                                                                             <  2:1068502/150‑1 (MQ=255)

                                                                      cGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAAAAGGGATATTGCTGGCAACGGCGATAAt                                                                            >  1:486351/1‑150 (MQ=255)

                                                                        tGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAAttt                                                                          >  1:1433602/1‑150 (MQ=255)

                                                                          cgcgTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTg                                                                        <  2:1098027/150‑1 (MQ=255)

                                                                           gcgTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGTAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGc                                                                       <  1:1143062/150‑1 (MQ=255)

                                                                           gcgTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGTAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGc                                                                       <  1:1148380/150‑1 (MQ=255)

                                                                           gcgTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGc                                                                       <  1:579789/150‑1 (MQ=255)

                                                                           gcgTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGc                                                                       <  2:778468/150‑1 (MQ=255)

                                                                            cgtcgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGcc                                                                      <  1:823012/150‑1 (MQ=255)

                                                                               cgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGtt                                                                   <  1:355745/150‑1 (MQ=255)

                                                                                gtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCAAAATTATCGCCGttg                                                                  <  1:438694/150‑2 (MQ=255)

                                                                                  gggACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGc                                                                >  2:801263/1‑150 (MQ=255)

                                                                                       cTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgacga                                                           >  2:370430/1‑150 (MQ=255)

                                                                                        tGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAg                                                          <  2:336484/150‑1 (MQ=255)

                                                                                        tGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAg                                                          >  1:76042/1‑150 (MQ=255)

                                                                                        tGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAg                                                          <  2:943869/150‑1 (MQ=255)

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

                                                                                           ttCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGAc                                                       <  2:588441/150‑1 (MQ=255)

                                                                                            tCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAgg                                                         <  2:142033/147‑1 (MQ=255)

                                                                                                    accgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAAt                                              <  2:129071/150‑1 (MQ=255)

                                                                                                    accgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAAt                                              <  2:732864/150‑1 (MQ=255)

                                                                                                    accgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAAt                                              <  2:563780/150‑1 (MQ=255)

                                                                                                      cgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTc                                            >  2:410208/1‑150 (MQ=255)

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

                                                                                                          aaTCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCg                                        <  2:76042/150‑1 (MQ=255)

                                                                                                           aTCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGa                                       <  1:370430/150‑1 (MQ=255)

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

                                                                                                                 tGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTcc                                 >  1:2252/1‑150 (MQ=255)

                                                                                                                    aCGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGgc                              <  2:332393/150‑1 (MQ=255)

                                                                                                                      gtgtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGgcgc                            <  2:793991/150‑1 (MQ=255)

                                                                                                                       tgtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGgcgta                           >  1:1425166/1‑148 (MQ=255)

                                                                                                                       tgtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCa                           >  1:106421/1‑150 (MQ=255)

                                                                                                                        gtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCAt                          >  2:823250/1‑150 (MQ=255)

                                                                                                                        gtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCAt                          >  2:1033529/1‑150 (MQ=255)

                                                                                                                           gATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACa                       >  2:166146/1‑150 (MQ=255)

                                                                                                                            aTCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTATAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACaa                      <  2:837642/150‑1 (MQ=255)

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

                                                                                                                                       gtGGGCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCt           >  2:229490/2‑150 (MQ=255)

                                                                                                                                       gtGGGCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCt           >  1:924148/2‑150 (MQ=255)

                                                                                                                                         gggCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACTGCGATAATTTTGCCGTTATCGACGAGGACATAGCAAATTGCCCGATCTTCCGGCGCATACAAATGCGCTCCGTgc         <  2:1138295/150‑1 (MQ=255)

                                                                                                                                         gggCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgc         <  2:2252/150‑1 (MQ=255)

                                                                                                                                           gCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCg                                                   <  1:202015/106‑1 (MQ=255)

                                                                                                                                           gCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCg                                                   >  2:202015/1‑106 (MQ=255)

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

                                                                                                                                                ggATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAggg  >  1:828123/1‑150 (MQ=255)

AGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGG  >  NC_000913/4559205‑4559498

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