BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,908,659	A→G	F261L (TTT→CTT)	pstA ←	phosphate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,908,659	0	A	G	100.0%	379.0 / NA	104	F261L (TTT→CTT)	pstA	phosphate ABC transporter permease
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (57/47); total (57/47)

Only 100 of 104 total aligned reads displayed.

TAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGT  >  NC_000913/3908519‑3908806

tAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATc                                                                                                                                           >  2:287137/1‑151 (MQ=255)

tAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATc                                                                                                                                           >  1:449010/1‑151 (MQ=255)

 aGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCg                                                                                                                                          >  3:866133/1‑151 (MQ=255)

  gggTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCAt                                                                                                                                            <  3:117598/148‑1 (MQ=255)

  gggTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGc                                                                                                                                         >  3:421035/1‑151 (MQ=255)

    gTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCa                                                                                                                                        >  3:164241/1‑150 (MQ=255)

     agAGCGGAAAACCGTGTCAACAATATCAACCGAGTTTATTCTTCGCAAAACCAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATGCCGCAAGCGGGCTCATCGCaa                                                                                                                                       <  3:1007317/149‑1 (MQ=255)

      gAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATAATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAAc                                                                                                                                     >  3:46607/1‑151 (MQ=255)

        ggggAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACt                                                                                                                                    <  3:891403/150‑1 (MQ=255)

        ggggAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACt                                                                                                                                    <  2:1246115/150‑1 (MQ=255)

         gggAAATCCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTa                                                                                                                                  <  1:147923/151‑1 (MQ=255)

         gggAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTa                                                                                                                                  <  2:1350242/151‑1 (MQ=255)

               aCCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTaa                                                                                                                                 >  2:1101035/1‑146 (MQ=255)

               aCCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATc                                                                                                                            <  2:1323557/151‑1 (MQ=255)

                  gtgtCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTc                                                                                                                         >  2:248318/1‑151 (MQ=255)

                   tgtCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTc                                                                                                                         >  2:1005190/1‑150 (MQ=255)

                    gtCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCAc                                                                                                                       <  2:659977/151‑1 (MQ=255)

                       aGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCgg                                                                                                                    >  3:25008/1‑151 (MQ=255)

                       aGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCgg                                                                                                                    >  2:1296220/1‑151 (MQ=255)

                         cAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGc                                                                                                                   >  3:583711/1‑150 (MQ=255)

                          aaTATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGc                                                                                                                   >  2:1018245/1‑149 (MQ=255)

                          aaTATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCa                                                                                                                  <  3:31142/150‑1 (MQ=255)

                          aaTATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAg                                                                                                                 <  3:67344/151‑1 (MQ=255)

                           atatCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAgg                                                                                                                >  2:1231926/1‑151 (MQ=255)

                                aaCCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGc                                                                                                           <  1:150138/151‑1 (MQ=255)

                                aaCCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGc                                                                                                           <  2:1165035/151‑1 (MQ=255)

                                 aCCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCg                                                                                                          >  2:126385/1‑151 (MQ=255)

                                  ccGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAgg                                                                                                                >  2:628389/1‑144 (MQ=255)

                                   cGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGAt                                                                                                        <  3:1238687/151‑1 (MQ=255)

                                     tgtTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCg                                                                                                      >  2:989651/1‑151 (MQ=255)

                                      gtTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCgg                                                                                                     >  2:79313/1‑151 (MQ=255)

                                      gtTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCgg                                                                                                     <  3:457874/151‑1 (MQ=255)

                                        ttATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGc                                                                                                    >  3:1298427/1‑150 (MQ=255)

                                        ttATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGc                                                                                                    >  3:802432/1‑150 (MQ=255)

                                         tATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCt                                                                                                   >  3:60817/1‑150 (MQ=255)

                                         tATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCt                                                                                                   <  3:1275855/150‑1 (MQ=255)

                                           ttcttcGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGc                                                                                                 >  3:76073/1‑150 (MQ=255)

                                            tcttcGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGca                                                                                                <  2:237481/150‑1 (MQ=255)

                                             cttcGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGcatc                                                                                              >  2:138848/1‑151 (MQ=255)

                                                cGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATg                                                                                           >  3:654291/1‑151 (MQ=255)

                                                cGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATg                                                                                           <  3:598764/151‑1 (MQ=255)

                                                cGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATg                                                                                           >  3:667665/1‑151 (MQ=255)

                                                 gCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGt                                                                                          >  2:1174630/1‑151 (MQ=255)

                                                 gCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGt                                                                                          >  3:744436/1‑151 (MQ=255)

                                                   aaaaaCAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTcc                                                                                        <  1:286961/151‑1 (MQ=255)

                                                   aaaaaCAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTcc                                                                                        >  3:354569/1‑151 (MQ=255)

                                                   aaaaaCAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGGTGGCGATCGGCTGCATCATGTcc                                                                                        >  3:747403/1‑151 (MQ=255)

                                                    aaaaCAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCg                                                                                       >  2:462844/1‑151 (MQ=255)

                                                     aaaCAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGt                                                                                      <  2:242307/151‑1 (MQ=255)

                                                         aacGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTc                                                                                  >  2:891605/1‑151 (MQ=255)

                                                          acGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTcc                                                                                 >  3:180874/1‑151 (MQ=255)

                                                           cgcgcgCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCa                                                                                <  2:686311/151‑1 (MQ=255)

                                                             cgcgcCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTcc                                                                                 >  3:321133/1‑148 (MQ=255)

                                                             cgcgcCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTcc                                                                                 >  3:544299/1‑148 (MQ=255)

                                                               cgcCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAAc                                                                            >  1:260893/1‑151 (MQ=255)

                                                                 ccAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTg                                                                          >  2:76656/1‑151 (MQ=255)

                                                                   aGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTgg                                                                         >  2:766388/1‑150 (MQ=255)

                                                                   aGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGt                                                                        >  1:253824/1‑151 (MQ=255)

                                                                    gAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGtt                                                                       >  3:963188/1‑151 (MQ=255)

                                                                     aaTGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTg                                                                      <  4:241193/151‑1 (MQ=255)

                                                                      aTGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTgg                                                                     <  2:702476/151‑1 (MQ=255)

                                                                          tCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTTGTTGgaga                                                                  <  4:200901/150‑1 (MQ=255)

                                                                          tCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGgag                                                                   >  1:240736/1‑149 (MQ=255)

                                                                           cagcagcagTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGagag                                                                 >  3:315934/1‑150 (MQ=255)

                                                                             gcagcagTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAgcgc                                                              <  2:153391/151‑1 (MQ=255)

                                                                              cagcagTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGcgcg                                                             <  1:183870/151‑1 (MQ=255)

                                                                              cagcagTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGcgcg                                                             >  4:339484/1‑151 (MQ=255)

                                                                                  agTACGCACAGGGTAATGATCAATACCGCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTa                                                          >  3:234401/1‑150 (MQ=255)

                                                                                     aCGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAAc                                                       <  2:292784/150‑1 (MQ=255)

                                                                                      cGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAAc                                                       >  3:348656/1‑149 (MQ=255)

                                                                                        cacaGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAAcagca                                                   <  4:249872/151‑1 (MQ=255)

                                                                                          caGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACagcagc                                                 <  2:164218/151‑1 (MQ=255)

                                                                                              gTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAgcggcg                                             <  2:1397674/151‑1 (MQ=255)

                                                                                              gTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAgcggcg                                             >  2:614918/1‑151 (MQ=255)

                                                                                               tAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAgcggcg                                             >  3:969446/1‑150 (MQ=255)

                                                                                                aaTGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGcgcg                                           >  2:1130768/1‑151 (MQ=255)

                                                                                                aaTGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGcgcg                                           <  2:529982/151‑1 (MQ=255)

                                                                                                      cAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTc                                      <  3:989692/150‑1 (MQ=255)

                                                                                                        aTACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCAcc                                   <  3:844284/151‑1 (MQ=255)

                                                                                                            cccGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCa                               <  2:866095/151‑1 (MQ=255)

                                                                                                             ccGGCCCAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCa                               >  3:617679/1‑150 (MQ=255)

                                                                                                                 cccAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACg                          <  3:628415/151‑1 (MQ=255)

                                                                                                                 cccAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACg                          >  3:401491/1‑151 (MQ=255)

                                                                                                                 cccAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACg                          >  2:459656/1‑151 (MQ=255)

                                                                                                                  ccAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACgg                         <  3:601347/151‑1 (MQ=255)

                                                                                                                  ccAGGCAAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGGTTCACCAGCAATACgg                         >  2:1040161/1‑151 (MQ=255)

                                                                                                                   cAGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCCGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACcgg                        >  3:322706/1‑151 (MQ=255)

                                                                                                                    aGGCCAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGc                       <  3:597549/151‑1 (MQ=255)

                                                                                                                        cAATTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATc                   >  3:398998/1‑151 (MQ=255)

                                                                                                                         aaTTGCTGCCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCaaac                   >  3:970106/1‑148 (MQ=255)

                                                                                                                             gctgcCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCcag              <  3:287162/151‑1 (MQ=255)

                                                                                                                             gctgcCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCcag              <  1:339487/151‑1 (MQ=255)

                                                                                                                             gctgcCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCca               <  3:138865/150‑1 (MQ=255)

                                                                                                                                gcCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCcagcag           <  2:802399/151‑1 (MQ=255)

                                                                                                                                gcCATTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCcagca            <  2:60806/150‑1 (MQ=255)

                                                                                                                                   aTTCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGAt        >  3:950121/1‑151 (MQ=255)

                                                                                                                                     tCCGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATa       <  2:579411/150‑1 (MQ=255)

                                                                                                                                       cGCAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCg    >  2:990059/1‑151 (MQ=255)

                                                                                                                                         cAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGt  <  2:315904/151‑1 (MQ=255)

                                                                                                                                         cAAGCGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGt  <  3:126401/151‑1 (MQ=255)

TAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGT  >  NC_000913/3908519‑3908806

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