BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,648,132	A→C	7.6%	noncoding (70/84 nt)	agrA →	inactive antisense sRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,648,132	0	A	C	7.6%	411.7 / 18.8	141	noncoding (70/84 nt)	agrA	inactive antisense sRNA
Reads supporting (aligned to +/- strand): ref base A (58/72); new base C (4/7); total (62/79)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 146 total aligned reads displayed.

TGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAAT  >  NC_000913/3647983‑3648281

tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaaa                                                                                                                                                        <  2:205483/150‑1 (MQ=255)

 ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAg                                                                                                                                                       >  2:704822/1‑150 (MQ=255)

    tAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCg                                                                                                                                                    <  2:627226/150‑1 (MQ=255)

    tAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGACACCAACTCAGGGCTGGAAAGTAAAAAACCTACGCAAAGTCg                                                                                                                                                    >  1:548150/1‑150 (MQ=255)

     aaTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCgg                                                                                                                                                   <  2:338709/150‑1 (MQ=255)

     aaTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCgg                                                                                                                                                   <  2:110136/150‑1 (MQ=255)

       tCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑tt                                                                                                                                                >  2:1443675/1‑150 (MQ=255)

        cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ttt                                                                                                                                               >  2:1031956/1‑150 (MQ=255)

        cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ttt                                                                                                                                               >  1:1601717/1‑150 (MQ=255)

         aTACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑tttt                                                                                                                                              <  2:1620906/150‑1 (MQ=255)

          aaCAGTACATGCAGGTTATAAACCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGG‑ttttt                                                                                                                                             <  1:835912/149‑1 (MQ=38)

           aCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGG‑tttttt                                                                                                                                            <  2:1137586/150‑1 (MQ=38)

            cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ttttttt                                                                                                                                           <  2:507601/150‑1 (MQ=255)

            cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ttttttt                                                                                                                                           >  1:1556300/1‑150 (MQ=255)

             aggacgtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTa                                                                                                                                          >  2:1228406/7‑150 (MQ=25)

             aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGG‑ttttttac                                                                                                                                          <  1:965724/150‑3 (MQ=38)

                acgtGCTGGTTTTATAACCTGCATGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACAt                                                                                                                                       >  1:1004183/4‑150 (MQ=25)

                    gCATGTTATAAAACCAGCACGTACTAGCAATAGTTTCAGTATGTTATTAGCATTGATGCGTTAGATTATGTCTATCTCACTCCAGTCAGATCCACCAACTCAGGGCTGGTAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCgg                                                                                                                                   <  1:491645/150‑1 (MQ=25)

                     cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACTTCCGGa                                                                                                                                  >  1:168310/1‑150 (MQ=255)

                     cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAGACGCAAAGTCGG‑TTTTTTTACTTCCTGa                                                                                                                                  >  2:1724003/1‑150 (MQ=255)

                      aGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGAt                                                                                                                                 >  2:1442688/1‑150 (MQ=255)

                        gTTATATAACCAGCACGTCGTTGCAATAGTTTCAGTTTGGTATTAGCATTGATGCGTTAGTTGATGGCTATCTCAATCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCTACGCAAAGTCGG‑TTTTTTTACATCCGGATTc                                                                                                                               <  2:829347/150‑1 (MQ=25)

                         ttATAAAACAAGCACGTCCTTTCAATAGTTTCATTATGGTATTAGCATTGATGAGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCg                                                                                                                              <  2:1102626/150‑1 (MQ=255)

                          tataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGTCTGGAAAGTAAAAAACCGACGCAAAGGCGG‑TTTTTTTACATCCGGATTCgg                                                                                                                             <  1:378462/150‑1 (MQ=255)

                           ataAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGa                                                                                                                            <  2:239507/150‑1 (MQ=255)

                               aaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACaac                                                                                                                        >  1:769993/1‑149 (MQ=255)

                               aaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAg                                                                                                                        <  1:865794/150‑1 (MQ=255)

                                aCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAgg                                                                                                                       <  1:1451113/150‑1 (MQ=255)

                                aCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAgg                                                                                                                       >  1:435664/1‑150 (MQ=255)

                                aCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAAGCAAAGACGG‑TTATTTTACATCCTGATTCGGACAAgg                                                                                                                       >  1:100494/1‑150 (MQ=255)

                                 ccAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGc                                                                                                                      <  2:133184/150‑1 (MQ=255)

                                   aGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCtt                                                                                                                    <  2:949718/150‑1 (MQ=255)

                                       cGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                                                                <  1:399095/150‑1 (MQ=255)

                                        gTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACGTCCAGATTCGTACATGGCTtattat                                                                                                               >  1:475610/1‑150 (MQ=25)

                                           cTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                                                                            <  2:901579/150‑1 (MQ=255)

                                           cTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                                                                            <  1:609141/150‑1 (MQ=255)

                                            ttGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                                                                           <  1:560295/150‑1 (MQ=255)

                                               cAATAGTTTCAGTATGTTTTTAGGATTGATGCGTTTGATGATGGCTATCTCACTCCAGTGAGAGCGAGCAACTGAGGGCTGGAAAGTTAATACCCGACGCAAAGTGGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                                                                        <  1:1707760/150‑1 (MQ=255)

                                               cAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGATGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAtt                                                                                                        >  2:1725951/1‑149 (MQ=255)

                                               cAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGATGGAAATTAAAAAACCGACGCAAAGACGG‑TTTTTTTACATCCGGATACGGACAAGCCTTAATATGAAGAta                                                                                                        >  2:1575566/1‑149 (MQ=255)

                                                aaTTGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                                                                       >  1:327410/1‑150 (MQ=255)

                                                aaTAGTTTCAGTTTGGGATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                                                                       <  1:311200/150‑1 (MQ=255)

                                                           tctGGTAGTAGCATTGATGCGTTATATGATGGCTATCTCACTCCAGGCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                                                                            <  2:1429535/148‑1 (MQ=255)

                                                           tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAGTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAAAATGACGATGACCCAGTGAAAg                                                                                            <  1:132093/150‑1 (MQ=255)

                                                             tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                                                                          <  1:668048/150‑1 (MQ=255)

                                                                tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCTGACAAGTCTTAATATGACGATGACCCAGTGAAAGtatat                                                                                       >  2:1210068/1‑150 (MQ=255)

                                                                 tttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCATCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                                                                      <  1:94019/149‑1 (MQ=255)

                                                                     gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                                                                  <  1:370525/150‑1 (MQ=255)

                                                                      cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCTGATACGTACAAGCCATATTATGACGAAGAGCCAGTGAAAGTATATAAATCg                                                                                 >  1:799741/1‑150 (MQ=21)

                                                                      cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                                                                 >  1:638782/1‑150 (MQ=255)

                                                                         tGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                                                                              <  2:655691/150‑1 (MQ=255)

                                                                          gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAGATCGTCAc                                                                             <  1:627678/150‑1 (MQ=255)

                                                                              cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                                                                         >  2:473920/1‑150 (MQ=255)

                                                                                   gatgatGGCTATCTCACTCCAGTCAGAGCCACCAAATCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCTG‑TTTTTTTACATACTGATTCGGACAAGGCTTAATATGAAGATGACCCAGTGAAAGTATATAAATCGTCAATCCGatata                                                                    >  2:1038786/1‑150 (MQ=25)

                                                                                      tatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc                                                                 <  1:1193181/149‑1 (MQ=255)

                                                                                       atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                                                                <  1:1172595/150‑1 (MQ=255)

                                                                                        tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACACAGTGAAAGTATATAAATAGACACTCACATATATAcag                                                               >  1:1326157/1‑148 (MQ=25)

                                                                                           cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg                                                            >  1:1661911/1‑150 (MQ=255)

                                                                                                  aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc                                                     <  1:704822/150‑1 (MQ=255)

                                                                                                  aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc                                                     <  2:1491964/150‑1 (MQ=255)

                                                                                                  aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc                                                     >  2:1155676/1‑150 (MQ=255)

                                                                                                  aCTCCAGTCAGAGCCACCAACTCAGGGATGGAAATTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc                                                     <  2:787067/150‑1 (MQ=255)

                                                                                                     ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc                                                  <  2:1596605/150‑1 (MQ=255)

                                                                                                         tCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCa                                              <  2:1661911/150‑1 (MQ=255)

                                                                                                            gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGcg                                           >  1:383406/1‑149 (MQ=255)

                                                                                                             agCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTg                                          <  1:1031956/150‑1 (MQ=255)

                                                                                                              gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAGGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTga                                         <  2:768371/150‑1 (MQ=255)

                                                                                                              gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTga                                         <  1:1700547/150‑1 (MQ=255)

                                                                                                                caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaag                                       <  1:1228406/150‑1 (MQ=255)

                                                                                                                  ccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTgaagaa                                     >  2:1005129/1‑150 (MQ=255)

                                                                                                                   caACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGATTAAAATTACGAAGAACCAGAGAAATTATATAAATCGTAAATGCGATATATACAGAAGTGCTCCCTAAGCCAGATAAAGaaa                                    >  2:895982/1‑150 (MQ=21)

                                                                                                                    aaCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGaaa                                    >  2:912135/1‑150 (MQ=255)

                                                                                                                     aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAAt                                   <  1:912135/150‑1 (MQ=255)

                                                                                                                     aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATc                                  <  1:79016/150‑1 (MQ=255)

                                                                                                                        cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCt                               <  1:1472988/150‑1 (MQ=255)

                                                                                                                        cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCt                               >  2:407220/1‑150 (MQ=255)

                                                                                                                          gggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGATGTGCTCCCTCCGCCAGCTGAAGAAATCTCTat                             >  2:85223/1‑149 (MQ=255)

                                                                                                                          gggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCaaa                             >  1:1282110/1‑150 (MQ=255)

                                                                                                                          gggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTaa                             >  2:56388/1‑150 (MQ=255)

                                                                                                                               ggAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCtt                        >  2:602636/1‑150 (MQ=255)

                                                                                                                               ggAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCtt                        >  2:542371/1‑150 (MQ=255)

                                                                                                                                gAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTg                       <  1:485389/150‑1 (MQ=255)

                                                                                                                                 aaaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                      <  1:1347599/150‑1 (MQ=255)

                                                                                                                                 aaaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                      <  1:1378144/150‑1 (MQ=255)

                                                                                                                                 aaaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                      <  1:1443675/150‑1 (MQ=255)

                                                                                                                                  aaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCa                     <  2:638782/150‑1 (MQ=255)

                                                                                                                                   gtaaaaaaaCCGACTTTGCGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACAGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCaa                    >  2:477191/4‑150 (MQ=255)

                                                                                                                                    taaaaaaaCCGACTTTGCGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAAt                   >  2:1686944/3‑150 (MQ=255)

                                                                                                                                     tAAAAAACCGACGCAAAGTCGG‑TTTTATTACATCCGGATTCGGACAAGGCTTAATATGACTATGACCCAGTGAAAGTATATACATCGTCACTGCTATATATACCGAATTGATCCCTACGCCAGCTGAAGAAATCGCTAATTATTGCAAtt                  >  2:510449/1‑149 (MQ=255)

                                                                                                                                     aaaaaaaCCGACTTTGCGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATg                  >  2:1507766/2‑150 (MQ=255)

                                                                                                                                      aaaaaaCCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGt                 <  2:1004183/150‑1 (MQ=255)

                                                                                                                                       aaaaaCCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGtt                >  1:1479846/1‑150 (MQ=255)

                                                                                                                                       aaaaaCCGACGCAAAGTCGG‑TTTTTTTAAATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTTCGATATATACCGAAGTGCTAACTCAGACAGCTGAAGAAATAGCTAATTATTCAAATGtt                >  1:1630686/1‑150 (MQ=255)

                                                                                                                                       aaaaaCCGACGCAAAGTCGG‑TTTATATACATACCGAGTAGGAGAAATATTAATATTACGATTAAACATTTAATGTATATAAATAGTAAATAATAAAAATACCGAAGTGCGCGCTCCGCCAGCTGAAGCAGTCGGTAATGGTGGCAATGtt                >  2:348758/1‑150 (MQ=255)

                                                                                                                                         aaaCCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCCTCACTGCGATATATACCGAAGTGCTAAATCCGCCAGCTGAAGAAAACGCTATTTCATGCAATGTTAg              >  1:1538568/1‑150 (MQ=255)

                                                                                                                                          aaCCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGc             <  2:400256/150‑1 (MQ=255)

                                                                                                                                                 gCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGc      <  2:1290126/150‑1 (MQ=255)

                                                                                                                                                  cAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCt     <  1:1442688/150‑1 (MQ=255)

                                                                                                                                                     aGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAAt  <  1:673694/150‑1 (MQ=255)

                                                                                                                                                     aGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAAt  >  1:525846/1‑150 (MQ=255)

TGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAAT  >  NC_000913/3647983‑3648281

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