BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,648,130	A→T	7.7%	noncoding (68/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,130	0	A	T	7.7%	403.5 / 15.9	141	noncoding (68/84 nt)	agrA	inactive antisense sRNA
Reads supporting (aligned to +/- strand): ref base A (58/72); new base T (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 144 total aligned reads displayed.

ACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTA  >  NC_000913/3647981‑3648279

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

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

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

   ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAg                                                                                                                                                     <  2:1559799/150‑1 (MQ=255)

     aTAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTc                                                                                                                                                   <  2:1725243/150‑1 (MQ=255)

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

       aaTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCgg                                                                                                                                                 <  2:110136/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)

            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)

               aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATTGTATTAGCATTGATGCGTTAGATGTTGGCTTTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTa                                                                                                                                        <  2:1286547/150‑1 (MQ=255)

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

               aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGG‑ttttttac                                                                                                                                        <  2:542259/150‑3 (MQ=38)

                  acgtGCTGGTTTTATAACCTGCATGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACAt                                                                                                                                     >  1:1004183/4‑150 (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)

                           ttttataaCCTGCATGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCg                                                                                                                            <  2:1683740/144‑1 (MQ=255)

                           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: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)

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

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

                                                  aaTAGTTTCAGTTTGGGATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                                                                     <  1:311200/150‑1 (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)

                                                                 gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                                                                      <  2:581652/150‑1 (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)

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

                                                                                     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)

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

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

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

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

                                                                                                    aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc                                                   <  2:1491964/150‑1 (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)

                                                                                                           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                             >  2:407220/1‑150 (MQ=255)

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

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

                                                                                                                               cTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAAtta                        >  2:423238/1‑149 (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:1443675/150‑1 (MQ=255)

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

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

                                                                                                                                   aaaGTAAAAAACCGAAGAAAAGTCGG‑TTTTTTTACATCCGGATTCGTACAAGGCCTAATAGGACGATGACCCAGTGAAAGTATATAAATCGACACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGc                    <  2:120165/150‑1 (MQ=255)

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

                                                                                                                                     gtaaaaaaaCCGACTTTGCGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACAGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCaa                  >  2:477191/4‑150 (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)

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

                                                                                                                                          aaaaCCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTa             <  2:1601717/150‑1 (MQ=255)

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

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

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

                                                                                                                                                     aaaGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTa  <  1:697751/150‑1 (MQ=255)

ACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTA  >  NC_000913/3647981‑3648279

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