BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,648,131	A→G	7.8%	noncoding (69/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,131	0	A	G	7.8%	394.8 / 19.0	141	noncoding (69/84 nt)	agrA	inactive antisense sRNA
Reads supporting (aligned to +/- strand): ref base A (57/71); new base G (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 = 8.02e-01

Only 100 of 143 total aligned reads displayed.

TGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTA  >  NC_000913/3647983‑3648279

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

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

   aTAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTc                                                                                                                                                   <  2:1725243/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                                                                                                                                             <  1:832535/150‑1 (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                                                                                                                                         <  2:654881/150‑1 (MQ=255)

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

            cAGAACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATAAGCAATGATGCGTTAGATGATGCCTATCTCACTCCAGGCAGAGCCACCAACTCAGCGCTGGAAAGTAAAAAACCGACTTTGCGTCGG‑tttttta                                                                                                                                         <  1:166128/150‑2 (MQ=255)

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

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

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

             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)

                  aTGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCTTTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGG‑ttttttactttcc                                                                                                                                   <  2:1096411/150‑8 (MQ=37)

                    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)

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

                         ttATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCg                                                                                                                            <  2:1009021/150‑1 (MQ=255)

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

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

                             aaaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACAAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTAAATCCGGATTCGGACa                                                                                                                        >  1:874510/1‑150 (MQ=255)

                              aaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACaa                                                                                                                       <  1:370692/150‑1 (MQ=255)

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

                                aCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAgg                                                                                                                     <  1:1451113/150‑1 (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)

                                         tCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATg                                                                                                            <  1:394163/150‑1 (MQ=255)

                                          ccTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                                                                           <  1:787592/150‑1 (MQ=255)

                                          ccTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATACGGACAAGGCTTAATATGa                                                                                                           >  2:1130042/1‑150 (MQ=255)

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

                                           cTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                                                                          >  1:1572403/1‑150 (MQ=255)

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

                                             tGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                                                                        >  1:963374/1‑150 (MQ=255)

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

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

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

                                                  tAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                                                                   <  1:463780/150‑1 (MQ=255)

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

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

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

                                                               gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                                                                      <  2:581652/150‑1 (MQ=255)

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

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

                                                                     gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                                                                >  2:1433461/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                   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)

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

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

                                                                                                                               ggAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCtt                      >  2:602636/1‑150 (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)

                                                                                                                                     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)

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

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

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

TGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTA  >  NC_000913/3647983‑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: `‑`