BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,378,859	T→C	7.1%	H29R (CAT→CGT) ‡	ampC ←	beta‑lactamase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,378,859	0	T	C	7.1%	754.9 / 33.8	252	H29R (CAT→CGT) ‡	ampC	beta‑lactamase
Reads supporting (aligned to +/- strand): ref base T (118/116); new base C (10/8); total (128/124)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 254 total aligned reads displayed.

TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTT  >  NC_000913/4378711‑4379008

tAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                                                                                      <  1:1785795/150‑1 (MQ=255)

 aaCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                                                                                     <  2:944368/150‑1 (MQ=255)

 aaCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                                                                                     >  1:320868/1‑150 (MQ=255)

   cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                                                                                   >  2:660666/1‑150 (MQ=255)

     cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACaa                                                                                                                                                 >  1:1132388/1‑150 (MQ=255)

     cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACaa                                                                                                                                                 <  1:445689/150‑1 (MQ=255)

      aaaCAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                                                                                >  1:433393/1‑150 (MQ=255)

        acaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCACCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                                                                              >  1:615397/1‑150 (MQ=255)

         caacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                                                                             >  1:1002154/1‑150 (MQ=255)

          aacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                                                                            <  1:1203874/150‑1 (MQ=255)

           acGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                                                                                           <  1:662673/150‑1 (MQ=255)

           acGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCATCGCACAattaccc                                                                                                                                           <  1:645341/150‑6 (MQ=255)

            cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                                                                          >  1:1973674/1‑150 (MQ=255)

             gTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                                                                         <  1:1619009/150‑1 (MQ=255)

                tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGATCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                                                                                      >  2:850061/1‑150 (MQ=255)

                    gtgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                                                                                  >  1:350800/1‑150 (MQ=255)

                     tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                                                                                 <  2:1633553/150‑1 (MQ=255)

                     tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                                                                                 <  1:1173955/150‑1 (MQ=255)

                     tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                                                                                 >  2:979595/1‑150 (MQ=255)

                     tatgACGGGCTGCTTTTTGGCGATGTCAGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                                                                                 <  2:1277195/148‑1 (MQ=255)

                      gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                                                                                >  1:898624/1‑150 (MQ=255)

                      gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                                                                                >  1:1538380/1‑150 (MQ=255)

                        gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                                                                              >  2:802532/1‑150 (MQ=255)

                          cGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                                                                            <  1:373415/150‑1 (MQ=255)

                          cGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                                                                               >  2:609529/1‑147 (MQ=255)

                               tgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                                                                                       <  2:165899/150‑1 (MQ=255)

                                gctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                                                                         <  1:1370822/131‑1 (MQ=255)

                                gctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTATACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                                                                                      >  2:1415513/1‑150 (MQ=255)

                                       ggCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATTAATTACCGCCACCGCCATACGCGGGAGCTGTTGTTGCGCTATAAGCGGGGTAATTGTGCGATGCACAAGATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                                                                               >  1:1251763/1‑150 (MQ=255)

                                       ggCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                                                                               <  2:370108/150‑1 (MQ=255)

                                       ggCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                                                                               >  1:1422602/1‑150 (MQ=255)

                                         cGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                                                                             >  2:151212/1‑150 (MQ=255)

                                         cGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                                                                             <  2:1345346/150‑1 (MQ=255)

                                          gATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                                                                            >  1:689441/1‑150 (MQ=255)

                                                 gCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                                                                                     >  2:1195330/1‑150 (MQ=255)

                                                  cATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                                                                                    <  2:1132388/150‑1 (MQ=255)

                                                      gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                                                                                      >  1:1328114/1‑144 (MQ=255)

                                                      gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                                                                                <  1:788398/150‑1 (MQ=255)

                                                       ccccAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                                                                               <  1:698595/150‑1 (MQ=255)

                                                        cccAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                                                                              >  2:1867358/1‑150 (MQ=255)

                                                         ccAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                                                                             <  2:320868/150‑1 (MQ=255)

                                                          cAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                                                                            >  1:603774/1‑150 (MQ=255)

                                                                  gTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                                                                                    >  1:475508/1‑150 (MQ=255)

                                                                   taataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                                                                                   >  1:406319/1‑150 (MQ=255)

                                                                    aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                  >  2:1457760/1‑150 (MQ=255)

                                                                    aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                  >  1:577349/1‑150 (MQ=255)

                                                                    aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGAGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                  <  1:1276149/150‑1 (MQ=255)

                                                                     ataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                                                                                 >  2:378712/1‑150 (MQ=255)

                                                                       aaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                                                                                        <  2:1583162/125‑1 (MQ=255)

                                                                        aGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt                                                                              >  1:1677506/1‑150 (MQ=255)

                                                                        aGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt                                                                              <  1:1926982/150‑1 (MQ=255)

                                                                          gTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATAGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                                                                            >  1:1615938/1‑150 (MQ=255)

                                                                             tACCCTGATAAATTACCGCCACCGCCATACACGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGAGATGCACAATATCGTTTATTTGTTGTGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAATGCGCATAGCGTCGttt                                                                         <  2:1930723/150‑1 (MQ=255)

                                                                              aCCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt                                                                        >  2:905776/1‑150 (MQ=255)

                                                                               cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                       >  2:631448/1‑150 (MQ=255)

                                                                               cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                       >  1:830439/1‑150 (MQ=255)

                                                                               accTGATAATTTACCGCCACGGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                       <  2:358811/149‑1 (MQ=255)

                                                                                ccTGATAAATTACCGCCACATCCATACCAGGGATCTTTTGTTGCTATATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTTGAGCAAGAGGAGGTAATTAATAAGGCGCAGAGCGTCGTTTTGa                                                                      <  2:1211480/150‑1 (MQ=255)

                                                                                 cTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa                                                                     >  1:1099392/1‑150 (MQ=255)

                                                                                       aaTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg                                                               <  2:1540950/150‑1 (MQ=255)

                                                                                        aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                                                                                           >  1:150027/1‑73 (MQ=255)

                                                                                        aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                              <  1:1160504/150‑1 (MQ=255)

                                                                                         ttACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGt                                                             >  1:771299/1‑150 (MQ=255)

                                                                                           accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCt                                                           <  2:1091398/150‑1 (MQ=255)

                                                                                               ccaccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAGTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGtt                                                       >  1:1657163/1‑150 (MQ=255)

                                                                                                caccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGttt                                                      <  1:573045/150‑1 (MQ=255)

                                                                                                 accgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTc                                                     <  2:433393/150‑1 (MQ=255)

                                                                                                  ccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTcc                                                    >  1:121878/1‑150 (MQ=255)

                                                                                                   cgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCa                                                   <  1:1021067/150‑1 (MQ=255)

                                                                                                   cgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCa                                                   >  1:1887131/1‑150 (MQ=255)

                                                                                                   cgccaCACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCa                                                   >  2:1296261/1‑150 (MQ=255)

                                                                                                    gccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCAt                                                  >  1:213917/1‑150 (MQ=255)

                                                                                                       aTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACa                                               <  2:238154/150‑1 (MQ=255)

                                                                                                         aCCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACaaa                                             <  1:1451327/150‑1 (MQ=255)

                                                                                                          cctgggATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGGGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACaaaa                                            >  2:1828420/4‑150 (MQ=255)

                                                                                                          cccGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACaaaa                                            >  2:1027723/1‑150 (MQ=255)

                                                                                                                  cTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGc                                    >  2:1372231/1‑150 (MQ=255)

                                                                                                                    tttGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCa                                  >  2:135820/1‑150 (MQ=255)

                                                                                                                      tgttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATa                                <  2:1002154/150‑1 (MQ=255)

                                                                                                                           ctctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCgg                           <  2:559811/150‑1 (MQ=255)

                                                                                                                             ctATAAGCTGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                                                                                    <  2:260929/75‑1 (MQ=255)

                                                                                                                             ctATAAGCTGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                                                                                    >  1:260929/1‑75 (MQ=255)

                                                                                                                             ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGc                         >  2:1763956/1‑150 (MQ=255)

                                                                                                                             ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGc                         >  2:333446/1‑150 (MQ=255)

                                                                                                                               ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCg                       <  1:52479/150‑1 (MQ=255)

                                                                                                                                taAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCgg                      >  1:5565/1‑150 (MQ=255)

                                                                                                                                   gCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                   >  1:1984637/1‑150 (MQ=255)

                                                                                                                                    caacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                  >  1:1405022/8‑150 (MQ=255)

                                                                                                                                    cGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                  <  2:1620948/150‑1 (MQ=255)

                                                                                                                                     ggggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa                 <  2:771299/150‑1 (MQ=255)

                                                                                                                                      acgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTcc                                                    <  2:453022/109‑1 (MQ=255)

                                                                                                                                       ggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTTAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACa               <  2:1977653/150‑1 (MQ=255)

                                                                                                                                       cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACa               >  2:483897/5‑150 (MQ=255)

                                                                                                                                        gTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAt              >  1:1750673/1‑150 (MQ=255)

                                                                                                                                        gTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAt              >  2:1809167/1‑150 (MQ=255)

                                                                                                                                         tAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATATGGTCTGGTTTCCATACAAAACGTCCCGCCATAGGCGGGCCGGATTTACAtt             >  1:293205/1‑150 (MQ=255)

                                                                                                                                                gCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATg      <  2:270959/150‑1 (MQ=255)

                                                                                                                                                gCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATg      <  2:694623/150‑1 (MQ=255)

                                                                                                                                                gCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATg      >  1:958287/1‑150 (MQ=255)

                                                                                                                                                  gATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCg    <  2:1984637/150‑1 (MQ=255)

TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTT  >  NC_000913/4378711‑4379008

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