BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,648,128	G→T	7.2%	noncoding (66/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,128	0	G	T	7.2%	375.0 / 12.1	142	noncoding (66/84 nt)	agrA	inactive antisense sRNA
Reads supporting (aligned to +/- strand): ref base G (59/72); new base T (4/7); total (63/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 = 9.86e-01

Only 100 of 147 total aligned reads displayed.

TAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGC  >  NC_000913/3647979‑3648277

tAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg                                                                                                                                                        >  1:1668979/1‑150 (MQ=255)

tAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg                                                                                                                                                        >  2:1501897/1‑150 (MQ=255)

 aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGc                                                                                                                                                       >  1:239507/1‑150 (MQ=255)

 aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGc                                                                                                                                                       >  2:1514042/1‑150 (MQ=255)

  aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAATAACCGACGCa                                                                                                                                                      <  1:1224651/150‑1 (MQ=255)

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

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

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

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

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

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

            cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTATAATTGATGCGTTAGATGATGCCTATCTCACTCCAGTCAGACCCACAAACTCAGGGCTGGAAAGTAAAAAACCTACGCAAAGTCGG‑ttt                                                                                                                                           >  2:290938/1‑150 (MQ=38)

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

            cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ttt                                                                                                                                           >  2:1031956/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                                                                                                                                       >  1:1556300/1‑150 (MQ=255)

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

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

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

                 aGTACATGCCCGTTATAAAAAGAGCACGTCCTTGCAATAGTTTCAGAAGGGTATTAGCATTGATCCGTTAGATCATCTCTCTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGTTGCGTCGG‑ttttttac                                                                                                                                      <  1:49118/150‑3 (MQ=18)

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

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

                     cgtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATc                                                                                                                                  >  1:1698611/3‑150 (MQ=255)

                     cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATc                                                                                                                                  >  1:133184/1‑150 (MQ=255)

                         cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAGACGCAAAGTCGG‑TTTTTTTACTTCCTGa                                                                                                                              >  2:1724003/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)

                                 aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACa                                                                                                                      <  1:204415/150‑1 (MQ=255)

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

                                   aaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAg                                                                                                                    <  1:865794/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                                                                                                        <  1:609141/150‑1 (MQ=255)

                                               cTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                                                                        <  2:901579/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)

                                                     aTAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                                                                  >  2:485389/1‑150 (MQ=255)

                                                       tgTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                                                                <  2:366938/149‑1 (MQ=255)

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

                                                               tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAGTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAAAATGACGATGACCCAGTGAAAg                                                                                        <  1:132093/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)

                                                                         gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                                                              >  2:1433461/1‑150 (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)

                                                                              gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGCCGCACAGTCGG‑GTTTTGTACAGCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAAGCGTCag                                                                         >  1:1534288/1‑149 (MQ=37)

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

                                                                                       gatgatGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                                                                >  1:546334/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)

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

                                                                                                      aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc                                                 <  1:704822/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)

                                                                                                            gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGCCCCAGTGAAAGTATATAAATCGTCACTGAGATATATACCGAAGTGCTCCCTCCGcc                                           >  1:643797/1‑150 (MQ=255)

                                                                                                            gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGcc                                           <  1:1433461/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)

                                                                                                                           tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGc                            >  1:1451130/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                 cTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAAtta                      >  2:423238/1‑149 (MQ=255)

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

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

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

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

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

                                                                                                                                      aaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCa                 >  2:1485029/1‑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‑TTTTTTTAAATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTTCGATATATACCGAAGTGCTAACTCAGACAGCTGAAGAAATAGCTAATTATTCAAATGtt            >  1:1630686/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)

                                                                                                                                                     gCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGc  >  1:575132/1‑150 (MQ=255)

TAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGC  >  NC_000913/3647979‑3648277

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: `‑`