BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	130,289	A→G	5.1%	G324G (GGT→GGC)	yacH ←	DUF3300 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	130,289	0	A	G	5.1%	633.9 / 14.1	193	G324G (GGT→GGC)	yacH	DUF3300 family protein
Reads supporting (aligned to +/- strand): ref base A (86/97); new base G (3/7); total (89/104)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

Only 100 of 196 total aligned reads displayed.

CATGATAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGTT  >  NC_000913/130141‑130437

           acacCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtg                                                                                                                                           <  1:96308/149‑1 (MQ=255)

            cacCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtgg                                                                                                                                          <  1:1569441/149‑1 (MQ=255)

             acCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtggt                                                                                                                                         >  3:1103116/1‑149 (MQ=255)

             acCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtggt                                                                                                                                         <  1:1337926/149‑1 (MQ=255)

             acCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCgtggt                                                                                                                                         <  3:1136668/149‑1 (MQ=255)

             acCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGGAAGATGCGCACCGGTGATACGGTTGAAAATTTTGACGTCTATTTTTATGGTGTCGCGGTTGTGTTGCCACCCACAACCGACGCGAGGGCTGCCACCgtggt                                                                                                                                         <  4:292240/149‑1 (MQ=255)

             acCATTACGGTAGTTTGGATTGTGCCACCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGGCGGTCATCgtgac                                                                                                                                         <  4:552939/149‑3 (MQ=255)

               cATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtga                                                                                                                                       <  2:893379/149‑1 (MQ=255)

                aTTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtgat                                                                                                                                      <  1:1650506/149‑1 (MQ=255)

                 ttACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtgatg                                                                                                                                     >  3:1591362/1‑149 (MQ=255)

                  tACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtgatga                                                                                                                                    <  4:748984/149‑1 (MQ=255)

                     ggTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATaa                                                                                                                                 <  1:579964/149‑1 (MQ=255)

                        aGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAatca                                                                                                                              <  2:1463998/149‑1 (MQ=255)

                          tttGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTTACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGGGGTGATGATAatcatc                                                                                                                            >  4:797865/1‑149 (MQ=255)

                          tttGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAatcatc                                                                                                                            <  1:1089717/149‑1 (MQ=255)

                          tttGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAatcatc                                                                                                                            >  2:819452/1‑149 (MQ=255)

                               aTTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGt                                                                                                                       >  3:88273/1‑149 (MQ=255)

                               aTTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGt                                                                                                                       <  3:1269903/149‑1 (MQ=255)

                                 tgtgCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCa                                                                                                                     <  2:597214/149‑1 (MQ=255)

                                    gccgccATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatga                                                                                                                  <  4:277434/149‑1 (MQ=255)

                                     ccgccATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat                                                                                                                 >  2:1261517/1‑149 (MQ=255)

                                          atgccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatga                                                                                                            >  3:195311/1‑149 (MQ=255)

                                          atgccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatga                                                                                                            <  1:1330389/149‑1 (MQ=255)

                                          atgccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatga                                                                                                            <  1:1542667/149‑1 (MQ=255)

                                           tgccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatgat                                                                                                           <  2:874170/149‑1 (MQ=255)

                                            gccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatgatg                                                                                                          >  1:875827/1‑149 (MQ=255)

                                            gccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATGACAatgatgattatc                                                                                                          <  4:39575/149‑2 (MQ=255)

                                             ccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatgatg                                                                                                          <  3:121936/148‑1 (MQ=255)

                                             ccatATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatgatg                                                                                                          >  4:121936/1‑148 (MQ=255)

                                               atatTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTc                                                                                                       >  1:1245323/1‑149 (MQ=255)

                                               atatTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTc                                                                                                       >  2:1473271/1‑149 (MQ=255)

                                                   tCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGa                                                                                                   <  4:181435/149‑1 (MQ=255)

                                                     ttATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtc                                                                                                 >  2:725171/1‑149 (MQ=255)

                                                     ttATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtc                                                                                                 >  4:907882/1‑149 (MQ=255)

                                                       aTCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATc                                                                                                                                              <  4:1253850/102‑1 (MQ=255)

                                                        tCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATTTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtc                                                                                              <  1:1317674/149‑1 (MQ=255)

                                                        tCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtc                                                                                              <  4:247792/149‑1 (MQ=255)

                                                        tCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtc                                                                                              <  4:500433/149‑1 (MQ=255)

                                                           gTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtcgtc                                                                                           >  2:189606/1‑149 (MQ=255)

                                                           gTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtcgtc                                                                                           >  4:1624208/1‑149 (MQ=255)

                                                           gTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtcgtc                                                                                           <  2:1365357/149‑1 (MQ=255)

                                                            tAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAtcgtc                                                                                              <  2:1611834/145‑1 (MQ=255)

                                                              aGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCgtcgtcgtc                                                                                        <  4:1260118/149‑1 (MQ=255)

                                                              aGATGCTCACCGGTGATACGGTTGAAACTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCgtcgtcgtc                                                                                        <  3:1161400/149‑1 (MQ=255)

                                                               gATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTcc                                                                                       >  3:498978/1‑149 (MQ=255)

                                                                    tCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTc                                                                                  <  2:1421418/149‑1 (MQ=255)

                                                                      aCCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGAt                                                                                                  <  4:856039/131‑1 (MQ=255)

                                                                      aCCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGa                                                                                <  2:1545683/149‑1 (MQ=255)

                                                                             aTACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATgctgct                                                                         <  3:423460/149‑1 (MQ=255)

                                                                                 ggTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAAt                                                                     <  2:393341/149‑1 (MQ=255)

                                                                                  gTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATa                                                                    >  2:565499/1‑149 (MQ=255)

                                                                                  gTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATa                                                                    >  4:1111671/1‑149 (MQ=255)

                                                                                  gTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATa                                                                    >  2:97792/1‑149 (MQ=255)

                                                                                  gTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATa                                                                    >  1:892805/1‑149 (MQ=255)

                                                                                   ttGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCAtt                                                                                                                         <  4:1144936/95‑1 (MQ=255)

                                                                                   ttGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAg                                                                   >  3:333688/1‑149 (MQ=255)

                                                                                       aaTTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGa                                                                       >  1:850782/1‑141 (MQ=255)

                                                                                        aTTGTTGACGTAGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGtgatga                                                                                                                                    <  4:116554/79‑1 (MQ=255)

                                                                                         ttgttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTa                                                             >  4:529642/1‑149 (MQ=255)

                                                                                            ttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgatgatgatg                                                                                                          <  2:798753/101‑1 (MQ=255)

                                                                                            ttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGAtg                                                                              >  2:128228/1‑129 (MQ=255)

                                                                                            ttgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGAtg                                                                              <  1:128228/129‑1 (MQ=255)

                                                                                             tgACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTg                                                         >  2:1103279/1‑149 (MQ=255)

                                                                                               aCGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat                                                                                                                 <  2:133362/91‑1 (MQ=255)

                                                                                                cGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTa                                                      <  1:565499/149‑1 (MQ=255)

                                                                                                cGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTa                                                      <  2:253937/149‑1 (MQ=255)

                                                                                                 gTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCatgat                                                                                                                 <  3:1254186/89‑1 (MQ=255)

                                                                                                 gTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAg                                                     >  2:673435/1‑149 (MQ=255)

                                                                                                  tCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGc                                                    >  3:268472/1‑149 (MQ=255)

                                                                                                     atatTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAAc                                                 >  1:1093922/1‑149 (MQ=255)

                                                                                                          gATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATTGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCa                                            >  2:507660/1‑149 (MQ=255)

                                                                                                          gATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCAATCGCGGCGATGATAATCATCATTGTCATGATGATGATAGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCa                                            >  4:1360049/1‑149 (MQ=255)

                                                                                                           aTGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCAt                                           <  4:826399/149‑1 (MQ=255)

                                                                                                           aTGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCAt                                           <  1:189606/149‑1 (MQ=255)

                                                                                                             gttgtCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGc                                         <  3:1624208/149‑1 (MQ=255)

                                                                                                                 tCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCAt                                           <  3:704057/143‑1 (MQ=255)

                                                                                                                 tCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCAt                                           >  4:704057/1‑143 (MQ=255)

                                                                                                                 tCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCtata                                     <  1:507660/149‑1 (MQ=255)

                                                                                                                 tCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCtata                                     >  4:1020361/1‑149 (MQ=255)

                                                                                                                    ccGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGcc                                  <  3:1020361/149‑1 (MQ=255)

                                                                                                                       ttGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTc                                                                                  <  1:1423501/98‑1 (MQ=255)

                                                                                                                          tgtTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTg                                                         <  1:1132578/120‑1 (MQ=255)

                                                                                                                          tgtTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTg                                                         >  2:1132578/1‑120 (MQ=255)

                                                                                                                            ttGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCcgc                          <  4:1103116/149‑1 (MQ=255)

                                                                                                                             tgatgattatcatcacCACGATGGCGGTCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCcgcg                         <  2:843739‑M1/133‑1 (MQ=255)

                                                                                                                             tGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCcgcg                         <  3:907882/149‑1 (MQ=255)

                                                                                                                              gccagccaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTa                                                      >  3:1510943/1‑119 (MQ=255)

                                                                                                                              gccagccaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTa                                                      <  4:1510943/119‑1 (MQ=255)

                                                                                                                               ccagccaTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTa                       >  1:1724190/1‑149 (MQ=255)

                                                                                                                                     aTTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGc                 >  3:466556/1‑149 (MQ=255)

                                                                                                                                        aCCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGt              <  1:819452/149‑1 (MQ=255)

                                                                                                                                         ccGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTc             >  2:1173375/1‑149 (MQ=255)

                                                                                                                                         ccGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTc             >  2:137822/1‑149 (MQ=255)

                                                                                                                                           gTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTcaa           <  3:853471/149‑1 (MQ=255)

                                                                                                                                            cCACGATGGCGGTCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGa                                                                                                   <  4:1428592‑M1/59‑1 (MQ=255)

                                                                                                                                            cCACGATGGCGGTCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGa                                                                                                   >  3:1428592‑M1/2‑60 (MQ=255)

                                                                                                                                                GATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGTTGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTTGCAACGCCCATGCTATCGCCGAATCCGCGTACAAAGCtgtcaacacac      <  4:432558‑M1/149‑4 (MQ=255)

                                                                                                                                                GATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATgctgct                                                                         >  1:756365‑M1/1‑82 (MQ=255)

                                                                                                                                                GATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAAc      <  1:1680567‑M1/149‑1 (MQ=255)

                                                                                                                                                    aCCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGtt  >  1:1160672/1‑149 (MQ=255)

CATGATAGGGCACACCATTACGGTAGTTTGGATTGTGCCGCCATGCCATATTCTTATCAGTAAGATGCTCACCGGTGATACGGTTGAAATTGTTGACGTCGATATTGATGTTGTCGCCGTTGTGTTGCCAGCCATTACCGTCACGATGACCGCCATCGTGGTGATGATAATCATCATTGTCATGATGATGATGGTCATGATCGTCGTCGTCCCAGTCGATGCTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGTT  >  NC_000913/130141‑130437

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