BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,648,129	C→T	7.7%	noncoding (67/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,129	0	C	T	7.7%	422.5 / 15.9	142	noncoding (67/84 nt)	agrA	inactive antisense sRNA
Reads supporting (aligned to +/- strand): ref base C (58/72); new base T (4/7); total (62/80)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

Only 100 of 145 total aligned reads displayed.

AACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCT  >  NC_000913/3647980‑3648278

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

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

 aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAATAACCGACGCa                                                                                                                                                       <  1:1224651/150‑1 (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)

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

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

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

           cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ttt                                                                                                                                            <  1:832535/150‑1 (MQ=255)

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

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

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

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

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

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

                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)

                    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)

                         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)

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

                                 aaaCCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACaa                                                                                                                      <  1:370692/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                                                                                                         <  1:609141/150‑1 (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)

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

                                                  cAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGATGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAtt                                                                                                     >  2:1725951/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)

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

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

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

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

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

                                                                                                               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)

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

                                                                                                                           cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCt                            >  2:407220/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)

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

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

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

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

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

                                                                                                                                     aaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCa                  >  2:1485029/1‑150 (MQ=255)

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

                                                                                                                                         aaaaaaCCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGt              <  2:1004183/150‑1 (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   <  2:1290126/150‑1 (MQ=255)

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

AACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCT  >  NC_000913/3647980‑3648278

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