BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,991,334	T→A	5.8%	I61N (ATT→AAT)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,334	0	T	A	5.8%	386.7 / 9.2	130	I61N (ATT→AAT)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base T (56/66); new base A (2/6); total (58/72)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 131 total aligned reads displayed.

CTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGT  >  NC_000913/3991186‑3991482

  ggATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTATTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTg                                                                                                                                                   <  1:1965952/150‑1 (MQ=255)

  ggATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTg                                                                                                                                                   <  2:1096551/150‑1 (MQ=255)

  ggATGCCATAAATCAATTGCGTGTCGATCTCGCTCTTCCTCCTATCAGTCCAAAAAAACAACAGGTCTACAGTCTACTGCCGCCATTGTTGCCCTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTg                                                                                                                                                   >  2:179931/1‑150 (MQ=255)

   gATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAGCGTTCCCAAAGGCATTTGc                                                                                                                                                  >  1:1017221/1‑150 (MQ=255)

     tGCCCTAAATCATTTTCTTGTTGATCGCGCGCTTGCTGCTTGTGGGCCTGCATTCCAACAGGTCTACATTCTACTGACTACCTTGTTGCACTATCAACATCCGCTAATGCCTTGTTACCTTGATGGTAACGTTCCCAAAGGGATTTGCCt                                                                                                                                                <  2:1166309/150‑1 (MQ=255)

      gCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCtt                                                                                                                                               <  1:437584/150‑1 (MQ=255)

        cATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTTATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                                                                                                             >  2:153451/1‑150 (MQ=255)

        cATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCAAACAGGTATACAGTCTACTGCCGACATTGTTGAACTATCACCATCCGCTAATGACGGCTTACCTTGATGTTAACGTTCCCAAAGGCATTTGCCTTTa                                                                                                                                             >  1:862308/1‑150 (MQ=255)

          tAAATCATTTGCGTGTGGATCGAGCGCTTGCTGCTATGGGTCCTGCATTCGAACAGGTCTACTGTATACTGCCGACATTGTTGCACTATCTCCATCCGCTAATGCCGGGTTACCTTGATGGTATCGTTCCCAAAGGCAAATGCCtttggg                                                                                                                                           <  1:1667207/150‑4 (MQ=255)

          tAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTaca                                                                                                                                           <  1:776049/150‑1 (MQ=255)

          tAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttggg                                                                                                                                           <  2:367921/150‑4 (MQ=255)

          tAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttggg                                                                                                                                           <  2:1845315/150‑4 (MQ=255)

          aaaaTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttggg                                                                                                                                           <  1:665144/149‑4 (MQ=255)

              tCAATTGCGTGTTGATCGCGCGCTTGCTGCTATGGGTCCTTCATTCCAACATGTATACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAATGCATTTGCCTTTACACGcc                                                                                                                                       <  2:795384/150‑1 (MQ=255)

              tCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGcc                                                                                                                                       <  1:1502045/150‑1 (MQ=255)

               cAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCCCTATCACCATCCGCTTATGCCGGGTTACCGTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACCCGCCt                                                                                                                                      >  1:1475570/1‑150 (MQ=255)

                aaTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtg                                                                                                                                     <  1:1230234/150‑1 (MQ=255)

                 aTTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtga                                                                                                                                    <  1:603654/150‑1 (MQ=255)

                  ttccGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCTTTCCAAGTGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgat                                                                                                                                   <  1:13970/147‑1 (MQ=255)

                   tGCGTGTGGATAGCGCGCTTGCTGCTATGGGGCCTGCATTCAAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgatg                                                                                                                                  <  1:1919508/150‑1 (MQ=255)

                          ggATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCa                                                                                                                           >  2:1572730/1‑150 (MQ=255)

                          ggATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCa                                                                                                                           >  2:1790557/1‑150 (MQ=255)

                          ggATCGCGCGCTTGCTGCTATGGGGCCTGCATGCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGAGGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCa                                                                                                                           <  2:59918/150‑1 (MQ=255)

                           gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCaa                                                                                                                          >  1:641917/1‑150 (MQ=255)

                            aTCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                                                                                                         >  1:550069/1‑150 (MQ=255)

                                      tgctgcTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCt                                                                                                               <  1:1079469/150‑1 (MQ=255)

                                       gctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTg                                                                                                              <  2:686888/150‑1 (MQ=255)

                                        ctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTgc                                                                                                             >  1:1566092/1‑149 (MQ=255)

                                         tgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCCCTACCTGaa                                                                                                            >  1:1021781/1‑150 (MQ=255)

                                           ctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTTCCGACATTGTTGCACTATCACCATCCTCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACg                                                                                                          <  1:158888/150‑1 (MQ=255)

                                            tATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGa                                                                                                         >  1:1095306/1‑150 (MQ=255)

                                               ggggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt                                                                                                      >  2:880615/1‑150 (MQ=255)

                                               ggggCCTGCATTCCAACAGGTCTACAGTATACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTATCATTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTAAACCACACTACCTGAACGAgat                                                                                                      >  1:112357/1‑148 (MQ=255)

                                               ggggCCTCCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt                                                                                                      <  2:451942/150‑1 (MQ=255)

                                                 ggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGACCGAGCTTg                                                                                                    >  2:1581382/1‑150 (MQ=255)

                                                 ggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg                                                                                                    >  1:1397908/1‑150 (MQ=255)

                                                 ggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg                                                                                                    >  2:39743/1‑150 (MQ=255)

                                                 ggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg                                                                                                    >  1:124937/1‑150 (MQ=255)

                                                 ggCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGAAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAAAGAGCtta                                                                                                    >  1:744124/1‑149 (MQ=255)

                                                   ccTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCCCTACCTGAACGAGCTTGaa                                                                                                  >  2:1589097/1‑150 (MQ=255)

                                                   ccTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGaa                                                                                                  >  1:1783305/1‑150 (MQ=255)

                                                      gCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTg                                                                                               <  1:800899/150‑1 (MQ=255)

                                                       cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGt                                                                                              >  1:621927/1‑150 (MQ=255)

                                                            cAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt                                                                                         >  1:757750/1‑150 (MQ=255)

                                                             aaCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCAAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTg                                                                                        <  1:924518/150‑1 (MQ=255)

                                                                    cTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc                                                                                 >  2:608987/1‑150 (MQ=255)

                                                                       cAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGt                                                                              >  2:708510/1‑150 (MQ=255)

                                                                       cAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGGTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGt                                                                              <  1:1291225/150‑1 (MQ=255)

                                                                         gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCAATTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTAc                                                                            >  1:79664/1‑150 (MQ=255)

                                                                                gCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATcc                                                                     >  2:1668465/1‑150 (MQ=255)

                                                                                gCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATcc                                                                     <  1:128222/150‑1 (MQ=255)

                                                                                gCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATcc                                                                     >  2:101462/1‑150 (MQ=255)

                                                                                  cGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTTCAGGATccgc                                                                   >  2:2028831/1‑150 (MQ=255)

                                                                                  cGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccgc                                                                   >  1:508561/1‑150 (MQ=255)

                                                                                  cGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccgc                                                                   <  2:1077409/150‑1 (MQ=255)

                                                                                   gACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCAAAAGGCATTTGCCTTTACCCGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATccgcc                                                                  >  1:939297/1‑150 (MQ=255)

                                                                                       ttgttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGaaa                                                              >  1:1506096/1‑150 (MQ=255)

                                                                                       ttgttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTACCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGaaa                                                              <  2:1560697/150‑1 (MQ=255)

                                                                                         gttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTCCCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAgg                                                            <  1:364383/150‑1 (MQ=255)

                                                                                         gttgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAgg                                                            >  2:83415/1‑150 (MQ=255)

                                                                                           tgCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTg                                                          <  1:1602703/150‑1 (MQ=255)

                                                                                            ccACTATCACCATCCGCTACTGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGa                                                         <  2:1106134/149‑1 (MQ=255)

                                                                                             cACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAg                                                        >  1:1112027/1‑150 (MQ=255)

                                                                                               cTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                      <  1:264457/150‑1 (MQ=255)

                                                                                               cTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCt                                                      <  2:1073415/150‑1 (MQ=255)

                                                                                                 aTCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGCAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTc                                                    >  1:1796469/1‑150 (MQ=255)

                                                                                                  tCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACTCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTcc                                                   >  1:2025443/1‑150 (MQ=255)

                                                                                                  tCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTcc                                                   <  2:1939322/150‑1 (MQ=255)

                                                                                                   cACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTAAGTACAGGATCCGCCGAAAGGTGAGATTCCa                                                  >  2:230742/1‑150 (MQ=255)

                                                                                                     ccATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCGTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAAt                                                <  1:1233496/150‑1 (MQ=255)

                                                                                                          cGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTg                                           >  1:928430/1‑150 (MQ=255)

                                                                                                             tAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGgtg                                        <  1:608987/150‑1 (MQ=255)

                                                                                                             tAAGGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTTAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAATGTGAGCTTCCAATTACTGgtg                                        <  2:1933775/150‑1 (MQ=255)

                                                                                                                tGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGtat                                     >  1:529154/1‑150 (MQ=255)

                                                                                                                    gggTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATacac                                 <  1:501638/150‑1 (MQ=255)

                                                                                                                     ggTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCATTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACAcc                                <  1:1969699/150‑1 (MQ=255)

                                                                                                                     ggTTACCTTGATGGTAACGTTCCCAAAGGAATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACAcc                                <  1:632505/150‑1 (MQ=255)

                                                                                                                         aCCTTGATGGTAACGTTCCCTAAGGCATTTGGCTTTAGACGCCTGCTGCTACGCATCGCCACTACCTGAACGAGCTTGATCTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCGACCAATTACTGGTGTATACACCATgg                            <  1:299459/150‑1 (MQ=255)

                                                                                                                          ccTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATggg                           >  1:1704290/1‑150 (MQ=255)

                                                                                                                          ccTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACAGTATCGTGGAAGGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATAACACATggg                           >  2:536669/1‑150 (MQ=255)

                                                                                                                             tGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGATACTCATCGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgcag                        >  1:103684/1‑150 (MQ=255)

                                                                                                                              gATGTTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCTCTACCTGAACGAGCTTGATCTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACTCCATGGgcagc                       <  2:112357/150‑1 (MQ=255)

                                                                                                                              gATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgcagc                       <  2:1106771/150‑1 (MQ=255)

                                                                                                                              gATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGgcagc                       >  1:1149210/1‑150 (MQ=255)

                                                                                                                                tGGTTACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCATTTACTGGTGTATACACCATGGGCAGCAc                     >  2:249800/1‑150 (MQ=255)

                                                                                                                                tGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAc                     >  1:1830237/1‑150 (MQ=255)

                                                                                                                                   tAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAACCTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTc                  >  2:1038811/1‑150 (MQ=255)

                                                                                                                                    aaCGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCg                 <  2:124937/150‑1 (MQ=255)

                                                                                                                                      cGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGCGCTTCCAATTACTGGTGTATACACCATGGTCAGCACCTCGTc               >  2:553955/1‑150 (MQ=255)

                                                                                                                                       gtccccAACGGCATTTGCCCTTACACGCCCGATGAAACTCACCCCCACAACCTGAACGCGCGTGAACTGTCTCGTGCAATGCCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTcc              <  1:1361365/147‑1 (MQ=255)

                                                                                                                                        ttCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCg             <  1:708511/150‑1 (MQ=255)

                                                                                                                                        ttCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCg             >  2:1740839/1‑150 (MQ=255)

                                                                                                                                         gtgtaaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAcc                    <  1:327362/138‑1 (MQ=255)

                                                                                                                                         gtgtaaaGGCAAATGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCAcc                    >  2:327362/5‑142 (MQ=255)

                                                                                                                                          cccAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTa           >  2:388695/1‑150 (MQ=255)

                                                                                                                                           ccAAAGGCCTTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCCTGGGCAGCACCTCGTCCGTAg          >  2:624709/1‑150 (MQ=255)

                                                                                                                                             aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAggg        <  1:1687713/150‑1 (MQ=255)

                                                                                                                                                  cATTTGCCTTTACACGCCTGATGAAACTCCACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATAAACCATGGGAAGCCCCTCGTCAGTAGTGCAAAg   >  1:1285946/1‑150 (MQ=255)

                                                                                                                                                  cATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAg   >  2:1987108/1‑150 (MQ=255)

                                                                                                                                                   aTTTGCCTTTACACGCGTGATGAAACTCAACGCCACTACCTGAACGAGCTTGATCTGTATCGTGGAATGTCTGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAgt  <  2:723608/150‑1 (MQ=255)

CTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGT  >  NC_000913/3991186‑3991482

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