BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,993,305	T→A	5.4%	L718Q (CTG→CAG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,305	0	T	A	5.4%	657.3 / 16.4	200	L718Q (CTG→CAG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base T (99/90); new base A (2/9); total (101/99)
Fisher's exact test for biased strand distribution p-value = 3.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 204 total aligned reads displayed.

CACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGG  >  NC_000913/3993156‑3993454

 aCCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTg                                                                                                                                                      >  2:77182/1‑150 (MQ=255)

    cGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCa                                                                                                                                                   >  1:1060267/1‑150 (MQ=255)

     gCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCTTCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAg                                                                                                                                                  >  1:468229/1‑150 (MQ=255)

       ctgcaacccGGGCTTTAAAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTg                                                                                                                                                <  1:1424200/143‑1 (MQ=255)

       cAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTg                                                                                                                                                >  2:454056/1‑150 (MQ=255)

           aaaCCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAgg                                                                                                                                            <  2:1243107/150‑1 (MQ=255)

            aaCCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACTGACAGGCCGTGCAgtt                                                                                                                                           <  1:215026/150‑3 (MQ=255)

            aaCCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGt                                                                                                                                           <  2:1292224/150‑1 (MQ=255)

               cGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGa                                                                                                                                        <  2:21413/150‑1 (MQ=255)

               cGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGa                                                                                                                                        <  1:390195/150‑1 (MQ=255)

                 ggCGTTTCAAGGCGCTGCGCGGTTCTGGTCAAACCTGGGGGTTGTTCTTCGTACGGCTGTATGTTTCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAACCTGCACTGACAGGCCGTGCAGTTTgtta                                                                                                                                      <  1:795991/150‑4 (MQ=255)

                  gCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAAc                                                                                                                                     <  1:1103930/150‑1 (MQ=255)

                   cGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACTGACAGGCCGTGCAGTTTgttatg                                                                                                                                    <  2:903166/150‑6 (MQ=255)

                     tttCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCaa                                                                                                                                  <  2:164618/150‑1 (MQ=255)

                      ttCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAAt                                                                                                                                 >  1:161604/1‑150 (MQ=255)

                       tCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAACCTGCACTGACAGGCCGTGCAGTTTgt                                                                                                                                        >  1:132219/1‑141 (MQ=255)

                       tCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAACCTGCACTGACAGGCCGTGCAGTTTgt                                                                                                                                        <  2:132219/142‑2 (MQ=255)

                       tCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATc                                                                                                                                <  2:799570/150‑1 (MQ=255)

                       tCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATc                                                                                                                                <  2:161604/150‑1 (MQ=255)

                        cAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCa                                                                                                                               >  1:632940/1‑150 (MQ=255)

                             cacTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCa                                                                                                                          >  2:469458/3‑150 (MQ=255)

                               cTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCaaa                                                                                                                        <  1:181194/150‑1 (MQ=255)

                                  cgcgTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTa                                                                                                                     >  1:66011/1‑150 (MQ=255)

                                   gcgTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTAc                                                                                                                    <  2:666020/150‑1 (MQ=255)

                                   gcgTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTAc                                                                                                                    >  2:147760/1‑150 (MQ=255)

                                    cgTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGaa                                                                                                                                       <  2:322253/130‑1 (MQ=255)

                                    cgTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGaa                                                                                                                                       >  1:322253/1‑130 (MQ=255)

                                    cgTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTAcc                                                                                                                   <  1:533135/150‑1 (MQ=255)

                                         cTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACcggcgg                                                                                                              >  2:659027/1‑150 (MQ=255)

                                           ggTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtg                                                                                                            >  1:247581/1‑150 (MQ=255)

                                           ggTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtg                                                                                                            >  1:225731/1‑150 (MQ=255)

                                            gTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtgg                                                                                                           >  1:846617/1‑150 (MQ=255)

                                            gTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtgg                                                                                                           >  1:786834/1‑150 (MQ=255)

                                            gTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtgg                                                                                                           >  2:416191/1‑150 (MQ=255)

                                             tCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtggt                                                                                                          <  2:432808/150‑1 (MQ=255)

                                             tCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtggt                                                                                                          >  2:1213023/1‑150 (MQ=255)

                                               aaaCCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGtggtgg                                                                                                        <  1:103466/150‑1 (MQ=255)

                                                  ccTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACg                                                                                                     >  1:1233334/1‑150 (MQ=255)

                                                  ccTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACg                                                                                                     >  1:329373/1‑150 (MQ=255)

                                                    tGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGc                                                                                                   <  1:1212874/150‑1 (MQ=255)

                                                     gggggTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGca                                                                                                  >  2:230912/1‑149 (MQ=255)

                                                         gttgttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGc                                                                                              <  2:744877/150‑1 (MQ=255)

                                                          ttgttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGcc                                                                                             <  1:454056/150‑1 (MQ=255)

                                                            gttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAg                                                                                           >  2:1072825/1‑150 (MQ=255)

                                                             ttcttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGc                                                                                          >  2:971096/1‑150 (MQ=255)

                                                              tcttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAAcc                                                                                                                                    >  1:1339776/1‑107 (MQ=255)

                                                              tcttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAAcc                                                                                                                                    <  2:1339776/107‑1 (MQ=255)

                                                              tcttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCg                                                                                         >  1:1348415/1‑150 (MQ=255)

                                                              tcttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCg                                                                                         <  1:1299803/150‑1 (MQ=255)

                                                              tcttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCg                                                                                         <  2:556640/150‑1 (MQ=255)

                                                                ttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGag                                                                                       >  1:614932/1‑149 (MQ=255)

                                                                 tcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAg                                                                                      >  2:834441/1‑150 (MQ=255)

                                                                   gAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAggg                                                                                    >  1:1002956/1‑150 (MQ=255)

                                                                   gAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAggg                                                                                    >  2:948690/1‑150 (MQ=255)

                                                                    aaCGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGa                                                                                   >  1:118178/1‑150 (MQ=255)

                                                                     aCGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGat                                                                                  <  1:983540/150‑1 (MQ=255)

                                                                     aCGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGat                                                                                  <  1:1154316/150‑1 (MQ=255)

                                                                        ccTGAATGTATTGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatcat                                                                               >  2:1182481/1‑150 (MQ=255)

                                                                        ccTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatcat                                                                               >  1:273676/1‑150 (MQ=255)

                                                                        ccTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatcat                                                                               >  2:1252127/1‑150 (MQ=255)

                                                                        ccTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatcat                                                                               >  1:409173/1‑150 (MQ=255)

                                                                           gAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCa                                                                            >  2:1104306/1‑150 (MQ=255)

                                                                           gAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCa                                                                            >  1:244780/1‑150 (MQ=255)

                                                                                  tCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                     <  2:852400/150‑1 (MQ=255)

                                                                                    ggTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttttt                                                                   >  1:986690/1‑150 (MQ=255)

                                                                                    ggTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttttt                                                                   >  1:592310/1‑150 (MQ=255)

                                                                                    ggTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttttt                                                                   <  1:652902/150‑1 (MQ=255)

                                                                                        cAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaa                                                               <  2:390279/150‑1 (MQ=255)

                                                                                         aGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaag                                                              <  2:921955/150‑1 (MQ=255)

                                                                                         aGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaag                                                              <  2:786834/150‑1 (MQ=255)

                                                                                          gAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaaga                                                             >  1:1401209/1‑150 (MQ=255)

                                                                                          gAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaaga                                                             >  1:198667/1‑150 (MQ=255)

                                                                                            aaCTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGaaa                                                           <  2:1196323/150‑1 (MQ=255)

                                                                                              cTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACg                                                         >  1:706971/1‑150 (MQ=255)

                                                                                                ggAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCa                                                       <  1:929280/150‑1 (MQ=255)

                                                                                                       gCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCtt                                                                      >  2:742854/1‑128 (MQ=255)

                                                                                                       gCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACgag                                                >  2:238712/1‑150 (MQ=255)

                                                                                                         cATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatca                                                                                <  1:1207919/116‑1 (MQ=255)

                                                                                                             gAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGc                                          >  1:1244707/1‑150 (MQ=255)

                                                                                                             gAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGc                                          <  2:624434/150‑1 (MQ=255)

                                                                                                                    aTGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAata                                   <  1:77182/150‑1 (MQ=255)

                                                                                                                     tGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAatat                                  <  1:416191/150‑1 (MQ=255)

                                                                                                                       gcgcGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTACCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCt                                >  1:1012993/1‑150 (MQ=255)

                                                                                                                        cgcgATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTa                               >  1:550687/1‑150 (MQ=255)

                                                                                                                        cgcgATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTa                               <  1:233899/150‑1 (MQ=255)

                                                                                                                              ttCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATtct                         >  1:81396/1‑150 (MQ=255)

                                                                                                                               tCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATtctc                        <  2:329373/150‑1 (MQ=255)

                                                                                                                                cGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcg                       >  2:615816/1‑150 (MQ=255)

                                                                                                                                   gtttcAACCTGCACTGACAGGCCGTGCAGTTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgacg                    >  1:598773/5‑150 (MQ=255)

                                                                                                                                   aTAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgacg                    >  2:1426998/1‑150 (MQ=255)

                                                                                                                                    tAACAAACTGCACGGCCTGTCAGTGAAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCGAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgagga                   >  2:574854/1‑150 (MQ=255)

                                                                                                                                      aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaaa                 >  1:722301/1‑150 (MQ=255)

                                                                                                                                       cAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAg                <  1:948690/150‑1 (MQ=255)

                                                                                                                                           cTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGNNNGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAAc            >  2:981809/1‑150 (MQ=255)

                                                                                                                                           cTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAAc            >  2:1142725/1‑150 (MQ=255)

                                                                                                                                            tGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAAcc           >  1:304254/1‑150 (MQ=255)

                                                                                                                                             gCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCg          >  2:992118/1‑150 (MQ=255)

                                                                                                                                               aCGGCCTGTCATTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCggg        <  1:495183/150‑1 (MQ=255)

                                                                                                                                                     tGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAGGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAgg  <  2:614932/150‑1 (MQ=255)

                                                                                                                                                     tGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAgg  <  1:112504/150‑1 (MQ=255)

CACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGG  >  NC_000913/3993156‑3993454

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