BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,742,795	A→T	5.1%	E138V (GAA→GTA) ‡	dgcN →	diguanylate cyclase DgcN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,742,795	0	A	T	5.1%	488.2 / 9.5	150	E138V (GAA→GTA) ‡	dgcN	diguanylate cyclase DgcN
Reads supporting (aligned to +/- strand): ref base A (75/67); new base T (3/5); total (78/72)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 153 total aligned reads displayed.

TTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAG  >  NC_000913/2742648‑2742944

tttCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGaaa                                                                                                                                                     <  1:135472/150‑1 (MQ=255)

 ttCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAAc                                                                                                                                                    >  2:2442856/1‑150 (MQ=255)

  tCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAAcc                                                                                                                                                   <  1:18462/150‑1 (MQ=255)

       tGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTTTCACCATTgt                                                                                                                                              <  1:2325623/150‑2 (MQ=255)

        gCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAAcc                                                                                                                                                   <  1:375666/144‑1 (MQ=255)

        gCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGc                                                                                                                                             <  2:863861/150‑1 (MQ=255)

         cAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTTTCACCATTgtga                                                                                                                                            <  1:1989120/150‑4 (MQ=255)

         cAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCg                                                                                                                                            >  1:522784/1‑150 (MQ=255)

         cAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCg                                                                                                                                            <  1:916644/150‑1 (MQ=255)

          aGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGa                                                                                                                                           >  1:1228903/1‑150 (MQ=255)

           gAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGaa                                                                                                                                          >  1:1369340/1‑150 (MQ=255)

                 cGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGc                                                                                                                                    >  1:938741/1‑150 (MQ=255)

                 cGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGc                                                                                                                                    <  1:318471/150‑1 (MQ=255)

                  gttaTAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATGATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCt                                                                                                                                   <  1:2122096/147‑1 (MQ=255)

                   tGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCtt                                                                                                                                  >  2:1628088/1‑150 (MQ=255)

                     aTAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTaa                                                                                                                                >  1:1554857/1‑150 (MQ=255)

                      tAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAAc                                                                                                                               >  1:1632471/1‑150 (MQ=255)

                        agcagcAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCg                                                                                                                             <  2:1369340/150‑1 (MQ=255)

                        agcagcAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCg                                                                                                                             >  2:1072465/1‑150 (MQ=255)

                         gcagcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGc                                                                                                                            >  1:2145446/1‑150 (MQ=255)

                          cagcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCt                                                                                                                           >  2:1795334/1‑150 (MQ=255)

                          cagcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCt                                                                                                                           >  2:2275035/1‑150 (MQ=255)

                            gcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGTTACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTcg                                                                                                                         <  1:2175772/150‑1 (MQ=255)

                                  tatTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAg                                                                                                                   >  2:365540/1‑150 (MQ=255)

                                      cTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCa                                                                                                               <  2:961857/150‑1 (MQ=255)

                                      cTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCa                                                                                                               >  1:428436/1‑150 (MQ=255)

                                         gCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATc                                                                                                            >  2:1151708/1‑150 (MQ=255)

                                             ccTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGcc                                                                                                        <  2:1437544/150‑1 (MQ=255)

                                                 gCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCAttt                                                                                                    <  1:1987336/150‑1 (MQ=255)

                                                       cacCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTg                                                                                              >  2:91227/1‑150 (MQ=255)

                                                         cccGTAAGGATCCATGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCACCATCATTCAGCCGATTCGTCAAAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGt                                                                                            >  2:1769026/1‑150 (MQ=255)

                                                         cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTg                                                                                              >  2:1279750/1‑148 (MQ=255)

                                                         cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTg                                                                                              <  1:1279750/148‑1 (MQ=255)

                                                         cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGt                                                                                            >  2:386661/1‑150 (MQ=255)

                                                         cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGACATTTTATCTGGt                                                                                            >  1:975894/1‑150 (MQ=255)

                                                            gTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGtttt                                                                                         <  1:2222674/150‑1 (MQ=255)

                                                               aGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgc                                                                                      <  1:2076946/150‑1 (MQ=255)

                                                               aGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgc                                                                                      >  1:302688/1‑150 (MQ=255)

                                                                   tCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGTCATTTTATCTGGTTTtcgctcgc                                                                                  <  2:1530427/150‑1 (MQ=255)

                                                                     cAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCg                                                                                <  2:1081052/150‑1 (MQ=255)

                                                                      aGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGt                                                                               <  1:1078605/150‑1 (MQ=255)

                                                                         cGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACt                                                                            >  1:488157/1‑150 (MQ=255)

                                                                          gATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTg                                                                           >  1:1616344/1‑150 (MQ=255)

                                                                               tttCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCgg                                                                      >  1:2479188/1‑150 (MQ=255)

                                                                               tttCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCgg                                                                      >  2:599072/1‑150 (MQ=255)

                                                                                ttCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGt                                                                     <  1:468197/150‑1 (MQ=255)

                                                                                  cAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTg                                                                   <  2:586104/150‑1 (MQ=255)

                                                                                   aGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGt                                                                  >  2:1352151/1‑150 (MQ=255)

                                                                                    gTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTa                                                                 >  1:2010322/1‑150 (MQ=255)

                                                                                    gTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTa                                                                 <  2:1134341/150‑1 (MQ=255)

                                                                                    gTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTa                                                                 >  1:237877/1‑150 (MQ=255)

                                                                                    gTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTa                                                                 <  1:1285892/150‑1 (MQ=255)

                                                                                      aaTTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTAtt                                                               >  2:721654/1‑150 (MQ=255)

                                                                                       aTTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCGGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTAttn                                                              <  1:121807/150‑2 (MQ=255)

                                                                                          tCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctgc                                                           >  1:889878/1‑150 (MQ=255)

                                                                                          tCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctgc                                                           >  2:2300589/1‑150 (MQ=255)

                                                                                           cATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctgct                                                          >  1:1004978/1‑150 (MQ=255)

                                                                                            aTAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctgctg                                                         >  1:2321762/1‑150 (MQ=255)

                                                                                            aTAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctgctg                                                         >  1:2352274/1‑150 (MQ=255)

                                                                                             tAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTgg                                                        <  1:2334821/150‑1 (MQ=255)

                                                                                             tAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCAGATTCGTCACAATGGTGAACCCATTGGCGAAGTAAGCTTAACCGCTCGCGACAGTTCAATAAGCCATTTTATATGGTTTTCGCTCGCCGTACTTACCGGTTGTATTCTGCTgg                                                        >  2:794306/1‑150 (MQ=255)

                                                                                               aGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCa                                                      <  1:1812808/150‑1 (MQ=255)

                                                                                                 ccACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATc                                                    <  1:72383/150‑1 (MQ=255)

                                                                                                   aCTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAg                                                  <  2:1216553/150‑1 (MQ=255)

                                                                                                    cTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAgg                                                 <  1:482554/150‑1 (MQ=255)

                                                                                                    cTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAgg                                                 >  2:705044/1‑150 (MQ=255)

                                                                                                      ggCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCa                                               <  1:192270/150‑1 (MQ=255)

                                                                                                        ctctTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATc                                             >  2:1241746/1‑150 (MQ=255)

                                                                                                         tctTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCg                                            <  1:2442856/150‑1 (MQ=255)

                                                                                                          ctTCCCTGCCTCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGc                                           <  2:2227836/150‑1 (MQ=255)

                                                                                                          ctTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATGTGGTTTTCGCTCGCCGTACTGCCCGGTTGTATTCTGCTGGCATCAGGCATCGc                                           >  2:2100539/1‑150 (MQ=255)

                                                                                                          ctTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGc                                           <  2:428436/150‑1 (MQ=255)

                                                                                                            tCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCAt                                              >  2:1508755/1‑145 (MQ=255)

                                                                                                            tCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCAt                                              <  1:1508755/145‑1 (MQ=255)

                                                                                                            tCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCaa                                         >  2:2476666/1‑150 (MQ=255)

                                                                                                             cccTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAAt                                        >  2:1992979/1‑150 (MQ=255)

                                                                                                              ccTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAAtt                                       >  2:1937949/1‑150 (MQ=255)

                                                                                                                  cccccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTAccc                                   >  2:819696/1‑150 (MQ=255)

                                                                                                                  cccccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTAccc                                   >  2:518074/1‑150 (MQ=255)

                                                                                                                  cccccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTAccc                                   >  2:687503/1‑150 (MQ=255)

                                                                                                                  cccccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTAccc                                   <  2:889878/150‑1 (MQ=255)

                                                                                                                    cccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTc                                 >  2:2263862/1‑150 (MQ=255)

                                                                                                                    cccATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTc                                 >  2:2164028/1‑150 (MQ=255)

                                                                                                                         catTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCAcccc                            >  2:1444423/1‑149 (MQ=255)

                                                                                                                             cAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCAt                        >  2:900809/1‑150 (MQ=255)

                                                                                                                                   agtCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGcac                  <  1:2008315/148‑1 (MQ=255)

                                                                                                                                     tCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACaa                <  1:609565/150‑1 (MQ=255)

                                                                                                                                       gTCACAATGGTGAAACCATTGGCGGAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATg              >  2:2185805/1‑150 (MQ=255)

                                                                                                                                         gccaATGGTTTCACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGc            >  2:377583/3‑150 (MQ=255)

                                                                                                                                         gccaATGGTTTCACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGc            >  1:517327/3‑150 (MQ=255)

                                                                                                                                         gccaATGGTTTCACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGc            >  2:2049877/3‑150 (MQ=255)

                                                                                                                                          acaATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGcc           <  1:1509275/150‑1 (MQ=255)

                                                                                                                                          acaATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGcc           <  2:495570/150‑1 (MQ=255)

                                                                                                                                             aTGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTgg        <  1:365540/150‑1 (MQ=255)

                                                                                                                                               tgTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTa      <  1:1628088/149‑1 (MQ=255)

                                                                                                                                               ggTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTa      >  1:849880/1‑150 (MQ=255)

                                                                                                                                               ggTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTa      <  2:1269797/150‑1 (MQ=255)

                                                                                                                                                 tGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGa    >  1:2264902/1‑150 (MQ=255)

                                                                                                                                                  gAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGaa   >  2:2079210/1‑150 (MQ=255)

                                                                                                                                                   aaaCCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAg  <  2:83157/150‑1 (MQ=255)

TTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAGAAG  >  NC_000913/2742648‑2742944

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