BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	361,974	G→C	6.2%	G402G (GGC→GGG)	lacY ←	lactose permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	361,974	0	G	C	6.2%	550.5 / 20.1	162	G402G (GGC→GGG)	lacY	lactose permease
Reads supporting (aligned to +/- strand): ref base G (78/74); new base C (3/7); total (81/81)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 173 total aligned reads displayed.

GCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGCAGAA  >  NC_000913/361825‑362120

gCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccgggg                                                                                                                                                    <  1:918865/150‑1 (MQ=255)

gCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAA‑GCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccggggc                                                                                                                                                   <  1:1516197/150‑1 (MQ=255)

    gCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGc                                                                                                                                                >  2:1343463/1‑150 (MQ=255)

        ttATTCTTTCGGACATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGCCCCGGCCCGCtttc                                                                                                                                            <  1:610607/150‑4 (MQ=255)

         tATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGc                                                                                                                                           >  2:261412/1‑150 (MQ=255)

         tATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGc                                                                                                                                           >  2:2412138/1‑150 (MQ=255)

          aTTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCg                                                                                                                                          <  2:477808/150‑1 (MQ=255)

          aTTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCGCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGCCCCGGCCCGCtttccc                                                                                                                                          <  1:529070/150‑6 (MQ=255)

           ttCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGt                                                                                                                                         >  1:183087/1‑150 (MQ=255)

           ttCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGCCCCGGCCCGCtttccct                                                                                                                                         <  2:455636/150‑7 (MQ=255)

            tCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTg                                                                                                                                        <  2:2442663/150‑1 (MQ=255)

            tCTTTCCGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTg                                                                                                                                        >  2:2371067/1‑150 (MQ=255)

               ttCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAac                                                                                                                                     >  1:1666462/1‑150 (MQ=255)

                  ggTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTGATTCACCGGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAGCACg                                                                                                                                  >  2:1366577/1‑150 (MQ=255)

                       ttGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAAt                                                                                                                             <  2:1559195/150‑1 (MQ=255)

                          gCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTaa                                                                                                                          >  1:790952/1‑150 (MQ=255)

                          gCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTaa                                                                                                                          <  1:1281993/150‑1 (MQ=255)

                           cATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAg                                                                                                                         >  2:68502/1‑150 (MQ=255)

                            aagaacAATACGATCACACCGTTATGAAATGTTGGTCTGATAAGGCGCTAGCGCCGCATACCACATTGATTGATTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAgg                                                                                                                        >  2:1102856/6‑150 (MQ=255)

                            aTGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAgg                                                                                                                        <  2:1006844/150‑1 (MQ=255)

                            aTGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAgg                                                                                                                        <  1:175409/150‑1 (MQ=255)

                              gTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTg                                                                                                                      <  1:192542/150‑1 (MQ=255)

                                 cAATGCGATCACTCCGTTATTATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAg                                                                                                                   <  2:94118/150‑1 (MQ=255)

                                 cAATGCGATCACTCCGTTATGATATGTTGGTCGGATAATGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTTAACACGGAAATTAAGGTGAAg                                                                                                                   <  2:2114496/150‑1 (MQ=255)

                                  aaTGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGc                                                                                                                  <  2:1516301/150‑1 (MQ=255)

                                      cGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAg                                                                                                              <  1:2427717/150‑1 (MQ=255)

                                      cGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAg                                                                                                              >  2:1308933/1‑150 (MQ=255)

                                      cGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAg                                                                                                              >  1:480233/1‑150 (MQ=255)

                                       gATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCGTCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCCAAGCGTGAACACGGAAATTAAGGTGAAGCCCAgc                                                                                                             >  1:2462655/1‑150 (MQ=255)

                                        aTCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAgcg                                                                                                            <  2:2326730/150‑1 (MQ=255)

                                            cTCCGTTATGATATGTTGGTCGGATACGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCAGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccac                                                                                                        >  1:733957/1‑150 (MQ=255)

                                             tCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccacc                                                                                                       <  2:2116644/150‑1 (MQ=255)

                                              ccGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAATGTGAAGCCCAGCGccacca                                                                                                      <  2:1967351/150‑1 (MQ=255)

                                               cGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGcc                                                                                                          >  1:2387314/1‑145 (MQ=255)

                                               cGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAg                                                                                                     >  2:2430456/1‑150 (MQ=255)

                                               cGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAg                                                                                                     <  2:365288/150‑1 (MQ=255)

                                                    tGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCa                                                                                                <  1:1210034/150‑1 (MQ=255)

                                                        atcaatgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAc                                                                                            >  1:1130432/7‑150 (MQ=38)

                                                         tGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAcc                                                                                           <  2:145509/150‑1 (MQ=255)

                                                         tGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAcc                                                                                           >  1:1563043/1‑150 (MQ=255)

                                                         tGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAcc                                                                                           <  2:1725140/150‑1 (MQ=255)

                                                          gTTGGTCTGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCa                                                                                          <  1:2093998/150‑1 (MQ=255)

                                                          gTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCa                                                                                          >  1:2021446/1‑150 (MQ=255)

                                                          gTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCa                                                                                          <  2:723316/150‑1 (MQ=255)

                                                            atgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGa                                                                                        >  1:1221517/3‑150 (MQ=255)

                                                            atgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCAACTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACCGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGa                                                                                        >  2:59630/3‑150 (MQ=255)

                                                               tCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAATCGACTTCATTCACCTGACGACGCAGAATGGAAAGCGGTCCGCGGCCGCTAAGCATGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATaa                                                                                     >  2:1058265/1‑150 (MQ=255)

                                                                cGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAg                                                                                    <  2:80447/150‑1 (MQ=255)

                                                                cGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAg                                                                                    <  1:1311484/150‑1 (MQ=255)

                                                                  gATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACACAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAgcg                                                                                  >  2:877581/1‑150 (MQ=255)

                                                                       ggCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTg                                                                             >  1:2444983/1‑150 (MQ=255)

                                                                         cgcTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGa                                                                           >  1:1707068/1‑150 (MQ=255)

                                                                             cgcgCCGCATCCGACGTTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAAcc                                                                       <  2:2237822/150‑1 (MQ=255)

                                                                                gccgcATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGAt                                                                    <  2:2192249/150‑1 (MQ=255)

                                                                                   gcATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCt                                                                 >  2:1229771/1‑150 (MQ=255)

                                                                                     aTCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGTGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCttt                                                               <  1:1666267/150‑1 (MQ=255)

                                                                                      tCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTc                                                              <  1:126765/150‑1 (MQ=255)

                                                                                       ccGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTca                                                             <  2:1130432/150‑1 (MQ=255)

                                                                                        cGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcat                                                            <  1:1343463/150‑1 (MQ=255)

                                                                                         gACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcata                                                           <  1:1660217/150‑1 (MQ=255)

                                                                                           cATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcataca                                                         >  2:2097810/1‑150 (MQ=255)

                                                                                           cATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcataca                                                         <  1:120118/150‑1 (MQ=255)

                                                                                            attgattgCTTAAGCGCCTTCATTCACCTGACGCCGCAGCAGGGAAAGCGCGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGCACCAGCACCAGATAACCGCCCTGGAACCCGAGGCTTTcatacat                                                        >  2:1576088/1‑150 (MQ=255)

                                                                                             ttgattgCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcatacata                                                       <  2:1291455/150‑1 (MQ=255)

                                                                                                  tgCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGcc                                                  <  1:2268667/150‑1 (MQ=255)

                                                                                                   gCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcatacata                                                       <  2:1671567/144‑1 (MQ=255)

                                                                                                   gCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGccc                                                 >  2:1087556/1‑150 (MQ=255)

                                                                                                   gCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGccc                                                 >  2:723272/1‑150 (MQ=255)

                                                                                                    cTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCg                                                >  1:1684241/1‑150 (MQ=255)

                                                                                                       aaGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCa                                             >  2:792760/1‑150 (MQ=255)

                                                                                                        aGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAATTTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAg                                            >  1:952239/1‑150 (MQ=255)

                                                                                                         gCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGc                                                                  >  1:71179/1‑127 (MQ=255)

                                                                                                         gCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGc                                                                  >  1:80270/1‑127 (MQ=255)

                                                                                                         gCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGc                                                                  <  2:80270/127‑1 (MQ=255)

                                                                                                           gACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAATTTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATTCTTTAATAGATATTGCTTGCCAGTAc                                         <  2:1919255/150‑1 (MQ=255)

                                                                                                            aCTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACa                                        >  1:698499/1‑150 (MQ=255)

                                                                                                                 attcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACAt                                   <  2:774293/150‑1 (MQ=255)

                                                                                                                  ttcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATa                                  >  2:2075816/1‑150 (MQ=255)

                                                                                                                   tcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATaa                                 >  1:764327/1‑150 (MQ=255)

                                                                                                                    caCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGACCACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATaaa                                >  1:270898/1‑150 (MQ=255)

                                                                                                                     aCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATaaaa                               >  2:1392493/1‑150 (MQ=255)

                                                                                                                        tGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAatc                            <  2:1855493/150‑1 (MQ=255)

                                                                                                                         gacgacGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAatca                           <  1:340685/150‑1 (MQ=255)

                                                                                                                          acgacgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAatcat                          <  1:1002119/150‑1 (MQ=255)

                                                                                                                           cgacgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAatcatc                         >  2:1672607/1‑150 (MQ=255)

                                                                                                                             acgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACAGGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGGCAGTACAGACATAAAAATCATCGc                       <  2:1661150/150‑1 (MQ=255)

                                                                                                                             acgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGc                       <  2:155762/150‑1 (MQ=255)

                                                                                                                              cgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCATCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGcc                      <  2:120446/150‑1 (MQ=255)

                                                                                                                                cagcagGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAg                    <  2:1922119/150‑1 (MQ=255)

                                                                                                                                    agGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGc                >  2:359231/1‑150 (MQ=255)

                                                                                                                                     gggAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCt               >  2:552442/1‑150 (MQ=255)

                                                                                                                                     gggAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCt               <  1:768558/150‑1 (MQ=255)

                                                                                                                                     gggAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCt               >  1:1425498/1‑150 (MQ=255)

                                                                                                                                     cgcttaGCGGCCCCGGCCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGACACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAACTCATCGCCAGTTGCt               >  1:359590/6‑150 (MQ=255)

                                                                                                                                     cgcttaGCGGCCCCGGCCCGCTAAGCGTGAACACGGAAATTAAG‑‑GAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTa             >  1:729062/6‑150 (MQ=255)

                                                                                                                                      ggAAAGCGGGCCGGGGCCGATAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCCCCAGACCCAGCACCAGATAAGCGCCCTGGACACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGATGCtt              >  1:1560826/1‑150 (MQ=255)

                                                                                                                                             gggccggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAaagaag       >  1:1954320/1‑150 (MQ=255)

                                                                                                                                             gggccggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAaagaag       >  1:635205/1‑150 (MQ=255)

                                                                                                                                                  ggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGCAGaa  <  1:274028/150‑1 (MQ=255)

                                                                                                                                                  ggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGCAGaa  <  1:46935/150‑1 (MQ=255)

GCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGCAGAA  >  NC_000913/361825‑362120

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